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Abstract 

Hurricanes,  powerful  storms  with  wind  speeds  that  can  exceed  80  m/s,  are  one  of  the 
most  destructive  natural  disasters  known  to  man.  While  current  satellite  technology 
has  made  it  possible  to  effectively  detect  and  track  hurricanes,  expensive  ‘hurricane¬ 
hunting’  aircraft  are  required  to  accurately  classify  their  destructive  power.  Here 
we  show  that  passive  undersea  acoustic  techniques  may  provide  a  promising  tool  for 
accurately  quantifying  the  destructive  power  of  a  hurricane  and  so  may  provide  a  safe 
and  inexpensive  alternative  to  aircraft-based  techniques. 

It  is  well  known  that  the  crashing  of  wind-driven  waves  generates  underwater 
noise  in  the  10  Hz  to  10  kHz  range.  Theoretical  and  empirical  evidence  are  combined 
to  show  that  underwater  acoustic  sensing  techniques  may  be  valuable  for  measuring 
the  wind  speed  and  determining  the  destructive  power  of  a  hurricane.  This  is  done 
by  first  developing  a  model  for  the  acoustic  intensity  and  mutual  intensity  in  an 
ocean  waveguide  due  to  a  hurricane  and  then  determining  the  relationship  between 
local  wind  speed  and  underwater  acoustic  intensity.  Acoustic  measurements  of  the 
underwater  noise  generated  by  hurricane  Gert  are  correlated  with  meteorological  data 
from  reconnaissance  aircraft  and  satellites  to  show  that  underwater  noise  intensity 
between  10  and  50  Hz  is  approximately  proportional  to  the  cube  of  the  local  wind 
speed.  Prom  this  it  is  shown  that  it  should  be  feasible  to  accurately  measure  the 
local  wind  speed  and  quantify  the  destructive  power  of  a  hurricane  if  its  eye  wall 
passes  directly  over  a  single  underwater  acoustic  sensor.  The  potential  advantages 
and  disadvantages  of  the  proposed  acoustic  method  are  weighed  against  those  of 
currently  employed  techniques. 

It  has  also  long  been  known  that  hurricanes  generate  microseisms  in  the  0.1  to 
0.6  Hz  frequency  range  through  the  non-linear  interaction  of  ocean  surface  waves. 
Here  we  model  microseisms  generated  by  the  spatially  inhomogeneous  waves  of  a 
hurricane  with  the  non-linear  wave  equation  where  a  second-order  aooustic  field  is 
created  by  first-order  ocean  surface  wave  motion.  We  account  for  the  propagation  of 
microseismic  noise  through  range-dependent  waveguide  environments  from  the  deep 
ocean  to  a  receiver  on  land.  We  compare  estimates  based  on  the  ocean  surface  wave 
field  measured  in  hurricane  Bonnie  with  seismic  me8isurements  from  Florida. 

Thesis  Supervisor:  Nicholas  C.  Makris 

Title:  Associate  Professor  of  Mechanical  Engineering 
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The  wind  bloweth  where  it  listeth,  and  thou  hearest  the  sound  thereof,  but  canst 
not  tell  whence  it  cometh,  and  whither  it  goeth;  so  is  every  one  that  is  born  of  the 

Spirit. 

John  3:8  (KJV) 
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total  destructive  power  of  the  hurricane  is  proportional  to  the  cube  of 
the  maximum  wind  speed,  which  occurs  in  the  eye  wall.  [53]  .  34 

2-2  Cross  section  of  the  stratified  ocean  waveguide  showing  the  geometry 
of  the  surface  noise  problem  (Not  to  scale).  On  the  surface  is  the  area 
covered  by  the  hurricane  and  surrounding  area  covered  5  m/s  ambient 
winds.  The  surface  noise  sources  are  modeled  as  a  plane  of  monopoles 
a  small  depth  Zq  below  the  surface  and  the  sound  field  is  measured  by 
receivers  and  receiver  arrays .  37 

2-3  Sound  speed  profiles  c(z)  for  the  North  Atlantic[82]  and  the  Bay  of 
Bengal[69,  114].  The  bottom  has  a  density  of  1.38  g/cm  and  an  attenu¬ 
ation  of  0.3  dB/A  corresponding  to  the  deep  silty  sediment  layers  of  the 
Bay  of  Bengal[123,  48]  and  the  North  Atlantic  Abyssal  plain[49,  125]. 

The  water  has  a  density  of  1  g/cm  and  an  attenuation  of  6  x  10“^  dB/A.  42 
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2-4  Noise  spectral  level  (dB  re  /re/)  in  the  North  Atlantic  (a)  and  (b) 
and  the  Bay  of  Bengal  (c)  and  (d)  for  n  =  3.  (a)  and  (c)  show  the 
level  as  a  function  of  range  at  a  depth  of  200  m  for  50,  400  and  3200 
Hz  frequencies.  Lv{p)  =  10  log  ( V'(p)/ Ke/)  is  plotted  for  comparison 
where  Ke/  =  5  m/s.  Ly  =  0  is  equivalent  to  V  =  5  m/s  and  Ly  =  10 
is  equivalent  to  K  =  50  m/s.  (b)  and  (d)  show  the  level  as  a  function  of 
range  and  depth  at  50  Hz.  In  both  waveguide  environments  the  noise 
level  closely  follows  the  local  wind  speed.  In  the  North  Atlantic  there  is 
a  convergence  zone  structure  corresponding  to  sound  that  propagates 
from  the  hurricane’s  eye  wall  via  refraction.  Note  the  convergence  zone 
near  the  surface  at  a  range  of  257  km  and  the  ray  vertex  depth  of  4.7 
km . 
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2-5  Error  induced  by  the  local  area  approximation  (Eq.  (2.9))  as  a  function 
of  local  source  area  radius  Riocai  for  a  single  sensor  under  the  maximum 
winds  in  the  eye  wall  of  a  hurricane.  Curves  are  shown  for  the  North 
Atlantic  and  the  Bay  of  Bengal  environments  used  in  this  paper  as  well 
as  for  infinite  half-space  and  shallow  water  continental  shelf  environ¬ 
ments.  Plots  are  given  for  sensor  depths  of  100  m  ((a)  and  (c))  and 
800  m  ((b)  and  (d))  and  for  frequencies  of  50  Hz  ((a)  and  (b))  and  400 
Hz  ((c)  and  (d)).  While  these  plots  are  given  for  n  =  3  the  difference 
for  values  n  =  1  to  4  is  less  than  0.1  dB.  The  North  Atlantic,  Bay  of 
Bengal  and  infinite  half-space  environments  are  very  similar.  In  these 
deep-water  environments,  for  the  shallow  100  m  sensor  depth,  we  see 
that  for  Riocai  over  roughly  300  m  the  error  from  this  approximation 
is  negligible.  For  the  deeper  800  m  sensor  depth,  the  Riocai  for  which 
the  error  is  negligible  is  roughly  2  km.  In  shallow  water  the  error  in 
the  local  area  approximation  is  higher  leading  to  a  larger  Riocai-  This 
is  likely  due  to  the  strong  reflection  of  sound  off  bottom.  In  deep 
water  environments  bottom  reflections  have  little  effect  and  most  of 
the  sound  measured  by  a  receiver  propagates  via  direct  path  from  the 
surface  source . 
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2-6  Simulated  noise  spectral  level  (dB  re  Iref)  in  the  North  Atlantic  for 
range-independent  winds  as  a  function  of  wind  speed  and  frequency 
including  attenuation  by  sea-surface  bubbles  assuming  n  =  3  from  Eq. 

(2.5).  Below  100  Hz  the  power-law  relationship  between  noise  intensity 
and  wind  speed  is  unaffected  by  bubble  attenuation  even  up  to  the  80 
m/s  wind  speeds  of  a  hurricane.  As  frequency  increases,  attenuation 
affects  the  noise  level  at  progressively  lower  wind  speeds.  For  a  given 
frequency  the  noise  level  increases  linearly  with  wind  speed,  peaks,  and 
then  decays  exponentially .  51 

2-7  (a)  Noise  spectral  level  (dB  re  /re/)  as  a  function  of  wind  speed  at 

several  frequencies,  assuming  n  =  3.  The  black  curves  show  the  atten¬ 
uation,  caused  by  bubbles,  at  50  Hz,  400  Hz  and  4  kHz.  The  range 
of  wind  speeds  typical  of  a  hurricane  is  also  shown,  (b)  Noise  spectral 
level  curves  as  a  function  of  frequency  for  typical  hurricane  wind  speeds 
of  30,  50  and  80  m/s.  The  black  curves  show  the  full  areal  integration 
from  Eq.  (2.5)  and  the  grey  curves  show  the  first-order  approximation 
of  the  field  given  by  Eq.  (2.19)  with  Eqs.  (2.15),  (2.20)  and  (2.21).  .  .  53 

2-8  The  percent  RMSE  i^rmse  (a)  and  percent  bias  i/bias  (b)  of  the  wind 
speed  estimate  V  where  attenuation  by  sea-surface  bubbles  is  insignif¬ 
icant,  evaluated  numerically  from  Eqs.  (2.5)  and  (2.27).  For  time- 
bandwidth  products  p  >  5  the  estimate  becomes  unbiased  and  the 
RMSE  attains  the  Cramer-Rao  Lower  Bound.  Piggott[102]  and  Perrone[96] 
have  measured  wind  noise  level  with  standard  deviations  less  than  1 
dB  which  corresponds  to  /x  >  19.  For  //  =  19  the  percent  RMSE  in 
the  wind  speed  estimate  ranges  from  6  to  25%  depending  on  n  which 
is  a  significant  improvement  over  the  primary  satellite  classification 
method .  57 
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2-9  The  percent  RMSE  u^mse  (a)  and  percent  bias  Vbia,  (b)  of  the  wind 
speed  estimate  V  including  the  effect  of  attenuation  calculated  numer¬ 
ically  from  Ekjs.  (2.5)  and  (2.27),  assuming  n  =  3,  at  /  =  400  Hz 
where  Vmax  =  58  m/s.  The  error  and  bias  increase  for  V  fa  Vmax  but 
for  /i  >  5  and  for  values  of  V  where  vi^at  <1%  the  percent  RMSE 
decreases  and  attains  the  Cramer-Rao  Lower  Bound.  For  spectral  in¬ 
tensity  measurements  with  /i  =  19  the  percent  RMSE  in  this  example 
is  between  9  and  20%  .  58 


2-10  Angular  spectral  density  level  \Qlog{B{Q.,  f))  (dB  re  )  at  100  Hz 
for  a  64-element  A/2-spaced  horizontal  broadside  array  as  a  function 
of  steering  angle  for  hurricane  generated  noise  in  the  North  Atlantic 
at  ranges  of  257  km  (a),  289  km  (b),  and  385  km  (c)  from  the  eye  of 
the  hurricane,  assuming  n  =  3.  Ranges  of  257  and  385  km  correspond 
to  the  fourth  and  sixth  convergence  zones  from  the  center  of  the  hur¬ 
ricane.  289  km  is  exactly  between  the  fourth  and  fifth  convergence 
zones.  Curves  are  shown  for  a  powerful  72  m/s  hurricajie,  a  medium 
50  m/s  hurricane  and  a  weak  33  m/s  hurricane.  The  angular  spectral 
density  level  from  ambient  noise  is  plotted  for  comparison.  A  steering 
angle  of  0°  corresponds  to  the  array  steered  toward  the  calm  eye  of  the 
hurricane  and  the  powerful  eye  wall  is  located  at  ^3°.  This  array  has 
an  angulaj  resolution  of  1.8°  which  at  a  range  of  257  km  corresponds 
to  a  8  km  spatial  resolution .  64 
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3-1  Spectrogram  of  the  hurricane  noise  in  dB  re  /iPa^/Hz,  from  1  to  50 
Hz,  received  by  the  NOAA  hydrophone  on  15  Sept.  Above  10  Hz  the 
maximum  between  13:00  and  16:00  correspond  to  the  powerful  winds 
of  the  hurricane’s  eye-wall.  At  1  Hz  we  see  noise  most  likely  due  to 
cable  strum.  There  seems  to  be  some  hurricane  related  noise  between 
1  and  10  Hz,  however  it  does  not  correlate  well  with  wind  speed  and 
may  be  caused  by  nonlinear  surface  wave  interaction  as  discussed  in 
the  next  chapter .  70 

3-2  Log  of  the  noise  spectral  intensity  I  (solid  blue)  in  Watts/m^Hz,  from 
10  to  50  Hz,  received  by  the  NOAA  hydrophone  on  15  Sept.  The  ex¬ 
pected  intensity  (red  dashed  line)  is  based  on  the  hindcast  wind  speed 
V  in  m/s  above  the  hydrophone  given  the  best  fit  relation  log[//(lWatt/m^Hz)] 
nlog[V"/(lm/s)]  -I-  h  where  n  =  3.36  and  b  =  —14.46  .  The  maxima  at 
13:30  and  15:30  correspond  to  the  powerful  winds  of  the  hurricane’s 
eye-wall  and  the  minimum  at  14:30  corresponds  to  the  hurricane’s  eye. 

The  percent  root-mean-square  error  between  the  two  curves  is  5%.  .  .  71 

3-3  Wind  speed  in  m/s  in  hurricane  Gert  as  a  function  of  position  relative 
to  the  center  or  ‘eye’.  As  hurricane  Gert  passed  over  the  acoustic 
sensor  the  sensor  effectively  sliced  a  path  through  the  hurricane.  The 
lines  show  the  paths  estimated  based  on  Satellite  (dash-dot).  Aircraft 
(dashed)  and  Acoustic  best  fit  (solid)  respectively .  74 

3-4  The  log  of  the  noise  spectral  intensity  I  in  Watts/m^Hz  versus  the  log 
of  the  surface  wind  speed  T  in  m/s  at  the  NOAA  hydrophone  (circles) 
based  on  the  best  fit  hurricane  track.  The  best  fit  linear  regression 
shows  a  log[//(lWatt/m^Hz)]  =  3.361og[V’/(lm/s)]  — 14.46  relationship 
between  intensity  and  wind  speed  with  a  percent  root-mean-square 
error  of  only  5% .  76 
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3-5  The  upper  plot  gives  the  root-mean-square  error  (RMSE)  of  the  lin¬ 
ear  regression  for  each  possible  hurricane  track.  The  y-axis  is  the 
north-south  distance  between  the  hydrophone  and  the  possible  hur¬ 
ricane  track  and  the  x-axis  is  the  speed  at  which  the  hurricane  was 
moving.  The  minimum  RMSE  of  0.76  dB,  indicated  y  the  -I-  symbol,  is 
taken  to  represent  the  ‘best  fit’  hurricane  track.  Note  that  the  RMSE 
is  small,  less  than  0.8  dB,  over  a  range  of  possible  tracks.  The  lower 
plot  gives  the  slope  n  of  the  linear  regression  for  each  possible  hurri¬ 
cane  track  with  the  same  x  and  y  axes  as  the  upper  plot.  For  the  best 
track  the  slope  n  =  3.36.  For  the  r^lnge  of  tracks  with  a  RMSE  less 
than  0.8  dB,  the  slope  of  the  regression  could  veiry  from  3  to  3.6.  ...  78 


3-6  Log  of  noise  intensity  I  versus  log  of  hindcast  local  wind  speed  V 
based  on  the  best-fit  track  (red  circles).  Also  plotted  is  the  log  of 
noise  intensity  1  versus  log  of  local  wind  speed  V  measured  by  Cato 
and  Tavener  [23]  (blue  circles)  at  lower  non-hurricane  wind  speeds 
off  Australia.  The  linear  regression  shows  a  log(//l  Watt/m^Hz)  = 
3.341og(F/l  m/s)  —  14.5  relationship  between  noise  intensity  I  and 
wind  speed  V .  79 
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4-1  (A)  Wind  speed  in  m/s  and  (B)  microseismic  source  level  (w,  -Zj  = 

0)  at  cj  =  1.26  rad/s  (/  =  0.2  Hz)  in  dB  re  Pa^s/m^  from  Eq.  (4.39) 
at  1200  on  24  August  1998  as  a  function  of  latitude  and  longitude. 

The  arrow  indicates  the  direction  hurricane  Bonnie  was  moving.  The 
letters  (a),  (b),  (c)  and  (d)  represent  features  of  interest;  (a)  indicates 
the  location  of  maocimum  wind  speed,  (b)  indicates  the  eye  where  wind 
speed  is  zero,  and  (c)  and  (d)  indicate  peaks  in  the  microseismic  source 
level.  This  figure  shows  that,  while  a  hurricjine  can  produce  significant 
microseismic  source  levels  (B),  these  level  do  not  directly  follow  the 
wind  speeds  (A)  in  the  hurricane .  87 


4-2  The  wave-height  power  spectral  level  (10  log Smm{K’) / (^^^Z {radP /rn^))) 
in  dB  re  m^/(rad^/m^)  at  the  locations  of  interest  (a),  (b),  (c),  and  (d) 
given  in  Fig.  4.2  at  1200  on  24  August  1998.  The  peak  in  the  spectra  is 
at  roughly  k  =  0.04  rad/m  or  cr  =  0.63  rad/s  which  corresponds  to  an 
acoustic  frequency  of  w  =  1.26  rad/s  or  /  =  0.2  Hz.  At  some  locations 
(b,  c,  and  d)  there  are  waves  propagating  with  opposing  wavenumber 
vectors  k,  while  at  other  locations  (a),  most  of  the  waves  propagate  in 
the  same  direction.  From  Eq.  (4.37)  we  expect  these  locations  with 
opposing  waves  to  produce  the  greatest  microseismic  source  levels  and 
shown  in  Fig.  4.2(B) .  88 
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4-3  Geometry  of  the  hurricane  wave  field  and  ocean  waveguide  (not  to 
scale).  In  our  model  the  waveguide  environment  may  be  range  depen¬ 
dent  and  in  this  paper  we  consider  the  example  of  an  upslope  propa¬ 
gation  from  the  deep  North  Atlantic  to  Florida.  The  range  and  ocean 
depth  parameters  R  and  d  are  given  in  Table  4.6.  The  depth  of  the 
receiver  below  the  ocean  bottom  of  162  m  corresponds  to  the  depth 
below  the  earth  surface  of  the  actual  seismometer  in  Florida.  The 
compression  wave  speeds  Cp  are  1500  and  5200  m/s  in  the  water  and 
bottom  respectively.  The  shear  wave  speed  in  the  bottom  is  3000 
m/s.  The  densities  p  in  the  water  and  bottom  are  1.0  and  2.5  g/cm 
respectively .  98 


4-4  The  ocean  depth  (solid  line)  between  hurricane  Bonnie  and  the  seis¬ 
mometer  in  Florida  at  noon  on  24  August.  Also  shown  is  the  ocean 
depth  for  the  idealized  up-sloping  environment  used  to  calculate  the 
Green  functions  in  Section  4.6  (dashed  line).  The  scale  of  the  figure 
makes  the  actual  slope  appear  to  change  rapidly,  however,  the  maxi¬ 
mum  slope  of  the  ocean  floor  is  roughly  3° .  110 
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4-5  Horizontal  (a-d)  and  vertical  (e-h)  velocity  spectra,  as  a  function  of  seis¬ 
mic  frequency  (/  =  27ra;),  modeled  (red)  and  measured  (blue)  at  the 
seismometer  in  Florida  at  four  times  on  24  August  1998.  The  dashed 
blue  lines  represent  portions  of  the  data  that  where  corrupted  by  non¬ 
hurricane  related  ambient  noise.  Note  that  the  peak  in  the  spectra  is 
at  roughly  /  =  0.2  Hz.  This  frequency  corresponds  to  the  peak  in  the 
wave-height  power  spectral  level  at  /c  =  0.04  rad/m  (surface  wave  fre¬ 
quency  /  =  27ra  =  0.1  Hz)  seen  in  Fig.  4.2  and  in  plots  (i-1).  The  sur¬ 
face  wave-height  power  spectral  level  (101ogS’,„m(/c)/(m^/(rad^/Tn^))) 
in  plots  (i-1)  is  taken  at  the  center  of  the  storm  and  averaged  over  wave 
propagation  direction .  Ill 

G-1  The  RHS  and  LHS  in  log  (upper  plot)  and  linear  (lower  plot)  domains 
from  13.33  to  13.4  seconds.  We  see  that  the  peaks  in  the  curves  occur  at 
the  same  times  but  that  the  amplitudes  of  the  peaks  differ  significantly. 

We  also  see  that  the  RHS  peaks  are  higher  for  the  earlier  low-order 
modes  while  the  LHS  peaks  are  higher  for  the  later  high-order  modes. 

Note  that  the  RHS  in  the  lower  plot  has  been  scaled  by  a  factor  of  five 
so  that  it  can  be  compared  against  the  LHS .  152 

G-2  The  RHS  and  LHS  (upper  plot)  over  a  longer  time  duration.  We 
continue  to  see  the  general  trend  where  the  peaks  of  the  RHS  are  higher 
for  low-order  early  arrivals  and  lower  for  later  high-order  arrivals.  This 
can  also  be  seen  in  the  lower  plot  which  shows  the  difference  between 
the  RHS  and  LHS .  153 

G-3  The  weighting  term  as  a  function  of  mode  number  m  at  several 

frequencies  in  an  isovelocity  waveguide.  The  source  depth  <<  A 
typical  of  ocean  surface  noise .  154 
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Chapter  1 


Introduction 


In  this  thesis  we  theoretically  and  experimentally  study  the  noise  generated  by  a  hur¬ 
ricane  and  show  that  passive  undersea  acoustic  techniques  may  provide  a  promising 
tool  for  accurately  quantifying  the  destructive  power  of  a  hurricane  and  may  provide  a 
safe  and  inexpensive  alternative  to  aircraft-based  techniques.  While  current  satellite 
technology  has  made  it  possible  to  effectively  detect  and  track  hurricanes,  expensive 
‘hurricane-hunting’  aircraft  are  required  to  accurately  classify  its  destructive  power. 

Quantification  of  a  hurricane’s  total  destructive  power,  which  is  proportional  to  the 
cube  of  the  hurricane’s  maximum  wind  speed  [52],  is  critical  for  hurricane  planning. 
The  destructive  power  of  a  tropical  cyclone  was  recently  demonstrated  by  hurricane 
Katrina  which  caused  over  1000  fatalities[29]  and  an  estimated  economic  loss  of  100 
billion  dollars[10]  and  in  1970  a  hurricane  killed  over  300,000  people  in  Bangladesh 
in [33].  Prior  to  Katrina  the  United  States  Commission  on  Ocean  Policy  emphasized 
the  need  for  more  accurate  quantifications  of  hurricane  destructive  power  to  improve 
disaster  planning[133].  Ingiccurate  quantifications  can  lead  to  poor  forecasting  and 
unnecessary  evacuations,  which  are  costly,  or  missed  evacuations,  which  can  result  in 
loss  of  life[36].  These  fatalities  and  costs  can  be  reduced  if  the  public  is  given  timely 
and  accurate  advanced  warning,  but  this  depends  on  the  ability  to  accurately  quantify 
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hurricane  wind  speed  while  it  is  still  far  from  land.  Current  classification  and  warning 
systems  save  an  average  $2.5  billion  each  year  in  the  United  States  alone[133]  and 
more  accurate  systems  could  save  even  more. 


1.1  Current  Hurricane  Classification  Techniques 

Currently  there  are  two  primary  tools  used  by  meteorologists  for  detecting  and  classi¬ 
fying  hurricanes  in  the  ocean;  satellites  and  aircraft.  Satellites  are  useful  for  detecting 
and  locating  hurricanes  but  their  poor  spatial  resolution  in  measuring  wind  speeds 
limits  their  ability  to  accurately  classify  the  hurricane.  Aircraft  are  much  better  for 
classifying  hurricanes,  however,  their  expense  means  that  their  use  is  not  feasible  for 
most  countries[33].  Also  their  limited  endurance  means  that  they  cannot  monitor  the 
hurricane  continuously. 

The  Dvorak  method  [31,  32,  131]  is  the  most  common  method  for  classifying  a 
hurricane’s  destructive  power.  This  method,  where  hurricane  cloud  features  observed 
in  satellite  images  are  interpreted  to  estimate  wind  speed  and  classify  destructive 
power,  can  yield  errors  in  wind  speed  estimates  exceeding  40%  [93,  43,  12,  4,  121] 
when  compared  to  the  best-estimate  wind  speed  from  aircraft  meeisurements.  For 
example,  of  the  eight  North  Atlantic  hurricanes  of  2000,  three  of  them  [93,  43,  12] 
experienced  Dvorak  errors  over  40%  and  three  more  [94,  72,  120]  experienced  Dvorak 
errors  over  20%.  Despite  these  errors,  the  Dvorak  method  is  still  the  primary  tech¬ 
nique  for  classifying  the  destructive  power  of  a  hurricane  from  satellite  measurements 
[44].  A  satellite-based  pattern-recognition  technique  similar  to  the  Dvorak  method 
using  SSM/I  satellite  microwave  (85  GHz)  instead  of  optical  and  infrared  images  has 
recently  been  developed  but  gives  similar  errors  as  the  Dvorak  method  [7]. 

Satellite  classification  of  hurricanes  with  Microwave  Sounding  Units  (MSU)  [130] 
is  secondary  to  the  primary  Dvorak  method  [44]  due  to  the  limited  spatial  resolution 
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of  the  unit.  The  55  GHz  microwave  radiation  given  off  by  warm  air  in  the  hurricane’s 
eye  is  used  to  estimate  temperature  and  then  infer  the  hurricane’s  power.  Because 
of  the  small  size  of  the  satellite  array  unit  the  spatial  resolution  of  the  measurement 
is  about  48  km  [65],  which  is  often  larger  than  the  diameter  of  the  eye,  resulting  in 
a  blurred  image  of  the  hurricane  and  potentially  leading  to  errors  in  the  estimate  of 
the  destructive  power  [130,  65]. 


Other  satellite  techniques  for  estimating  hurricane  wind  speed  and  destructive 
power  are  currently  being  studied.  For  an  overview  see  the  article  by  Katsaros, 
Vachon,  Liu  and  Black  [59].  These  techniques,  however,  are  still  under  development 
and  are  not  yet  being  used  operationally  for  hurricane  classification  and  disaster 
planning  [41]. 


More  accurate  classification  can  be  achieved  by  flying  specitdized  ‘hurricane-hunting’ 
aircraft,  like  the  Air  Force’s  WC-130  and  NOAA’s  WP-3,  through  the  high  winds  of 
a  hurricane  [41].  Using  on-board  sensors  and  expendable  dropsondes,  accurate  wind 
speed  estimates  with  errors  less  than  5  m/s  can  be  obtained  [41].  Unfortunately 
the  expense  of  these  aircraft  prohibits  their  routine  use  outside  of  the  United  States 
[33].  For  example,  the  cost  to  purchase  a  WC-130  aircraft  is  roughly  $82  million  [2] 
(adjusted  for  inflation  to  year  2005  dollars)  and  the  deployment  cost  is  $155,000  per 
flight  [9].  Given  these  limitations  in  current  capability,  the  United  States  Commission 
on  Ocean  Policy  has  recommended  that  future  ocean  observing  systems  be  used  to 
improve  weather  related  warnings  [133].  We  experimentally  and  theoretically  demon¬ 
strate  that  underwater  acoustic  measurements  of  noise  intensity  may  provide  a  useful 
technique  for  hurricane  wind  speed  estimation. 
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1.2  Acoustic  Field  from  Spatially  Uncorrelated  Sur¬ 
face  Noise  Sources 

Sea-surface  agitation  from  the  action  of  wind  and  waves  is  a  dominant  source  of 
ambient  noise  in  the  ocean  [67,  136]  in  the  10  Hz  to  10  kHz  frequency  range.  We  show 
that  it  may  be  practical  to  inexpensively  determine  local  wind  speed  and  quantify  the 
destructive  power  of  a  hurricane  by  measuring  this  noise  which  can  be  described  as  a 
sum  of  fields  radiated  from  many  random  sources  on  the  sea  surface  [26,  76,  70,  19, 
95,  81].  If  the  surface  noise  sources  have  the  same  statistical  distribution,  Ingenito 
and  Wolf  have  shown  that  the  wind-generated  noise  spectral  intensity  is  the  product 
of  two  separate  factors,  a  waveguide  propagation  factor  and  a  “universal  ambient 
noise”  [57]  source  factor  which  is  a  function  of  wind  speed  but  otherwise  is  expected 
to  be  effectively  independent  of  horizontal  position. 

Shaw,  Watts  and  Rossby[118]  first  considered  the  concept  of  using  underwater 
sound  to  estimate  wind  speed  for  spatially  uniform  wind  speed  distributions.  They 
found  sound  pressure  level  in  dB  to  be  linearly  related  to  the  log  of  the  wind  speed 
and  the  idea  of  a  universal  ambient  noise  source  factor  was  implicit  in  their  approach. 
We  will  show  that  the  slope  of  their  linear  relationship  corresponds  to  the  universal 
ambient  noise  factor  and  the  intercept  to  the  waveguide  propagation  factor.  Evans, 
Watts,  Halpern  and  Bourassa  demonstrate  that  these  wind  speed  estimates  could  be 
made  to  within  ±1  m/s  at  low  5  to  10  m/s  wind  speeds  [38],  which  is  much  less  than 
is  experienced  in  a  hurricane. 

Many  experiments  have  been  conducted  to  determine  the  relationship  between 
local  wind  speed  and  underwater  noise  intensity  £is  noted  in  Ref.  [18].  A  common 
difficulty  in  these  experiments  has  been  contamination  from  shipping  noise  [18,  61]. 
This  typically  leads  to  poorer  correlation  and  greater  variance  in  estimates  of  the 
relationship  between  wind  speed  and  noise  intensity  [18].  Two  experimental  stud- 


26 


ies  conducted  over  many  months  that  minimized  this  contamination  show  that  a 
consistent  high-correlation  power-law  relationship  exists  [102,  38].  They  also  show 
underwater  noise  intensity  to  be  linearly  proportional  to  wind  speed  to  a  frequency 
dependent  power,  ranging  from  two  to  four,  for  wind  speeds  between  5  and  20  m/s. 
It  is  possible  that  at  higher  frequencies  attenuation  by  bubbles  could  cause  a  roll-off 
in  the  power-law  relationship  [139].  This  attenuation,  however,  is  insignificant  at  low 
frequencies  and  can  be  accurately  measured  and  modeled  at  high  frequencies. 

In  this  thesis  we  present  the  first  published  data  relating  ambient  noise  and  wind 
speed  in  a  hurricane.  In  1999  an  autonomous  underwater  acoustic  sensor  package  in 
the  North  Atlantic  recorded  the  underwater  noise  from  crashing  wind-driven  waves  as 
hurricane  Gert  passed  overhead.  By  correlating  this  noise  with  meteorological  data 
from  reconnaissance  aircraft  and  satellites  we  show  that  underwater  noise  intensity 
between  10  and  50  Hz  is  approximately  proportional  to  the  cube  of  the  local  wind 
speed.  Passive  underwater  acoustic  intensity  measurements  from  a  single  sensor  may 
then  be  used  to  estimate  hurricane  wind  speed  to  within  a  5%  error  and  from  this 
accurately  quantify  the  destructive  power  of  the  hurricane. 

Based  on  theoretical  and  empirical  evidence  we  find  that  it  may  be  possible  to 
estimate  local  hurricane  wind  speed  by  generalizing  the  appro8ich  of  Shaw,  Watts  and 
Ross  by  [118].  We  show  that  the  wind-generated  noise  received  by  a  single  underwater 
acoustic  sensor  in  a  hurricane  can  be  well  approximated  by  sea-surface  contributions 
so  local  that  wind  speed  and  surface  source  intensity  can  be  taken  as  nearly  constant. 
With  these  findings,  noise  intensity  can  be  well  approximated  zis  the  product  of  a 
local  universal  ambient  noise  source  factor  and  a  waveguide  propagation  factor  even 
for  the  range-dependent  wind  speeds  of  a  hurricane. 

At  low  frequencies,  below  roughly  100  Hz,  we  show  that  attenuation  by  wind- 
induced  bubbles  in  the  upper-ocean  boundary  layer  should  be  insignificant  even  in 
hurriczme  conditions.  Temporal  variations  in  underwater  noise  intensity  should  then 
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be  primarily  caused  by  the  universal  ambient  noise  source  factor  which  is  expected 
to  depend  on  local  wind  speed  and  will  vary  as  a  hurricane  advects  over  a  fixed 
receiver.  By  extrapolating  known  relationships  [102]  between  wind  speed  and  noise 
level  in  this  frequency  range,  the  ambient  noise  level  should  increase  monotonically 
with  wind  speed,  and  it  should  be  possible  to  directly  estimate  local  wind  speed  from 
measured  noise  level. 

At  higher  frequencies  temporal  variations  in  underwater  noise  intensity  may  also 
be  caused  by  attenuation  due  to  scattering  from  bubbles  in  the  upper-ocean  boundary 
layer.  This  attenuation  increases  with  wind  speed  and  acoustic  frequency.  Farmer 
and  Lemon  [40]  experimentally  show  that  this  leads  to  a  frequency  dependent  peak 
in  noise  level  versus  wind  speed  at  frequencies  above  8  kHz  and  wind  speeds  above 
15  m/s.  We  analytically  show  that  such  a  peak  may  also  exist  for  frequencies  above 
100  Hz  in  typical  hurricane  wind  speeds.  Since  the  shape  of  the  ambient  noise  versus 
wind  speed  curve  and  the  location  of  its  peak  vary  strongly  with  frequency,  we  show 
that  wind  speed  may  still  be  unambiguously  estimated  from  broadband  ambient  noise 
measurements  in  hurricane  conditions  above  100  Hz  once  the  corresponding  universal 
source  dependence  is  empirically  determined. 

The  accuracy  of  underwater  acoustic  wind  speed  estimates  depends  on  the  signal 
to  noise  ratio  (SNR)  of  the  underwater  ambient  noise  intensity  measurements  upon 
which  they  cire  based.  Piggott  [102]  and  Perrone  [97]  have  consistently  measured 
wind  noise  with  a  standard  deviation  of  less  than  one  dB  and  in  measurements  of 
hurricane  Gert  we  find  standard  deviations  of  0.7  dB,  as  expected  from  theory  where 
the  variance  of  the  intensity  measurement  can  be  reduced  by  stationary  averaging 
[99,  80,  81].  For  previously  measured  power-law  relationships  that  range  from  quartic 
to  square  [102,  38],  a  one  dB  standard  deviation  in  sound  pressure  level  corresponds 
to  a  6%  to  12%  respective  error  in  estimated  wind  speed.  The  3.3  power  law  and 
0.7  dB  standard  deviation  measured  in  hurricane  Gert  results  in  only  a  5%  error  in 
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estimated  wind  speed. 

Ocean  acoustics  then  has  serious  potential  for  providing  accurate  and  inexpensive 
hurricane  classification  estimates.  Since  a  single  hydrophone  effectively  measures 
only  the  local  surface  noise,  it  will  effectively  cut  a  swath  through  the  hurricane 
yielding  local  wind  speed  estimates  as  the  storm  passes  over.  At  low  frequencies, 
current  evidence  suggests  a  simple  power-law  relationship  between  noise  intensity 
and  wind  speed.  At  higher  frequencies,  a  frequency-dependent  roll-off  is  expected 
in  the  relationship  due  to  attenuation  by  bubbles.  Wind  speed  can  still  be  uniquely 
estimated,  however,  by  making  broad-band  measurements  at  higher  frequency. 

In  Chapter  2  we  review  models  for  range-dependent  noise  in  the  ocean  and  develop 
a  model  for  wind  generated  noise  from  a  hurricane  for  both  single  sensors  and  arrays. 
We  use  this  model  to  demonstrate  the  potential  usefulness  of  quantifying  hurricane 
wind  speed  with  underwater  acoustic  sensors.  We  also  review  past  experiments  where 
the  relationship  between  underwater  noise  intensity  and  wind  speed  was  measured. 

In  Chapter  3  we  present  data  from  an  autonomous  underwater  acoustic  sensor 
package  in  the  North  Atlantic  that  recorded  the  underwater  noise  from  crashing  wind- 
driven  waves  as  hurricane  Gert  passed  overhead.  This  data  further  demonstrates  the 
potential  usefulness  of  classifying  hurricanes  with  underwater  acoustic  sensors. 


1.3  Acoustic  Field  from  the  Interaction  of  Surface 
Gravity  Waves  on  the  Ocean  Surface 

Hurricanes  also  generate  seismic  noise,  commonly  referred  to  as  microseisms,  in  the 
0.1  to  0.6  Hz  frequency  range.  In  Chapter  4  we  describe  the  microseisms  generated 
by  the  spatially  inhomogeneous  waves  in  a  hurricane.  Using  the  ocean  surface  di¬ 
rectional  wave  spectrum  in  hurricane  Bonnie  [142,  88]  we  hindcast  the  microseismic 
field  and  compare  it  with  seismic  measurements  from  Florida.  Previously  hurricane 
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surface  directional  wave  spectra  had  not  been  adequately  measured  or  modeled  so 
that  researchers  had  to  rely  on  assumed  spectra  in  order  to  model  hurricane  micro¬ 
seisms  [78].  This  analysis  is  useful  because  microseisms  are  a  primary  cause  of  noise 
in  seismic  measurements  [74,  107,  105]  that  raise  the  detection  threshold  for  moni¬ 
toring  earthquakes  [135]  and  tsunamis.  Historically  microseismic  measurements  have 
also  been  used  to  track  hurricanes  [107,  45]  although  that  task  is  now  accomplished 
using  satellites.  We  find  that,  due  to  the  nonlinear  nature  of  microseism  generation, 
measurements  of  the  microseismic  field  will  probably  not  be  useful  for  hurricane  wind 
speed  estimation. 

We  present  an  analytic  expression  for  microseism  generation  by  the  spatially  inho¬ 
mogeneous  waves  typical  in  a  hurricane  based  on  the  non-linear  wave  equation  where 
a  second-order  seismo-acoustic  field  is  generated  by  a  source  distribution  which  de¬ 
pends  on  the  first-order  ocean  surface  wave  motion.  The  seismo-acoustic  field  at  a 
receiver  can  then  be  expressed  as  the  integral  over  the  source  distribution  multiplied 
by  the  waveguide  Green  function.  This  expression  is  ideal  for  hurricane  generated 
microseisms  since  it  can  be  used  to  calculate  the  acoustic  field  due  to  spatially  in¬ 
homogeneous  surface  waves.  Also,  this  expression  may  be  used  in  range-dependent 
waveguide  environments  as  is  the  case  when  a  hurricane  at  sea  generates  microseisms 
that  propagate  up  the  continental  margin  to  a  receiver  on  land. 

Based  on  the  wave-height  spectra  in  hurricane  Bonnie,  we  calculate  the  microseis¬ 
mic  source  levels  generated  by  the  nonlinear  interaction  of  the  ocean  surface  waves. 
Our  derivation  shows  that  microseisms  are  generated  by  the  non-linear  interaction 
of  ocean  surface  waves  with  roughly  the  same  wavelength  but  nearly  opposing  prop¬ 
agation  directions.  This  is  in  agreement  with  earlier  works  [86,  78].  Recent  mea¬ 
surements  [142,  132]  and  models  [88]  of  surface  directional  wave  height  spectra  in 
hurricane  Bonnie  show  complex  patterns  with  the  opposing  surface  waves  necessary 
to  generate  microseisms. 
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Chapter  2 


Ocean  Acoustic  Hurricane  Wind 
Speed  Quantification  from 
Spatially  Uncorrelated  Surface 
Noise 


2.1  Introduction 

A  case  is  made  that  it  may  be  practical  to  inexpensively  determine  local  wind 
speed  and  quantify  the  destructive  power  of  a  hurricane  by  measuring  its  under¬ 
water  acoustic  noise  intensity.  The  dominant  source  of  ambient  noise  in  the  ocean  is 
sea-surface  agitation  from  the  action  of  wind  and  waves  [67,  136].  This  noise  can  be 
described  as  a  sum  of  fields  radiated  from  many  random  sources  on  the  sea  surface 
[26,  76,  70,  19,  95,  81].  If  the  surface  noise  sources  have  the  same  statistical  distri¬ 
bution,  Ingenito  and  Wolf  have  shown  that  wind-generated  noise  spectral  intensity  is 
the  product  of  two  separate  factors,  a  waveguide  propagation  factor  and  a  “universal 
ambient  noise”  [57]  source  factor  which  is  a  function  of  wind  speed  but  otherwise  is 
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expected  to  be  effectively  independent  of  horizontal  position. 


Shaw,  Watts  and  Rossby  [118]  found  sound  pressure  level  in  dB  to  be  linearly 
related  to  the  log  of  the  wind  speed  and  from  this  developed  the  concept  of  using 
underwater  sound  to  estimate  wind  speed  for  spatially  uniform  wind  speed  distribu¬ 
tions.  We  show  that  the  slope  of  their  linear  relationship  corresponds  to  the  universal 
ambient  noise  factor  described  by  Ingenito  and  Wolf  [57]  and  the  intercept  to  the 
waveguide  propagation  factor.  As  discussed  in  Chapter  1,  it  was  later  demonstrated 
that  these  estimates  could  be  made  to  within  ±1  m/s  in  the  5  to  10  m/s  wind  speed 
range  [38]. 


By  generalizing  the  approach  of  Shaw,  Watts  and  Rossby  [118]  we  find  that  it  may 
be  possible  to  estimate  local  hurricane  wind  speed.  The  wind-generated  noise  received 
by  a  single  underwater  acoustic  sensor  in  a  hurricane  can  be  well  approximated  by 
sea-surface  contributions  so  local  that  wind  speed  and  surface  source  intensity  can 
be  taken  as  nearly  constant.  Noise  intensity  can  then  be  well  approximated  as  the 
product  of  a  local  universal  ambient  noise  source  factor  and  a  waveguide  propagation 
factor  even  for  the  range-dependent  wind  speeds  of  a  hurricane. 


In  this  chapter  we  review  models  for  the  spatial  wind  speed  dependence  of  a 
hurricane  that  will  be  used  to  model  ambient  noise.  We  also  review  past  experiments 
that  measured  the  relationship  between  underwater  noise  intensity  and  wind  speed. 
We  review  models  for  range-dependent  noise  in  the  ocean  and  develop  a  model  for 
wind  generated  noise  from  a  hurricane  for  both  single  sensors  and  arrays.  We  use  this 
model  to  demonstrate  the  potential  usefulness  of  quantifying  hurricane  wind  speed 
with  underwater  acoustic  sensors. 
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2.2  Hurricane  Structure  and  Current  Classifica¬ 


tion  Techniques 

Hurricanes  are  severe  storms  characterized  by  surface  winds  from  33  to  over  80  m/s 
[33]  that  circulate  around  a  central  low  pressure  zone  called  the  eye.  Holland[52]  gives 
an  analytic  model  for  the  surface  wind  speed  profile  as  a  function  of  range  from  the 
eye  since  hurricanes  are  typically  cylindrically  symmetric, 


V  = 


1 


AB(p„  -  Pc) 


exp{7#} 

PaT^ 


(2.1) 


where  V  is  wind  speed  at  a  height  of  10  m  above  the  sea  surface,  Pc  and  Pn  are 
the  atmospheric  pressure  in  the  eye  and  outside  the  hurricane  respectively,  Pa  is  the 
density  of  the  air,  and  A  and  B  are  empirical  values.  Using  this  model,  the  surface 
wind  speed  profile  for  a  moderate  hurricane  is  given  in  Fig.  2-1,  where  wind  speed  in 
the  eye  is  zero  and  rapidly  increases  to  a  maximum  of  50  m/s  at  what  is  known  as 
the  eye  wall.  Outside  of  the  eye  wall,  which  is  on  the  order  of  ten  kilometers  thick, 
wind  speed  slowly  decreases  to  the  edge  of  the  hurricane  which  is  typically  hundreds 
of  kilometers  from  the  eye.  Most  of  a  hurricane’s  destructive  power  then  comes  from 
the  high  winds  in  the  eye  wall  since  this  power  is  roughly  proportional  to  the  cube  of 
the  maximum  wind  speed. [53] 

The  standard  approach  for  classifying  a  hurricane’s  destructive  power,  the  Dvo¬ 
rak  method, [31,  32,  131]  is  effectively  a  pattern-recognition  technique  where  satellite 
images,  in  the  visible  and  infra-red  spectrum,  are  used  to  classify  the  hurricane  based 
on  features  like  the  size  and  the  geometry  of  cloud  patterns.  As  discussed  in  Chapter 
1,  this  method  often  yielded  wind  speed  estimates  with  errors  of  over  40%  in  sev¬ 
eral  recent  hurricanes.  [93,  43,  12,  4,  121]  Despite  these  errors,  the  Dvorak  method  is 
still  the  primary  technique  for  classifying  the  destructive  power  of  a  hurricane  from 
satellite  measurements.  [44]  A  satellite-based  pattern- recognition  technique  similar  to 


33 


Figure  2-1:  Hurricane  wind  speed  as  a  function  of  distance  from  the  hurricane  center 
based  on  Holland’s  model[52]  with  parameters  A  =  72.44,  B  =  1.86,  Pc  —  96300  Pa, 
Pn  =  100500  Pa  and  Pajr  =  1-15  kg/m^.  The  zero  wind  speed  region  at  the  center  of 
the  hurricane  (0  km)  is  called  the  eye  and  the  high  wind  speed  region  (10  km)  is  the 
eye  wall.  The  total  destructive  power  of  the  hurricane  is  proportional  to  the  cube  of 
the  maximum  wind  speed,  which  occurs  in  the  eye  wall. [53] 
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the  Dvorak  method  using  SSM/I  satellite  microwave  (85  GHz)  instead  of  optical  and 
infrared  images  has  recently  been  developed  but  gives  similar  errors  as  the  Dvorak 
method.  [7] 


Satellite  classification  of  hurricanes  with  Microwave  Sounding  Units  (MSU)[130] 
is  secondary  to  the  primary  Dvorak  method[44]  due  to  the  limited  spatial  resolution 
of  the  unit.  The  55  GHz  microwave  radiation  given  off  by  warm  air  in  the  hurricane’s 
eye  is  used  to  estimate  temperature  and  then  infer  the  hurricane’s  power.  Because 
of  the  small  size  of  the  satellite  array  unit  the  spatial  resolution  of  the  measurement 
is  about  48  km, [65]  which  is  often  larger  than  the  diameter  of  the  eye,  resulting  in 
a  blurred  image  of  the  hurricane  and  potentially  leading  to  errors  in  the  estimate  of 
the  destructive  power.  [130,  65] 


Other  satellite  techniques  for  estimating  hurricane  wind  speed  and  destructive 
power  are  currently  being  studied.  For  an  overview  see  the  article  by  Katsaros, 
Vachon,  Liu  and  Black.  [59]  These  techniques,  however,  are  still  under  development 
and  are  not  yet  being  used  operationally  for  hurricane  classification  and  disaster 
planning. [41] 


To  overcome  the  limitations  of  satellite  techniques,  specially  equipped  aircraft, 
like  the  Air  Force’s  WC-130s  and  NOAA’s  WP-3s,  are  flown  through  the  center  of 
a  hurricane.  [41]  Using  on-board  sensors  and  expendable  dropsondes,  accurate  wind 
speed  estimates  with  errors  less  than  5  m/s  can  be  obtained. [41]  Unfortunately,  as 
discussed  in  the  introduction,  these  aircraft  are  expensive  to  purchaise  and  operate 
and  are  currently  only  used  by  the  United  States. [33] 
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2.3  Wind  Generated  Surface  Noise  from  Uncorre¬ 
lated  Surface  Sources 

Here  we  develop  a  model  for  the  surface  generated  noise  intensity  and  mutual  intensity 
from  a  hurricane  received  by  a  hydrophone  or  hydrophone  array  submerged  in  an 
ocean  waveguide.  The  geometry  of  the  problem  is  shown  in  Fig.  2-2.  The  hurricane 
is  centered  at  the  origin  and  is  surrounded  by  ambient  winds  all  of  which  cause  local 
sea-surface  agitation.  This  agitation  leads  to  sound  sources  with  amplitude  dependent 
on  the  local  wind  speed  that  are  modeled  as  a  sheet  of  monopoles  on  a  source  plane 
at  a  depth  zq  within  a  quarter  wavelength  of  the  free  surface  following  ocean  acoustic 
noise  modeling  convention. [70,  19,  95]  Intensity  and  mutual  intensity  are  determined 
by  directly  integrating  the  surface  source  contributions  using  the  waveguide  Green 
function. 

Several  previous  authors  have  addressed  similar  surface  noise  problems,  however, 
their  derivations  are  intertwined  with  approximations  or  parameterizations  that  are 
not  suitable  for  modeling  hurricane  noise.  Kuperman  and  Ingenito[70]  developed  a 
widely  used  surface  noise  model,  however,  embedded  in  their  derivation  is  the  as¬ 
sumption  that  the  source  field  is  range  independent.  This  is  not  true  for  hurricane 
generated  noise  where  the  wind  speed  and  source  level  change  drastically  with  posi¬ 
tion. 

Using  an  adiabatic  normal  mode  formulation  Perkins,  Kuperman,  Ingenito,  Fi- 
alkowski,  and  Glattetre[95]  extended  the  model  of  Kuperman  and  Ingenito  to  range- 
dependent  source  fields  and  mildly  range-dependent  waveguides.  They  did  this  by 
dividing  the  surface  area  into  smaller  sub-areas  over  which  the  source  field  could  be 
considered  constant.  They  used  far-held  approximations  for  each  sub-area.  These 
were  coupled  with  the  further  approximation  that  the  cross-spectral  density  for  each 
sub-area  could  be  expressed  as  a  single  sum  over  modes.  This  approximation  is  only 
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Region  of 
Hurricane  Winds 


Region  of 
Ambient  Winds 


Figure  2-2;  Cross  section  of  the  stratified  ocean  waveguide  showing  the  geometry 
of  the  surface  noise  problem  (Not  to  scale).  On  the  surface  is  the  area  covered  by 
the  hurricane  and  surrounding  area  covered  5  m/s  ambient  winds.  The  surface  noise 
sources  are  modeled  as  a  plane  of  monopoles  a  small  depth  ^o  below  the  surface  and 
the  sound  field  is  measured  by  receivers  and  receiver  arrays. 
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valid  when  the  inverse  of  the  difference  between  the  horizontal  wavenumber  of  the 
modes  is  much  less  than  the  dimension  of  the  sub-area.[108,  109]  For  the  highly  range- 
dependent  winds  of  a  hurricane  in  an  otherwise  range-independent  waveguide,  this 
approach  proves  to  be  less  accurate,  more  cumbersome  and  less  efficient  to  implement 
than  direct  integration. [140]  Carey,  Evans,  Davis  and  Botseas[19]  have  developed  a 
computational  approach  based  on  the  parabolic  equation  approximation  for  calculat¬ 
ing  range-dependent  surface  noise.  We  find  that  steep  angle  contributions  dominate 
the  intensity  measured  by  a  single  sensor  and  so  require  direct  integration  of  local 
noise  sources  with  a  full-field  model  for  the  Green  function  rather  than  an  elevation- 
ajigle-restricted  parabolic  approximation. 

It  is  useful  to  briefly  derive  the  direct  integration  approach  used  here  since  it  has 
not  explicitly  appeared  in  the  previous  literature  even  though  many  essential  elements 
are  implicit  in  the  work  of  Perkins,  Kuperman,  Ingenito,  Fialkowski  and  Glattetre.[95] 
For  uncorrelated  sources  the  cross-spectral  density  of  the  noise  field  can  be  written 
as 


C'(ri,r2,/)  =  ^o!Vo^^^^^2^^p(ri|ro,/)(/*(r2]ro,/)  (2.2) 

as  shown  in  Appendix  A  where  Sqq{V {po),  f)  is  the  source  power-spectral  density, 
which  is  a  function  of  wind  speed  V  and  frequency  /,  AA  is  a  smedl  area  increment  of 
integration  at  least  the  size  of  the  horizontal  coherence  area  of  the  source  distribution, 
and  ^(rj|ro,/)  is  the  waveguide  Green  function.  Throughout  this  paper  a  cylindrical 
coordinate  system  in  used  where  r  =  {p,z)  =  {p,0,z),  p  is  the  horizontal  location 
vector,  p  is  distance  from  the  origin,  9  is  azimuth  angle  and  z  is  depth  measured 
with  positive  downward  from  the  surface.  The  locations  ri  and  r2  are  receivers  and 
To  is  the  source.  Green  functions  are  calculated  by  a  combination  of  wavenumber 
integration  at  short  ranges  and  the  normal  mode  approximation  at  long  ranges.  The 
integration  over  surface  source  area  is  computed  numerically.  This  expression  is  valid 
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for  range-dependent  source  fields  and  environments. 

The  source  depth  Zq  is  taJcen  to  be  a  quarter  wavelength  for  all  simulations  in  the 
present  paper.  This  follows  noise  modeling  convention [70,  19,  95]  since  source  depths 
of  a  quarter  wavelength  or  less  lead  to  a  downward-directed  dipole  source  radiation 
pattern.  Hamson  has  shown  that  on  average  wind  generated  noise  in  the  ocean 
radiates  with  a  downward  directed  pattern  that  closely  fits  a  dipole  for  wind  speeds 
between  5  and  20  m/s  and  frequencies  from  400  Hz  to  3.2  kHz.  [50]  This  is  true 
even  for  average  source  depths  greater  than  a  quarter  wavelength  and  sea-surface 
roughness  much  larger  than  the  wavelength [50,  101]  as  in  a  hurricane  where  wave 
heights  may  exceed  10  m.  This  is  understandable  since  surface  noise  is  believed  to 
arise  from  many  monopole  sources,  in  particular  bubbles,  randomly  distributed  near 
the  sea  surface.  All  of  these,  by  the  method  of  images,  have  main  downward  directed 
lobes  which  reinforce  and  varying  horizontally  directed  side-lobes  which  cancel. 

As  discussed  in  the  introduction,  the  source  power-spectral  density  has  been  shown 
to  follow 


S„(V'./)  =  So(/)V'”'^>  (2.3) 

for  certain  frequency  and  wind  speed  ranges.  While  experiments[102]  at  wind  speeds 
below  20  m/s  give  n  =  3. 1^0. 3,  values  in  the  broader  n  =  1  to  n  =  4  range  will  be  used 
here  for  illustrative  purposes.  If  it  is  later  found  that  wind  speed  and  noise  intensity 
are  related  by  some  other  function,  the  power-law  relationships  considered  here  will 
provide  a  basis  for  piecewise  construction  of  this  more  complicated  dependence. 

Farmer  hcis  shown  experimentally  that  clouds  of  bubbles  near  the  ocean  surface 
may,  through  scattering  and  absorption,  lower  ambient  noise  levels  at  frequencies 
above  8  kHz  and  wind  speeds  above  15  m/s.[40]  While  such  attenuation  has  never 
been  observed  at  lower  frequencies,  we  will  consider  its  possibility  in  the  high  winds 
of  a  hurricane. 
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Attenuation,  in  dB/m,  can  be  written  as  a  =  101og(e)(T7a„,  where  o  is  the  extinc¬ 
tion  cross  section  of  an  individual  bubble  and  n„  is  the  number  of  bubbles  per  unit 
volume. [139]  Using  this  expression,  Weston[139]  provides  a  model  for  attenuation  by 
sea  surface  bubble  clouds,  based  on  the  extinction  cross  section  and  spatial  distri¬ 
bution  of  wind-generated  bubbles  as  a  function  of  wind  speed  and  frequency.  This 
attenuation  can  then  be  included  in  the  Green  function  in  Eq.  2.2  to  determine  its 
effect  on  the  underwater  noise  field.  This  is  done  by  calculating  the  Green  function 
for  a  waveguide  with  an  effective  attenuation  in  dB/m  of 


'  9.35  X 

2.44  X  10^^ fV^ 


f  <  1.5kHz 
f  >  1.5kHz 


in  a  layer  at  the  sea  surface  as  given  by  Weston. [139] 


(2.4) 


2.4  Single  Hydrophone  Analysis 

Here  it  is  shown  that  the  noise  intensity  measured  by  a  single  sensor  in  a  hurricane 
is  dominated  by  local  sea-surface  sources  rather  than  sound  propagating  from  longer 
ranges.  Underwater  acoustic  intensity  can  then  be  used  to  estimate  the  wind  speed 
within  a  local  resolution  area  since  wind  speed  in  a  hurricane  is  also  found  to  be 
effectively  constant  over  this  scale. 

Beginning  with  the  cross-spectral  density  of  the  noise  field  in  a  hurricane,  Eq. 
(2.2),  the  spectral  intensity  of  the  sound  field  received  at  r  can  written  as 


/(r./)  = 


C(r,r,/)  _  r  S„{V{po),f) 


Pw^ 


/oo 

-OO 


PyjCAA 


|p(r|ro,/,U(po))P  (2.5) 


where  the  total  instantaneous  intensity  is  given  by 
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(2.6) 


/(r)=  ri(T,f)df 

JQ 

The  Green  function  t/(r|ro, /,  V(po))  depends  on  local  wind  speed  V'(po)  because 
it  includes  attenuation  due  to  wind-generated  sea-surface  bubbles.  We  show  that 
this  wind  speed  dependence  is  negligible  at  frequencies  less  than  100  Hz  for  typical 
hurricane  wind  speeds,  but  needs  to  be  accounted  for  at  higher  frequencies.  Surface 
wind  speed  V  is  given  by  the  Holland  model  of  Fig.  2-1  for  a  hurricane,  while  the 
surrounding  ambient  wind  speed  is  taken  to  be  5  m/s. 

Two  hurricane-prone  ocean  environments  surrounded  by  densely  populated  coastal 
communities,  the  North  Atlantic  and  the  Bay  of  Bengal,  are  considered.  Their  sound 
speed  profiles  are  shown  in  Fig.  2-3.  The  difference  in  water  depth  between  these  two 
environments  leads  to  fundamental  differences  in  propagation.  Typical  near-surface 
sound  sources  will  lead  to  refractive  propagation  with  excess  depth  in  the  North 
Atlantic  but  not  in  the  Bay  of  Bengal.  In  the  former,  sound  may  propagate  efficiently 
to  long  ranges  via  the  deep-sound  channel,  while  in  the  latter,  it  will  multiply  reflect 
off  the  lossy  bottom  leading  to  far  greater  transmission  loss.  Although  hurricanes 
decrease  the  temperature  of  the  local  sea  surface  by  roughly  1°C  near  the  eye  wall  to 
roughly  35  m  depth,  the  corresponding  small  change  in  sound  speed[55]  of  roughly 
4  m/s  is  also  local  and  so  has  a  negligible  effect  on  the  curvature  of  both  local  and 
long-range  sound  paths. 

The  spectral  intensity  level,  given  by 


in  dB  re  Irefif),  of  hurricane  generated  noise  is  computed  by  the  direct  integration  of 
Eq.  2.5  as  a  function  of  receiver  range  p  and  depth  z  from  an  origin  at  the  center  of 
the  hurricane  on  the  sea-surface.  For  convenience  in  the  present  paper  the  reference 
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Figure  2-3:  Sound  speed  profiles  c{z)  for  the  North  Atlantic[82]  and  the  Bay  of 
Bengal[69,  114].  The  bottom  has  a  density  of  1.38  g/cm  and  an  attenuation  of 
0.3  dB/A  corresponding  to  the  deep  silty  sediment  layers  of  the  Bay  of  Bengal[123,  48] 
and  the  North  Atlantic  Abyssal  plain[49,  125].  The  water  has  a  density  of  1  g/cm 
and  an  attenuation  of  6  x  10“^  dB/A. 
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level  Ire f if)  is  taken  to  be  the  spectral  intensity  at  a  reference  depth  Zref  =  200  m 
for  a  reference  10-m  altitude  wind  speed  of  Vref  =  5  m/s  over  the  entire  ocean 


/„/(/)  =  (2.8) 

where  T^ej  =  (p,  2re/)-  Noise  intensity  has  been  measured  for  5  m/s  wind  speed  in 
many  ocean  environments  and  at  similar  depths.[102,  24,  97,  20]  In  an  experimental 
scenario  other  reference  values  could  be  chosen. 

Spectral  intensity  level  is  shown  in  Fig.  2-4  for  frequencies  of  50,  400  and  3200  Hz, 
spaced  three  octaves  apart,  using  Eqs.  (2.3)-(2.8)  and  assuming  n  —  Z.  The  choice 
of  n  =  3  is  within  measured  power-laws[102]  and  has  been  chosen  out  of  convenience 
since  it  is  linearly  related  to  the  power  of  the  wind.  [34]  The  wind  speed  profile  of  the 
hurricane  and  surroundings  based  on  the  Holland  model  at  an  altitude  of  10  m  from 
the  sea  surface  is  also  plotted  with  the  spectral  intensity  level  at  a  depth  of  200  m. 
The  most  apparent  feature  in  Fig.  2-4(a)  and  (c)  is  the  effectively  linear  relationship 
at  low  frequency,  50  Hz,  between  spectral  intensity  level  Lj  and  the  log  of  the  wind 
speed.  This  is  roughly  independent  of  depth  as  can  be  seen  in  Figs.  2-4(b)  and  (d). 
At  the  higher  frequencies  shown  sea-surface  bubbles  significantly  attenuate  sound  in 
the  high  wind  speed,  eye- wall  region  of  the  hurricane  but  the  noise  still  follows  local 
wind  speed  with  a  more  complicated  nonlinear  dependence  as  will  be  shown  in  the 
next  section.  The  small  increase  in  level  in  the  North  Atlantic  outside  the  hurricane 
at  ranges  of  193  and  257  km  and  at  a  depth  of  4.7  km  is  caused  by  convergence 
zone  propagation  from  the  powerful  sources  in  the  eye  wall.  This  convergence  zone 
structure  indicates  an  efficient  mechanism  exists  for  the  long  range  propagation  of 
hurricane  noise  in  this  environment  that  will  be  considered  in  Sec.  2.5. 
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Figure  2-4:  Noise  spectral  level  (dB  re  Irej)  in  the  North  Atlantic  (a)  and  (b)  and  the 
Bay  of  Bengal  (c)  and  (d)  for  n  —  3.  (a)  and  (c)  show  the  level  cis  a  function  of  range 
at  a  depth  of  200  m  for  50,  400  and  3200  Hz  frequencies.  Lv{p)  =  10  log  {V {p)jVrej) 
is  plotted  for  comparison  where  Vref  =  5  m/s.  Ly  =  0  is  equivalent  to  V  =  5  m/s  and 
Lv'  =  10  is  equivalent  to  V  =  50  m/s.  (b)  and  (d)  show  the  level  as  a  function  of  range 
and  depth  at  50  Hz.  In  both  waveguide  environments  the  noise  level  closely  follows 
the  local  wind  speed.  In  the  North  Atlantic  there  is  a  convergence  zone  structure 
corresponding  to  sound  that  propagates  from  the  hurricane’s  eye  wall  via  refraction. 
Note  the  convergence  zone  near  the  surface  at  a  range  of  257  km  and  the  ray  vertex 
depth  of  4.7  km. 
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2.4.1  Local  Noise  Dominates 


The  effectively  linear  relationship  between  the  log  of  the  local  wind  speed  and  un¬ 
derwater  acoustic  spectral  intensity  shown  in  Fig.  2-4  suggests  a  possible  simplifying 
approximation  to  our  formulation.  In  particular  the  areal  integral  of  Eq.  (2.5)  can 
be  approximated  by  integrating  only  over  the  local  sources  in  the  hurricane.  These 
fall  within  a  disc  of  area  A  —  ttR^  centered  at  the  horizontal  location  of  the  receiver 
p  which  provides  the  dominant  contribution  in  the  exact  integral  (Eq.  (2.5)).  The 
spectral  intensity  can  then  be  approximated  as 


I(Tj)  «  (2.9) 

where  Pq  —  Po  ~  p.  Such  a  simplification  can  potentially  lead  to  errors  if  R  is  too 
small. 

To  quantify  the  potential  error  of  this  local  approximation,  the  approximate  equa¬ 
tion  (2.9)  is  evaluated  for  a  receiver  under  the  eye  wall  of  the  hurricane  where  wind 
speed  varies  most  drastically.  When  compared  to  the  exact  result  of  Eq.  (2.5),  we 
take  the  error  induced  by  the  local  approximation  to  be  negligible,  less  than  or  equal 
to  1  dB,  for  R  greater  than  a  minimum  length  Riocai-  The  error  as  a  function  of 
R  is  given  in  Fig.  2-5  where,  for  deep-water  environments,  Riocai  =  300  fo  2000  m 
depending  on  sensor  depth. 

It  is  noteworthy  that  the  deep-ocean  North  Atlantic  and  Bay  of  Bengal  error  curves 
closely  match  those  of  the  infinite  half-space.  This  shows  that  bottom  reflections  and 
variations  in  sound  speed  profile  do  not  have  a  significant  effect  on  Riocai  in  deep  water. 
For  a  bottom  mounted  sensor  in  a  typical  shallow  water  environment  Riocai  =  2  to 
3  km  in  the  50  to  400  Hz  range.  Our  computations  also  show  that  Riocai  does  not 
change  significantly  for  the  expected  source  power-spectral  densities  and  attenuations 
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Figure  2-5:  Error  induced  by  the  local  area  approximation  (Eq.  (2.9))  as  a  function 
of  loczil  source  area  radius  Riocai  for  a  single  sensor  under  the  maximum  winds  in  the 
eye  wall  of  a  hurricane.  Curves  are  shown  for  the  North  Atlantic  and  the  Bay  of 
Bengal  environments  used  in  this  paper  as  well  as  for  infinite  half-space  and  shallow 
water  continental  shelf  environments.  Plots  are  given  for  sensor  depths  of  100  m  ((a) 
and  (c))  and  800  m  ((b)  and  (d))  and  for  frequencies  of  50  Hz  ((a)  and  (b))  and  400 
Hz  ((c)  and  (d)).  While  these  plots  are  given  for  n  =  3  the  difference  for  values  n—1 
to  4  is  less  than  0.1  dB.  The  North  Atlantic,  Bay  of  Bengal  and  infinite  half-space 
environments  are  very  similar.  In  these  deep-water  environments,  for  the  shsillow 
100  m  sensor  depth,  we  see  that  for  Rioccu  over  roughly  300  m  the  error  from  this 
approximation  is  negligible.  For  the  deeper  800  m  sensor  depth,  the  Riocai  for  which 
the  error  is  negligible  is  roughly  2  km.  In  shallow  water  the  error  in  the  local  area 
approximation  is  higher  leading  to  a  larger  Riocai-  This  is  likely  due  to  the  strong 
reflection  of  sound  off  bottom.  In  deep  water  environments  bottom  reflections  have 
little  effect  and  most  of  the  sound  measured  by  a  receiver  propagates  via  direct  path 
from  the  surface  source. 
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considered  in  this  paper. 

The  wind  speeds  in  a  hurricane  do  not  change  significantly  over  Riocai  and  can 
be  approximated  as  constant  in  Eq.  (2.9).  This  leads  to  less  than  0.2  dB  additionail 
error  in  the  spectral  intensity  level,  which  can  then  be  approximated  as 


H^.f) 


S^{V{p),  /) 

Pw 


/  /  Podp'odO'o\9{T\To,f,V{p))\‘^ 

Jo  Jo 


S„(Vip),f)W{Tj,V{p)) 


(2.10) 


where  only  the  local  wind  speed  V{p)  directly  above  the  receiver  has  a  significant 
effect  on  both  the  source  factor  5„(V(p),/)  and  the  waveguide  propagation  factor 
W{r,  f,V{p)).  The  source  factor  is  universal  in  that  it  does  not  depend  on  propa¬ 
gation  parameters  and  should  be  the  same  for  any  waveguide  environment  so  long 
as  the  ocean  depth  greatly  exceeds  the  ocean-atmosphere  boundary  layer.  While 
the  propagation  factor  does  depend  on  the  environment,  ocean  waveguides  typically 
change  gradually  with  horizontal  position.  The  wind-speed-independent  functionality 
of  W(r,  /,  V(p))  should  then  be  effectively  constant  over  Riocai  and  over  the  horizon¬ 
tal  extent  of  a  hurricane,  on  the  order  of  100  km.  Both  factors  may  be  characterized 
numerically  or  empirically  to  develop  a  set  of  curves  to  estimate  wind  speed  from 
acoustic  intensity.  In  the  next  section  we  find  that  it  is  possible  to  simplify  these 
factors  and  develop  an  approximate  analytic  equation  for  wind  speed  estimation. 

The  approximate  Eq.  (2.10)  for  range-dependent  sources  and  potentially  range- 
dependent  waveguides  is  simileir  to  Kuperman  and  Ingenito’s[70]  exact  Eq.  (30)  for 
range-independent  sources  and  waveguides  in  that  spectral  intensity  is  the  product 
of  a  “universal  ambient  noise”  source  factor,  following  Ingenito  and  Wolf[57]  and 
here  defined  as  5„(V (p),  /),  and  a  waveguide  propagation  factor  W (r,  /,  V (p)).  The 
implicit  cissumption  of  formulations  of  this  kind  is  that  variations  in  source  depth  can 
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be  accounted  for  as  equivalent  variations  in  Sqg{V{p),  f).  This  is  consistent  with  the 
measured  dipole  behavior  of  ambient  noise  in  the  ocean.  [50] 

Taking  the  log  of  Eq.  (2.10)  leads  to  a  useful  approximate  equation  for  spectrzil 
intensity  level 


L,(rJ)  «  Ls(V(p),S)  +  Lw(tJ,V(p))  (2.11) 

in  dB  re  Ire/ if)  where 

MV(p),/)  =  .0log(^)^).  (2.12) 

L„(T,f,V(p))  =  (2.13) 


and 


Jrefif)  =  S,,{VrejJ)W{rref,f,Vref)  =  S,,{Vref ,  f)Wref{f)  (2.14) 

Here  Ls{V{p),  f)  is  a  universal  ambient  noise  source  term  that  is  independent  of 
waveguide  propagation  parameters,  while  f,  y{p))  is  a  waveguide  propagation 

term.  The  functional  dependencies  of  the  first  term  cam  be  determined  empirically  in 
any  waveguide  where  the  ocean  depth  greatly  exceeds  the  ocean-atmosphere  bound¬ 
ary  layer,  while  the  functional  dependencies  of  the  second  term  should  be  locally 
determined. 

If  Sqq{V{p),  f)  follows  a  powcr-law,  such  as  Eq.  (2.3),  then  universal  ambient 
noise  source  level  is  linearly  related  by 
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is(V'(p),/)  =  10n(/)log(^) 


(2.15) 


to  the  log  of  wind  speed.  The  slope  of  this  linear  relationship  lOn(f)  has  been 
previously  measured  in  the  13  Hz  to  14.5  kHz  frequency  range  and  1  to  20  m/s  wind 
speed  range. [118,  38,  102,  24,  20] 

To  estimate  wind  speed  from  ambient  noise  measurements  using  Eq.  (2.11),  the 
dependence  of  Z/w(r,  /,  V{p))  on  wind-dependent  attenuation  by  sea-surface  bubbles 
needs  to  be  established.  This  may  be  done  empirically,  numerically  or  analytically  as 
in  the  next  section. 

2.4.2  Separating  the  Effect  of  Attenuation  by  Bubbles  from 
Local  Waveguide  Propagation 

Analytic  expressions  are  derived  to  show  how  attenuation  can  be  separated  from  other 
waveguide  propagation  effects  so  that  Lw(r, /,  V(p))  can  be  split  into  a  universal 
wind-speed-dependent  attenuation  term  and  a  local  waveguide  calibration  term  that  is 
wind-speed  independent.  These  analytic  expressions  also  demonstrate  the  uniqueness 
of  a  wind  speed  estimate  based  on  broad-band  underwater  noise  measurements.  They 
also  enable  analytic  expressions  for  estimation  error  to  be  obtained  in  some  important 
cases. 

Underwater  spectral  intensity  level  is  calculated  over  a  range  of  wind  speeds  and 
frequencies  relevant  to  hurricane  quantification  as  illustrated  in  Fig.  2-6  using  the 
full  areal  integration  of  Eq.  (2.5).  The  spectral  intensity  level  exhibits  a  maxima  that 
depends  on  wind  speed  and  frequency.  For  wind  speeds  and  frequencies  below  this 
maxima,  attenuation  by  bubbles  is  negligible  so  that  Lw(r, /,  F(p))  is  only  a  func¬ 
tion  of  the  local  waveguide  environment  and  spectral  intensity  level  L/(r,/)  should 
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depend  on  the  log  of  wind  speed  only  through  Eq.  (2.15)  given  the  power-law  n  =  3 
assumption  of  the  simulation.  For  higher  wind  speeds  and  frequencies,  attenuation 
by  bubbles  is  significant  and  eventually  leads  to  a  roll-off  in  the  spectral  intensity 
so  that  Lw(r, /,  y(p))  is  a  separable  function  of  both  wind-speed-dependent  and 
wind-speed-independent  terms. 

While  the  dependence  of  spectral  intensity  on  wind  speed  and  frequency  including 
attenuation  by  bubbles  can  be  calculated  exactly  using  the  full  areal  integration  of  Eq. 
(2.5)  or  the  local  integral  approximations  of  Eqs.  (2.9)  or  (2.10),  a  useful  first-order 
approximation  leads  to  the  analytic  result 


where 


W(tJ,V{p))  =  Wc(t, 


'  -y? 


(2.16) 


Wo{tJ)  =  W{T,/,0) 


(2.17) 


and 


A{V{p),f,kr  =  0)  = 


sm{kzo) 

2*’r2oi^cos(fcL)  + 


27rui  „-ikL 


(2.18) 


is  the  downward  plane-wave  amplitude  for  a  source  in  an  attenuating  sea-surface 
bubble  layer  following  the  Pekeris  solution. [58]  The  complex  wavenumber  k  =  -V- 

i °2^iog(e) ^  (2-18)  where  a(y{p),f)  is  given  in  Eq.  (2.4). 

The  spectral  intensity  level  of  Eq.  (2.11)  can  then  be  approximated  as 


L;(r,  /)  «  Ls{V{p),  f)  +  lAiVip),  f)  -h  Lwo(r,  /)  (2.19) 


where 
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Figure  2-6;  Simulated  noise  spectral  level  (dB  re  Iref)  in  the  North  Atlantic  for 
range-independent  winds  as  a  function  of  wind  speed  and  frequency  including  at¬ 
tenuation  by  sea-surface  bubbles  assuming  n  =  3  from  Eq.  (2.5).  Below  100  Hz 
the  power-law  relationship  between  noise  intensity  and  wind  speed  is  unaffected  by 
bubble  attenuation  even  up  to  the  80  m/s  wind  speeds  of  a  hurricane.  As  frequency 
increases,  attenuation  affects  the  noise  level  at  progressively  lower  wind  speeds.  For 
a  given  frequency  the  noise  level  increases  linearly  with  wind  speed,  peaks,  and  then 
decays  exponentially. 
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i^(V'(p),/)  =  201og( 


2^\A(V(p),fX  =  0)\ 


Zo 


(2.20) 


and 


i..(r,/)  =  101og(^)  (2.21) 

The  approximation  of  Eq.  (2.19)  is  in  agreement  with  the  full  areal  integration  of  Eq. 
(2.5)  to  within  1  dB  for  frequencies  below  500  Hz  even  at  hurricane  wind  speeds  as 
shown  in  Fig.  2-7. 

By  splitting  the  local  waveguide  and  bubble  attenuation  effects  of  Lw(r, /,  V’(p)) 
into  two  terms,  La{V (p),  /)  and  Two(r,  /),  wind  speed  can  now  be  estimated  from  am¬ 
bient  noise  using  Eq.  (2.19),  where  L^(V(p),/)  is  a  universal  attenuation  term  that 
depends  on  local  wind  speed  but  like  Ls{V{p),  f)  is  also  independent  of  waveguide 
parameters.  The  last  term  of  Eq.  (2.19),  Lwo(r, /),  is  a  local  waveguide  calibration 
that  is  independent  of  wind  speed. 

At  frequencies  below  100  Hz  where  attenuation  a  due  to  bubbles  is  negligible  at 
hurricane  wind  speeds  LA{V{p),f)  goes  to  zero,  as  expected  from  Fig  2-6.  In  this 
important  case,  if  Sqq{V{p),  f)  follows  a  power  law,  Eq.  (2.19)  reduces  to  a  linear 
equation  in  the  log  of  wind  speed 


L,(r,  /)  «  10n{/)  >og(^)  +  /)  (2.22) 

where  10n(/)  is  a  universal  empirically  determined  slope  and  Lwo(r,/)  is  a  local 
calibration  intercept.  The  log  of  wind  speed  can  be  then  found  from  measurements 
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Figure  2-7:  (a)  Noise  spectral  level  (dB  re  Irej)  as  a  function  of  wind  speed  at  several 
frequencies,  assuming  n  —  3.  The  black  curves  show  the  attenuation,  caused  by 
bubbles,  at  50  Hz,  400  Hz  and  4  kHz.  The  range  of  wind  speeds  typical  of  a  hurricane 
is  also  shown,  (b)  Noise  spectral  level  curves  as  a  function  of  frequency  for  typical 
hurricane  wind  speeds  of  30,  50  and  80  m/s.  The  black  curves  show  the  full  areal 
integration  from  Eq.  (2.5)  and  the  grey  curves  show  the  first-order  approximation  of 
the  field  given  by  Eq.  (2.19)  with  Eqs.  (2.15),  (2.20)  and  (2.21). 
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of  ambient  noise  level  by  standard  linear  least  squares  estimation,  as  has  been  done 
in  Refs.  [118]  and  [38]  at  low  wind  speed. 

As  frequency  increases,  bubble-layer  thickness  exceeds  a  quarter  wavelength  and 
the  La{V{p),  f)  term  can  be  approximated  as 

La{V{p),  /)  «  -a(V(p).  f)L  (2.23) 

If  we  use  for  illustrative  purposes  the  L  =  1.2  m  layer  thickness  given  by  Weston, [139] 
then  Eqs.  (2.20)  and  (2.23)  agree  to  within  1  dB  above  300  Hz  and  to  within  2  dB 
between  100  and  300  Hz.  While  Weston  notes  that  the  assumption  of  a  bubble  layer 
of  constant  thickness  may  be  poor  at  high  wind  speeds,  any  future  improvements  in 
our  knowledge  of  the  parameter  L  can  be  incorporated  in  Eqs.  (2.18)  and  (2.23). 

The  locations  of  maxima  in  noise  spectral  level  correspond  to  the  ridge  in  Fig. 
2-6.  These  can  now  be  approximated  analytically  by  substituting  Eqs.  (2.15),  (2.21) 
and  (2.23)  into  Eq.  (2.19)  and  taking  the  derivative  with  respect  to  wind  speed  to 
obtain 


’  (1/(2.15  X  10-^Lv7))i/3 
(l/(5.63  X  10-^L/))i/3 


300  </  <  1.5kHz 
f  >  1.5kHz 


here  assuming  n  —  3  and  a{V{p),j)  from  Eq.  (2.4). 


(2.24) 


2.4.3  Accuracy  of  Underwater  Acoustic  Wind  Speed  Esti¬ 
mates 

By  standard  stationary  averaging,  it  should  be  possible  to  reduce  the  variance  of 
an  underwater  acoustic  wind  speed  estimate  enough  to  be  useful  for  meteorological 
purposes.  Given  the  relationship  V  —  H{I)  between  the  true  wind  speed  V  and  true 
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ambient  noise  intensity  /,  the  maximum  likelihood  estimate  (MLE)  of  the  wind  speed 
V  given  a  measurement  of  ambient  noise  intensity  /  is  V  =  H{f)  by  the  invariance 
of  the  MLE.  [60]  The  function  H  can  be  found  either  numerically  from  the  excict 
integration,  Eq.  (2.5),  or  analytically  from  one  of  the  approximations,  Eqs.  (2.11), 
(2.19)  or  (2.22).  We  define  the  percent  root-mean-square  error  (RMSE)  of  the  wind 
speed  estimate  V  as 


^rmae 


100 


v/<  \v  -  vp  > 


<  V  > 


and  the  percent  bias  as 


(2.25) 


t^bias  =  100 


I  <  V>  >  -T| 
V 


(2.26) 


given 


<  >=  r  H^ihplhdi  (2.27) 

where  p(/)  is  the  probability  density  function  of  the  measured  intensity  I.  For  the 
hurricane  noise  measurements  considered  here,  where  the  contributions  from  a  large 
number  of  independent  sources  are  received  simultaneously,  the  acoustic  field  is  ex¬ 
pected  to  be  a  circular  complex  Gaussian  random  variable.  The  time  averaged  mea¬ 
sured  intensity  I  is  then  expected  to  follow  a  Gamma  distribution[80,  81] 


=  - fM - 


(2.28) 
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where  fj.  is  the  time-bandwidth  product  and  /  is  the  mean  of  the  noise  measurement. 

Prom  the  full  areal  integration  of  Eqs.  (2.5)  we  can  numerically  find  the  percent 
RMSE  and  percent  bias  of  the  wind  speed  estimate  V.  For  frequencies  below  100 
Hz,  where  attenuation  o;  is  insignificant,  we  find  that  the  percent  RMSE  and  percent 
bias  are  functions  of  n  and  fj,  as  shown  in  Fig.  2-8.  At  higher  frequencies,  where 
attenuation  is  significant,  the  percent  RMSE  and  percent  bias  are  also  functions  of 
frequency  and  wind  speed.  This  is  illustrated  in  Fig.  2-9  at  a  frequency  of  400  Hz 
assuming  n  =  3. 

Following  the  standard  practice  of  stationary  averaging,  the  variance  of  noise 
measurements  is  reduced  by  inverse  the  number  of  stationary  samples,  1/fi.  In  typ¬ 
ical  ocean  acoustic  applications,  such  as  matched  filtering,  /i’s  in  excess  of  100  are 
common.  [128,  25,  17]  For  example  Piggott[102]  and  Perrone[96]  have  obtained  mea¬ 
surements  of  wind  noise  level  with  standard  deviations  less  than  1  dB  corresponding[99, 
80,  81]  to  /i  >  19. 

Given  a  spectral  intensity  measurement  with  n  >  19,  underwater  acoustic  wind 
speed  estimates  with  errors  similar  to  the  6%  to  15%  errors  of  hurricane- hunting 
aircraft[41]  are  possible.  For  example,  at  low  frequencies  where  attenuation  is  insignif¬ 
icant,  a  measurement  of  noise  spectral  level  with  ^  =  19  would  yield  a  corresponding 
percent  RMSE  in  estimated  wind  speed  of  6  to  25%  for  the  range  of  published  values 
for  n  as  shown  in  Fig.  2-8.  For  the  higher  frequency  400  Hz  example  in  Fig.  2-9, 
where  attenuation  is  significant,  a  spectral  intensity  measurement  with  fi  =  19  will 
yield  percent  RMSEs  from  9  to  20%.  Even  larger  errors  are  common  for  remote  satel¬ 
lite  techniques,  as  high  as  40%  as  noted  in  the  introduction.  From  this  error  analysis 
we  find  that  underwater  acoustic  measurements  may  be  worthwhile  for  estimating 
hurricane  wind  speed.  Additional  errors  related  to  the  practical  application  of  the 
underwater  acoustic  technique  will  be  discussed  in  Sec.  2.4.4. 

At  low  frequencies,  less  than  100  Hz,  where  attenuation  a  from  bubbles  becomes 
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Figure  2-8:  The  percent  RMSE  t'rmse  (a.)  and  percent  bias  i/b^as  (b)  of  the  wind 
speed  estimate  V  where  attenuation  by  sea-surface  bubbles  is  insignificant,  evaluated 
numerically  from  Eqs.  (2.5)  and  (2.27).  For  time-bandwidth  products  n  >  5  the 
estimate  becomes  unbiased  and  the  RMSE  attains  the  Cramer-Rao  Lower  Bound. 
Piggott[102]  and  Perrone[96]  have  measured  wind  noise  level  with  standard  deviations 
less  than  1  dB  which  corresponds  to  //  >  19.  For  /i  =  19  the  percent  RMSE  in  the 
wind  speed  estimate  ranges  from  6  to  25%  depending  on  n  which  is  a  significant 
improvement  over  the  primary  satellite  classification  method. 
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Figure  2-9:  The  percent  RMSE  Vrmse  (a)  and  percent  bias  Ubias  (b)  of  the  wind  speed 
estimate  V  including  the  effect  of  attenuation  calculated  numerically  from  Eqs.  (2.5) 
and  (2.27),  assuming  n  =  3,  at  /  =  400  Hz  where  Vmax  =  58  m/s.  The  error  and  bias 
increase  for  V  ^  Vmax  but  for  /u  >  5  and  for  values  of  V  where  Ubias  <1%  the  percent 
RMSE  decreases  and  attains  the  Cramer-Rao  Lower  Bound.  For  spectral  intensity 
measurements  with  //  =  19  the  percent  RMSE  in  this  example  is  between  9  and  20% 
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insignificant,  the  moments  of  V  can  also  be  evaluated  einalytically  from  the  first-order 
approximation  of  Eq.  (2.22)  to  illustrate  the  fundamental  parameters  affecting  a  wind 
speed  estimate.  The  mean  of  the  wind  speed  estimate  can  then  be  written  as 


Tip)  \SQWofi} 
and  the  standard  deviation  as 
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At  these  low  frequencies  the  percent  bias  can  then  be  approximated  as 


(2.30) 


^bias 


100 


Tjp+l/n) 

TiiAn'i" 


(2.31) 


and  the  percent  RMSE  as 
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These  analytic  expressions  for  the  percent  RMSE  and  percent  bias  match  those  cal¬ 
culated  numerically  from  Eqs.  (2.5)  and  (2.27)  and  shown  in  Fig.  2-8  to  within 
1%. 

At  low  frequencies,  where  attenuation  is  insignificant,  the  Cramer-Rao  Lower 
Bound  can  be  derived  from  the  first-order  approximation,  Eq.  (2.22),  as  shown  in 
App.  B.  This  provides  a  straightforward  analytic  method  for  calculating  the  percent 
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RMSE  as 


-irifsy^^^O’Oymptoticiy)  mn^  /oooN 

t'rmae  ~  ^ -  =  100 - —  (2.33) 

<V  >  riy/Jl 

which  matches  the  numerically  computed  value  in  Fig.  2-8  for  /x  >  5.  This  is 
expected  since  the  Cramer-Rao  Lower  Bound  is  the  asymptotic  variance  for  large  /x. 
The  Cramer-Rao  Lower  Bound  can  also  be  used  to  calculate  the  percent  RMSE  at 
frequencies  above  300  Hz  from  the  first-order  approximation  in  Eqs.  (2.19)  with  Eqs. 
(2.15),  (2.21)  and  (2.23)  yielding 
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which  matches  the  numerical  results  in  Fig.  2-9  when  /x  >  5  and  <  !%• 


2.4.4  Practical  Issues 

We  have  shown  that  a  single  underwater  acoustic  sensor  provides  significant  potential 
as  a  measurement  tool  to  accurately  estimate  local  wind  speed  in  a  hurricane.  There 
are  practical  issues,  however,  to  consider  when  deploying  such  sensors  to  monitor 
a  hurricane.  While  this  is  not  a  definitive  discussion  of  all  the  issues  that  might 
be  involved,  we  will  attempt  to  illustrate  how  an  underwater  acoustic  measurement 
system  might  be  implemented.  For  example,  how  would  one  deploy  these  sensors, 
how  many  sensors  would  be  needed  to  fully  characterize  a  hurricane,  and  how  much 
would  it  cost. 
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One  possible  scenario  would  be  to  deploy  multiple  sonobouys,  similar  to  those 
used  in  weather  quantification  experiments  by  Nystuen  and  Selsor,[90]  from  aircraft 
or  ships  in  the  path  of  an  oncoming  hurricane.  As  the  hurricane  passes  over  each 
sonobouy  the  sensor  would  cut  a  swath  through  the  storm  recording  the  wind  speeds 
overhead.  The  swaths  from  multiple  sonobouys  could  give  a  fairly  complete  measure¬ 
ment  of  the  wind  speeds  in  the  hurricane.  This  is  similar  to  the  current  measurements 
made  by  hurricane-hunting  aircraft  which  fly  through  the  storm  cutting  a  swath  and 
measuring  wind  speed.  For  both  methods,  sonobouys  or  hurricane-hunting  aircraft, 
the  sensors  must  pass  through  the  eye  wall  of  the  hurricane  where  the  winds  are 
strongest.  For  aircraft  this  means  actively  piloting  the  plane  through  the  storm, 
whereas  with  stationary  sonobouys,  one  would  deploy  many  sensors  along  a  line  that 
crosses  the  expected  path  of  the  hurricane  to  insure  that  at  least  one  sonobouy  cuts 
through  the  eye  wall.  For  example,  a  line  of  20  sonobouys  spaced  5  km  apart  aeross 
the  hurricane’s  path  would  span  almost  100  km  assuring  several  measurements  of  the 
wind  speed  in  the  eye  wall.  Sonobouys  may  be  ideal  for  surviving  in  a  hurricane 
since  most  of  the  sensor  floats  at  a  depth  of  several  hundred  meters  and  only  a  small 
surface  expression  would  be  exposed  to  hurricane  winds. 

The  advantage  of  deploying  sonobouys  in  advance  of  a  hurricane  is  that  the  ship 
or  aircraft  never  has  to  enter  the  storm  and  would  not  need  to  be  as  expensive  as 
the  specialized  hurricane-hunting  aircraft  used  today.  The  cost  of  a  typical  hurricane¬ 
hunting  aircraft  such  as  the  WC-130  is  $78  million  (inflation  adjusted  to  year  2003 
dollars)  [2]  and  the  cost  of  a  single  flight[9]  is  roughly  $155,000.  Between  two  and  eight 
aircraft  flights  are  made  per  day  [41]  for  potentially  landfalling  hurricanes  in  the  North 
Atlantic  where  the  lifespan  of  a  hurricane  can  be  several  weeks.  Twenty  sonobouys, 
at  $500  each, [85]  could  be  deployed  from  inexpensive  non-specialized  ships  or  aircraft 
in  the  path  of  8in  oncoming  hurricane  well  before  conditions  are  dangerous  for  roughly 
$10,000. 
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An  alternative  scenario  would  be  to  deploy  hundreds  of  permanent  shore-cabled 
hydrophone  systems,  at  $10,000  to  $20,000  each  depending  on  cable  length,  in  strate¬ 
gic  hurricane  prone  areas  for  a  few  million  dollars.  As  noted  before,  this  is  much  less 
than  the  purchase  price  of  a  WC-130  hurricane-hunting  aircraft. 

Such  underwater  acoustic  systems  would  likely  be  used  in  conjunction  with  apriori 
location  estimates  from  satellites.  Satellites  would  determine  the  path  of  the  hurricane 
relative  to  the  hydrophone  and  show  whether  the  sensor  passed  through  the  high 
winds  of  the  eye  wall.  The  underwater  acoustic  measurement  would  then  provide  an 
estimate  of  the  wind  speeds  for  the  portions  of  the  hurricane  that  passed  overhead.  If  a 
hydrophone  does  not  pass  through  the  powerful  eye  wall  but  rather  through  the  weaker 
surrounding  winds  it  would  still  provide  a  lower  bound  or  threshold  measurement  of 
wind  speed  and  it  may  be  possible  to  extrapolate  these  lower  wind  speeds  to  determine 
the  higher  wind  speeds  of  the  eye  wall. 


2.5  Hydrophone  Array  Analysis 

The  analysis  in  the  previous  sections  demonstrates  how  omnidirectional  sensors  may 
be  used  to  accurately  measure  the  local  winds  and  quantify  the  destructive  power  of 
a  hurricane  as  it  passes  overhead.  It  may  be  possible  to  use  arrays  of  hydrophones 
to  beamform  on  the  acoustic  field  from  a  hurricane  at  long  range.  For  illustrative 
purposes  we  will  consider  horizontal  linear  arrays  of  the  type  that  might  be  towed 
from  an  oceanographic  or  naval  vessel,  however,  other  array  configurations,  such  as 
moored  arrays,  might  also  be  useful.  Arrays  might  also  be  useful  for  directionally 
filtering  out  other  noise  sources,  such  as  ships  and  surf,  in  local  measurements. 

Using  the  expression  for  the  cross-spectral  density  of  the  noise  field  of  Eq.  (2.2) 
we  find  the  angular  spectral  density  of  the  noise  received  by  an  TV  element  array,  or 
beamformed  output,  to  be 


62 


(2.35) 


■'''  m=l  n=l 

with  units  of  /iPa^/sr^Hz,  where  k  =  ^in,  in  is  a  unit  vector  in  the  steering  direction 
f2,  and  Tm  is  the  position  of  the  hydrophone  on  the  array. 

We  define  the  hurricane  wind-generated  noise  source  area  to  include  sources  within 
200  km  of  the  hurricane’s  center  as  shown  in  Fig.  2-1  and  the  ambient  noise  source 
area  to  include  sources  generated  by  the  5  m/s  winds  surrounding  the  hurricane.  To 
shov/  how  an  array  might  be  able  to  measure  the  destructive  power  of  a  hurricane,  the 
angular  spectral  density  of  the  noise  will  be  calculated  for  a  hurricane  as  a  function 
of  majcimum  wind  speed. 

The  angular  spectral  density  of  Eq.  (2.35)  in  the  direction  of  the  hurricane  in¬ 
creases  with  maximum  wind  speed,  as  shown  in  Fig.  2-10  for  an  array  at  200-m  depth 
at  a  range  far  from  the  hurricane  eye.  The  difference  in  spectral  density  between  the 
strong  72-m/s-wind-speed  and  weak  33-m/s- wind-speed  hurricanes  of  Fig.  2-10  is 
roughly  10  dB  given  the  assumption  here  that  n  =  3.  The  difference  in  spectral 
density  would  be  greater  for  larger  n. 

A  practical  horizontal  array  can  resolve  the  important  features  of  the  hurricane 
such  as  the  eye  wall,  which  has  dimensions  of  tens  of  kilometers,  when  placed  in  a 
convergence  zone  as  in  Figs.  2-10(a)  and  (c).  This  is  not  possible  for  an  array  just 
outside  the  convergence  zone  as  shown  in  Fig.  2-10(b).  In  the  former  case,  the  length 
L  of  an  array  oriented  at  broadside  to  the  hurricane,  would  have  to  be 

L  >  RX/l  (2.36) 

where  R  is  the  range  from  the  array  to  the  hurricane  and  I  is  the  size  of  the  eye- 
wall.  Typical  linear  arrays[128]  have  lengths  L  on  the  order  of  100^.  In  the  example 
of  Fig.  2-10,  a  broadside  array  with  L  =  32A,  similar  to  the  ONR  FORA  array, [1] 
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Angle  from  Broadside  (deg) 


Figure  2-10:  Angular  spectral  density  level  10log{B{^,  f))  (dB  re  at  100  Hz 

for  a  64-element  A/2-spaced  horizontal  broadside  array  as  a  function  of  steering  angle 
for  hurricane  generated  noise  in  the  North  Atlantic  at  ranges  of  257  km  (a),  289  km 
(b),  and  385  km  (c)  from  the  eye  of  the  hurricane,  assuming  n  —  3.  Ranges  of  257  and 
385  km  correspond  to  the  fourth  and  sixth  convergence  zones  from  the  center  of  the 
hurricane.  289  km  is  exactly  between  the  fourth  and  fifth  convergence  zones.  Curves 
are  shown  for  a  powerful  72  m/s  hurricane,  a  medium  50  m/s  hurricane  and  a  weak 
33  m/s  hurricane.  The  angular  spectral  density  level  from  ambient  noise  is  plotted 
for  comparison.  A  steering  angle  of  0°  corresponds  to  the  array  steered  toward  the 
calm  eye  of  the  hurricane  and  the  powerful  eye  wall  is  located  at  ^3°.  This  array 
has  an  angular  resolution  of  1.8°  which  at  a  range  of  257  km  corresponds  to  a  8  km 
spatial  resolution.  64 


images  the  hurricane  with  10  km  resolution  at  a  range  of  320  km.  The  width  of  the 
convergence  zone  must  also  be  sufficiently  small  in  rcinge  to  resolve  the  eye  wall.  For 
the  given  environment  and  ranges  considered,  this  condition  is  satisfied  because  the 
convergence  zone  width  in  range  is  roughly  5  km,  which  is  less  than  the  width  of  the 
eye  wall. 

A  horizontal  array  oriented  at  end-fire  to  the  hurricane  has  the  advantage  that 
it  discriminates  against  local  surface  noise  coming  from  near  broadside  in  favor  of 
sound  that  travels  from  long  distances  at  shallow  angles  in  the  waveguide.  This  could 
potentially  lead  to  longer  hurricane  detection  ranges.  Unfortunately,  at  end-fire,  the 
length  of  the  array  must  satisfy 


L  >  2X{R/lf  (2.37) 

to  resolve  the  eye-wall.  For  example  an  impractically  long  L  —  2000A  end-fire  array 
would  be  needed  to  achieve  10  km  resolution  at  a  range  of  320  km. 

The  analysis  presented  here  for  the  North  Atlantic  shows  that  it  may  be  possible 
to  image  the  features  of  a  hurricane  using  linear  broadside  arrays  of  sufficient  length. 
Waveguides  that  lack  excess  depth,  such  as  the  Bay  of  Bengal,  do  not  exhibit  the 
convergence  zone  structure  seen  in  the  North  Atlantic.  This  probably  makes  it  ex¬ 
tremely  difficult  to  even  detect  hurricanes  using  practical  linear  arrays  at  long  ranges 
in  these  environments. 


2.6  Conclusions 

We  have  shown  that  the  wind-generated  noise  received  by  a  single  underwater  acous¬ 
tic  sensor  in  a  hurricane  can  be  well  approximated  by  sea-surface  contributions  so 
local  that  wind  speed  and  surface  source  intensity  can  be  taken  as  nearly  constant. 
Two  terms  with  empirically  and  analytically  determined  dependencies  may  be  used 
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to  estimate  wind  speed  from  measured  ambient  noise  spectral  level  (1)  a  universal 
ambient  noise  source  term  and  (2)  a  local  waveguide  calibration  term.  At  low  fre¬ 
quencies,  current  evidence  suggests  a  simple  power-law  relationship  exists  between 
noise  intensity  and  wind  speed  so  that  the  log  of  wind  speed  may  be  estimated  accu¬ 
rately  from  spectral  ambient  noise  level  by  linear  least  square  estimation.  At  higher 
frequencies,  a  non-linear  relationship  is  expected  but  we  show  that  it  should  be  pos¬ 
sible  with  broadband  measurements  to  make  unambiguous  low-variance  wind  speed 
estimates  from  broad-band  noise  measurements. 
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Chapter  3 


Quantifying  Wind  Speed  in 
Hurricane  Gert  with  Undersea 
Sound 

3.1  Introduction 

In  this  chapter  we  demonstrate,  based  on  experimental  measurements,  that  inexpen¬ 
sive  ocean  acoustic  sensors  may  be  used  to  accurately  quantify  the  destructive  power 
of  a  hurricane  by  measuring  the  wind  speeds  to  within  5%  accuracy.  This  is  similar  to 
the  accuracy  achieved  by  hurricane-hunting  aircraft.  [41]  To  do  this  we  correlate  the 
undersea  noise  generated  by  hurricane  Gert, [71]  recorded  by  an  autonomous  acoustic 
sensor, [37]  with  meteorological  data  from  reconnaissance  aircraft  and  satellites.  From 
this  correlation  we  find  that  the  intensity  of  low  frequency  sound  generated  by  the 
hurricane  is  approximately  proportional  to  the  cube  of  the  local  wind  speed. 

Accurate  hurricane  quantification  is  critical  for  disaster  planning  in  order  to  miti¬ 
gate  the  destructive  effect  of  these  storms.  The  destructive  power  of  a  tropical  cyclone 
was  recently  demonstrated  by  hurricane  Katrina  which  caused  over  1000  fatalities[29] 
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and  an  estimated  economic  loss  of  100  billion  dollars[10].  Prior  to  Katrina,  the  United 
States  Commission  on  Ocean  Policy  emphasized  the  need  for  more  ciccurate  quantifi¬ 
cations  of  hurricane  destructive  power  to  improve  disaster  planning.  [133] 

As  discussed  in  Chapter  1,  the  standard  method  for  satellite  hurricane  classifica¬ 
tion  is  the  Dvorak  method  [31]  where  hurricane  cloud  features  observed  in  satellite 
images  are  interpreted  to  estimate  wind  speed  and  classify  destructive  power.  This 
method  can  yield  errors  in  wind  speed  estimates  exceeding  40%  [93,  43,  12,  4,  121] 
(when  compared  to  the  best-estimate  wind  speed  from  aircraft  measurements).  For 
example,  of  the  eight  North  Atlantic  hurricanes  of  2000,  three  of  them  [93,  43,  12] 
experienced  Dvorak  errors  over  40%  and  three  more  [94,  72,  120]  experienced  Dvorak 
errors  over  20%.  Several  satellite  microwave  techniques  show  some  promise  for  mea¬ 
suring  hurricane  wind  speed  [59]  but  because  of  resolution  and  accuracy  issues  the 
Dvorak  method  is  still  the  standard  for  satellite  hurricane  classification  [44] .  More  ac¬ 
curate  classification  can  be  achieved  by  flying  specialized  ‘hurricane-hunting’  aircraft 
through  the  high  winds  of  a  hurricane,  but  their  expense  prohibits  their  routine  use 
outside  of  the  United  States  [33].  For  example,  the  cost  to  purchase  a  WC-130  aircraft 
is  roughly  $82  million  [2]  (adjusted  for  inflation  to  year  2005  dollars)  and  the  deploy¬ 
ment  cost  is  $155,000  per  flight  [9].  Given  these  limitations  in  current  capability, 
the  United  States  Commission  on  Ocean  Policy  has  recommended  that  future  ocean 
observing  systems  be  used  to  improve  weather  related  warnings  [133].  We  experi¬ 
mentally  demonstrate  that  underwater  acoustic  measurements  of  noise  intensity  may 
provide  a  useful  technique  for  hurricane  wind  speed  estimation  and  quantification. 

The  physical  mechanisms  for  ocean  noise  generation  by  wind  are  not  entirely  un¬ 
derstood.  The  entrainment  of  oscillating  bubbles  by  wind-generated  waves  is  believed 
to  be  the  dominant  natural  source  mechanism  at  frequencies  between  10  Hz  and  10 
kHz  with  a  broad  peak  between  300  and  600  Hz. [67,  136,  102,  23]  The  non-linear 
interax:tions  of  ocean  waves  has  been  proposed  as  the  dominant  mechanism  between 
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0.1  and  10  Hz  with  a  peak  at  roughly  0.3  Hz. [134,  22]  In  the  10  to  50  Hz  frequency 
range  considered  here,  it  is  likely  that  the  oscillation  of  wave-entrained  bubbles  is 
the  dominant  effect.  Regardless  of  the  specific  physical  mechanism,  the  relationship 
between  noise  intensity  and  wind  speed  has  been  extensively  studied  experimentally 
at  low  non-hurricane  wind  speeds, [118,  38]  but  never  previously  at  hurricane  wind 
speed  to  our  knowledge  due  to  the  extraordinary  difficulty  of  conducting  experiments 
at  sea  in  hurricane  conditions. 

The  entrainment  of  bubbles  may  play  a  role  in  the  attenuation  of  sound  in  the 
ocean,  particularly  at  high  frequencies  and  wind  speeds. [40]  At  the  low  frequen¬ 
cies  and  hurricane  wind  speeds  considered  here,  attenuation  by  entrained  bubbles 
is  expected  [141]  and  found  to  be  negligible. 


3.2  Methods 

In  1999  the  United  States  National  Oceanographic  and  Atmospheric  Administra¬ 
tion’s  (NOAA)  Pacific  Marine  Environmental  Laboratory  (PMEL)  moored  an  under¬ 
sea  acoustic  sensor  at  a  depth  of  roughly  800  m  in  the  4.7  km  deep  North  Atlantic  at 
17.7°N  49.5°W  with  the  purpose  of  recording  seismic  events  along  the  mid-Atlantic 
ridge. [37]  On  September  15th,  1999  this  hydrophone  recorded  the  underwater  sound 
generated  by  hurricane  Gert[71]  as  it  passed  overhead  (Fig.  3-1). 

The  measured  acoustic  intensity  in  the  10  to  50  Hz  band,  shown  in  Fig.  3-2, 
changed  with  time  first  reaching  a  maximum,  then  a  minimum,  and  then  another 
maximum.  This  variation  is  consistent  with  the  advection  of  the  characteristic  ge¬ 
ometry  of  the  hurricane  over  the  hydrophone  with  the  first  maximum  in  the  sound 
intensity  corresponding  to  the  leading-edge  eye  wall  winds,  the  minimum  by  the  low 
wind  speeds  of  the  eye,  and  the  later  maximum  by  the  trailing  edge  of  the  eye  wall. 

In  order  to  reduce  the  uncertainty  of  these  noise  measurements,  we  average  the 
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Spectragram  at  NOAA  hydrophone  during  hurricane  Gert 


08:10  14:45  21:19 

Time(HH;MM) 


Figure  3-1:  Spectrogram  of  the  hurricane  noise  in  dB  re  //Pa^/Hz,  from  1  to  50  Hz, 
received  by  the  NOAA  hydrophone  on  15  Sept.  Above  10  Hz  the  maximum  between 
13:00  and  16:00  correspond  to  the  powerful  winds  of  the  hurricane’s  eye-wall.  At  1  Hz 
we  see  noise  most  likely  due  to  cable  strum.  There  seems  to  be  some  hurricane  related 
noise  between  1  and  10  Hz,  however  it  does  not  correlate  well  with  wind  speed  and 
may  be  caused  by  nonlinear  surface  wave  interaction  as  discussed  in  the  next  chapter. 
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Figure  3-2:  Log  of  the  noise  spectral  intensity  I  (solid  blue)  in  Watts/m^Hz,  from  10 
to  50  Hz,  received  by  the  NOAA  hydrophone  on  15  Sept.  The  expected  intensity  (red 
dashed  line)  is  based  on  the  hindcast  wind  speed  V  in  m/s  above  the  hydrophone 
given  the  best  fit  relation  log[//(lWatt/m^Hz)]  =  nlog[K/(lm/s)]  -1-6  where  n  =  3.36 
and  6  =  —14.46  .  The  maxima  at  13:30  and  15:30  correspond  to  the  powerful  winds 
of  the  hurricane’s  eye-wall  and  the  minimum  at  14:30  corresponds  to  the  hurricane’s 
eye.  The  percent  root-mean-square  error  between  the  two  curves  is  5%. 
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received  intensity  over  the  10  to  50  Hz  band  and  over  10  minute  consecutive  time 
windows.  This  leads  to  a  time-bandwidth  product  /i  =  24000  which  can  be  thought 
of  as  the  number  of  independent  samples  in  the  average.  By  the  method  of  stationary 
averaging[80]  the  standard  deviation  at  of  the  log  L  =  log  (//Watt/m^Hz)  of  the 
averaged  intensity  /  can  be  given  as  at  —  log(e)Y/l/^  for  //  >>  1.  For  the  log  of  the 
averaged  intensities  shown  in  Fig.  1,  the  standard  deviation  is  then  gl,  =  0.003  dB 
re  Watt/m^Hz. 

No  direct  measurements  of  wind  speed  were  taken  at  the  sensor’s  location  so  we 
synthesize  the  local  wind  speed  record  based  on  available  data  from  United  States 
Air  Force  (USAF)  hurricane-hunting  aircraft  and  NOAA  hurricane  tracking  satel¬ 
lites.  We  begin  by  estimating  the  hurricane’s  track  to  determine  its  location  relative 
to  the  acoustic  sensor  over  time.  NOAA  satellites  gave  location  estimates  of  Gert 
every  six  hours  with  error  margins[115]  of  db  28  km  and  USAF  aircraft  gave  location 
estimates[30]  with  error  margins[41]  of  ±  11  km.  These  two  sets  of  location  esti¬ 
mates  yield  two  different  possible  tracks  hurricane  Gert  could  have  taken  past  the 
acoustic  sensor.  Based  on  satellite  data  hurricane  Gert  passed  11  km  to  the  South  of 
the  acoustic  sensor  at  14:15  on  15  Sept.,  however,  based  on  aircraft  data  the  storm 
passed  32  km  to  the  South  of  the  acoustic  sensor  at  13:49.  To  resolve  this  ambiguity 
we  consider  an  ensemble  of  possible  tracks  hurricane  Gert  could  have  taken  past  the 
sensor  spanning  the  error  margins  in  the  measurements,  with  distances  from  50  km 
to  the  South  to  20  km  to  the  North  of  the  sensor  and  times  from  12:00  to  16:00.  Also, 
we  consider  a  range  of  typical  values  for  hurricane  translation  speed,  from  5  to  15 
m/s. 

In  addition  to  providing  hurricane  location  estimates,  USAF  hurricane-hunting 
aircraft  were  used  to  measure  the  wind  speed  structure  in  the  hurricane, [30]  normalized[104] 
to  an  altitude  of  10  m.  The  measurements  taken  on  16  Sept.,  a  day  after  the  hurricane 
passed  the  acoustic  sensor  and  the  first  day  aircraft  measurements  were  made,  show 
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a  wind  speed  structure  typical  of  a  hurricane  with  a  low  wind  speed  eye  surrounded 
by  the  high  47  m/s  wind  speeds  of  the  eye  wall  (Fig.  3-3).  To  reconstruct  the  wind 
speed  record  at  the  undersea  sensor  location  as  Gert  passed  overhead,  we  advect  this 
wind  speed  structure  backward  in  time  in  the  directions  determined  by  the  ensemble 
of  possible  hurricane  tracks.  Even  though  the  aircraft  wind  speed  measurement  was 
taken  nearly  a  day  after  the  acoustic  measurement,  it  should  be  possible  to  hindcast 
the  wind  speed  record  since  evidence  indicates  that  the  hurricane  wind  speeds  re¬ 
mained  roughly  constant  during  that  time.  This  can  be  seen  from  satellite  estimates, 
which  indicate  a  less  than  10%  increase  in  maximum  hurricane  wind  speed[71],  and 
from  NO  A  A  meteorological  forecast  models,  which  predict  no  change  in  maximum 
wind  speed[115],  over  the  24-hour  period. 

For  each  possible  hurricane  track  the  power-law  hypothesis  is  tested  by  calculating 
the  linear  regression  between  the  log  of  the  intensity  I  and  the  log  of  the  local  wind 
speed  V  (Eq.  (3.2)).  The  best-fit  hurricane  track  is  taken  as  the  one  where  the  root- 
mean-square  error  (RMSE)  of  the  linear  regression  is  minimized.  For  this  best-fit 
track,  hurricane  Gert  passed  29  km  to  the  South  of  the  NOAA  undersea  acoustic 
sensor  at  14:15  moving  at  a  speed  of  about  12  m/s.  This  is  within  the  error  margins 
of  the  aircraft  and  satellite  location  estimates.  Given  this  best-fit  track  the  log  of  the 
local  wind  speed  above  the  sensor  may  be  plotted  over  time  (Fig.  1). 


3.3  Results 

From  this  regression  analysis,  wind  speed  V  and  undersea  noise  intensity  /,  which 
is  a  scalar  magnitude  proportional  to  the  magnitude  square  of  the  complex  acoustic 
pressure,  follow  a  simple  power  law  relationship 

1{V)  =  V’^W  (3.1) 


73 


Wind  speed  (m/s)  for  Hurricane  Gert,  Sept  16  1999 
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Figure  3-3:  Wind  speed  in  m/s  in  hurricane  Gert  as  a  function  of  position  relative  to 
the  center  or  ‘eye’.  As  hurricane  Gert  passed  over  the  acoustic  sensor  the  sensor  effec¬ 
tively  sliced  a  path  through  the  hurricane.  The  lines  show  the  paths  estimated  based 
on  Satellite  (dcish-dot),  Aircraft  (dashed)  and  Acoustic  best  fit  (solid)  respectively. 
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as  shown  in  Fig.  3-4,  where  is  a  waveguide  propagation  factor.  [141]  Taking  the 
logarithm  this  becomes  a  linear  relationship 


'°«(i  wlu/lm.)  =  "'“*(1^) 

where  n  =  3.36  and  b  —  —14.46  in  the  low  frequency  band  of  the  acoustic  mea¬ 
surements  (Fig.  2).  The  slope  n  is  universal  and  independent  of  measurement 
position,  while  the  intercept  6  is  a  calibration  that  depends  on  local  waveguide 
environment .[  1 4 1  ] 

This  correlation  between  undersea  noise  and  local  wind  speed  demonstrates  that 
acoustic  sensors  can  be  used  as  anemometers  for  estimating  hurricane  wind  speed. 
This  is  illustrated  in  Fig.  1  where  the  hurricane  wind  speed  closely  follows  the  log  of 
the  noise  intensity.  Additionally,  an  acoustic  sensor  deployed  to  pass  under  a  hurri¬ 
cane’s  eye  wall  could  estimate  the  storms  total  destructive  power  because  hurricane 
power  is  proportional  to  the  cube  of  the  maximum  wind  speed.  [53] 

The  log  of  wind  speed  can  be  then  found  from  measurements  of  ambient  noise 
level  by  standard  linear  least  squares  estimation,  as  has  been  done  by  Shaw,  Watts 
and  Rossby[118]  and  Evans,  Watts,  Halpern  and  Bourassa[38]  for  low  (<  10  m/s) 
spatially-uniform  wind  speeds.  The  empirical  fit  in  Fig.  2  shows  that  wind  speed  can 
be  determined  from  undersea  noise  intensity  to  within  a  fractional  error  of  0.05,  or 
5%,  where  fractional  error  is  defined  as 


^/<\v-v\-‘> 

u  =  ^ - - -  (3.3) 

and  where  V  is  the  estimated  wind  speed  from  undersea  acoustic  measurements  and 
V  is  the  actual  wind  speed.  This  is  similar  to  the  5  m/s  errors  typical  of  aircraft 
estimates. [41]  It  should  be  noted  that  satellite  estimates[116]  of  Gert’s  wind  speed 
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Log(V  / 1  m/s) 


Figure  3-4:  The  log  of  the  noise  spectral  intensity  I  in  Watts/m^Hz  versus  the  log 
of  the  surface  wind  speed  V  in  m/s  at  the  NOAA  hydrophone  (circles)  based  on  the 
best  fit  hurricane  track.  The  best  fit  linear  regression  shows  a  log[//(lWatt/m^Hz)]  = 
3.361og[V/(lm/s)]  —  14.46  relationship  between  intensity  and  wind  speed  with  a  per¬ 
cent  root-mean-square  error  of  only  5%. 
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were  much  less  accurate  and  overestimated  the  wind  speed  by  roughly  40%  when 
compared  with  aircraft  estimates[30]  made  near  the  time  of  the  acoustic  measurement. 

The  results  of  Figs.  3-2  and  3-4  are  based  on  the  ‘best  fit’  track  which  is  the 
hurricane  track  for  which  the  RMSE  of  the  linear  regression  is  smallest.  The  other 
possible  tracks  that  hurricane  Gert  could  have  taken  must  also  be  considered.  Figure 
3-5  shows  the  RMSE  and  slope  n  of  the  linear  regression  for  a  range  of  other  possible 
hurricane  tracks.  While  the  minimum  RMSE  of  0.76  dB  corresponds  to  a  slope 
n  =  3.36,  there  is  a  range  of  tracks  with  small  RMSEs,  less  than  0.8  dB,  which 
correspond  to  slopes  n  between  3  and  3.6.  We  see  then  that,  due  to  the  uncertainty 
in  our  knowledge  of  wind  speed  at  the  hydrophone  location,  there  is  some  uncertainty 
in  the  slope  of  the  relationship  between  log  of  wind  speed  and  log  of  acoustic  intensity. 

This  uncertainty  can  be  somewhat  alleviated  be  comparing  the  data  from  hurri¬ 
cane  Gert  with  data  gathered  by  other  researchers  at  low  non-hurricane  wind  speeds. 
Figure  3-6  shows  the  linear  regression  between  the  log  of  intensity  and  log  of  wind 
speed  including  low-wind-speed  data  gathered  by  Cato  and  Tavener  [23].  With  this 
additional  data  the  slope  of  the  regression  becomes  n  =  3.34  which  is  effectively  the 
same  as  that  given  in  Fig.  3-4. 

It  should  be  noted  that  in  Fig.  3-6  the  variance  or  spread  in  the  data  from 
Cato  and  Tavener  [23]  is  significantly  greater  than  in  the  data  from  hurricane  Gert. 
There  are  three  factors  that  may  contribute  to  this  increased  variance.  First,  the 
time-bandwidth  product  for  the  data  samples  from  Gert  is  significantly  greater,  by  a 
factor  of  more  than  100,  than  for  the  data  from  Cato  and  Tavener.  This  results  in 
greater  stationary  averaging  of  the  data  from  Gert.  Second,  the  acoustic  sensor  used 
for  Gert  was  at  a  depth  of  800  m  compared  to  40  m  for  the  measurements  of  Cato  and 
Tavener.  This  leads  to  increased  spatial  averaging  over  more  surface  noise  sources 
in  the  hurricane  data.  Last,  the  low-wind-speed  measurements  may  be  corrupted 
by  other  non-wind-related  noise  sources.  These  three  effects  may  contribute  to  the 
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Figure  3-5:  The  upper  plot  gives  the  root-mean-square  error  (RMSE)  of  the  linear 
regression  for  each  possible  hurricane  track.  The  y-axis  is  the  north-south  distance 
between  the  hydrophone  and  the  possible  hurricane  track  and  the  x-axis  is  the  speed 
at  which  the  hurricane  was  moving.  The  minimum  RMSE  of  0.76  dB,  indicated  y  the 
-f  symbol,  is  taken  to  represent  the  ‘best  fit’  hurricane  track.  Note  that  the  RMSE 
is  small,  less  than  0.8  dB,  over  a  range  of  possible  tracks.  The  lower  plot  gives  the 
slope  n  of  the  linear  regression  for  each  possible  hurricane  track  with  the  same  x  and 
y  axes  as  the  upper  plot.  For  the  best  track  the  slope  n  =  3.36.  For  the  range  of 
tracks  with  a  RMSE  less  than  0.8  dB,  the  slope  of  the  regression  could  vary  from  3 
to  3.6. 
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Figure  3-6:  Log  of  noise  intensity  I  versus  log  of  hindcast  local  wind  speed  V  based 
on  the  best-fit  track  (red  circles).  Also  plotted  is  the  log  of  noise  intensity  /  ver¬ 
sus  log  of  local  wind  speed  V  measured  by  Cato  and  Tavener  [23]  (blue  circles) 
at  lower  non-hurricane  wind  speeds  off  Australia.  The  linear  regression  shows  a 
log(//l  Watt/m^Hz)  =  3.341og(F/l  m/s)  —  14.5  relationship  between  noise  intensity 
/  and  wind  speed  V. 
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increase  accuracy  of  the  acoustic  measurements  in  hurricane  Gert. 


3.4  Conclusions 

Our  theoretical  calculations[141]  show  that  an  omnidirectional  sensor  at  the  depth 
of  the  NOAA  sensor  (800  m)  in  a  deep-water  environment  like  the  North  Atlantic 
effectively  measures  the  noise  generated  by  sea-surface  winds  within  a  horizontal 
radius  (?a  2  km)  that  is  much  smaller  than  the  typical  length  scales  of  a  hurricane 
making  it  ideal  for  hurricane  wind  speed  measurements.  A  number  of  distributed 
point  sensors,  similar  to  those  used  in  weather  quantification  experiments  by  Nystuen 
and  Selsor,[90]  could  be  scattered  from  aircraft  or  ships  in  the  path  of  an  oncoming 
hurricane.  Each  sensor  would  then  cut  a  swath  through  the  storm  recording  the 
wind  speeds  overhead.  This  is  analogous  to  the  procedure  followed  by  hurricane¬ 
hunting  aircraft  which  fly  an  ‘X-pattern’  through  a  hurricane  cutting  two  swaths 
while  measuring  wind  speed. [41] 

The  advantage  of  deploying  sonobouys  in  advance  of  a  hurricane  is  that  the  ship  or 
aircraft  never  has  to  enter  the  storm  and  would  not  need  to  be  as  expensive  as  the  spe¬ 
cialized  hurricane-hunting  aircraft  used  today.  Several  sonobouys,  at  $500  each  (M. 
May,  personal  communication),  could  be  deployed  from  inexpensive  non-specialized 
ships  or  aircraft  in  the  path  of  an  oncoming  hurricane  well  before  conditions  are  dan¬ 
gerous.  While  it  is  unlikely  that  undersea  acoustic  techniques  would  render  aircraft 
obsolete,  this  technique  would  provide  inexpensive  hurricane  monitoring  capability  for 
areas  of  the  world  not  currently  covered  by  specialized  hurricane-hunting  aircraft  such 
as  the  Paciflc  and  Indian  Oceans.  In  the  North  Atlantic  where  specialized  aircraft 
axe  already  in  use,  acoustic  techniques  may  make  it  possible  to  reduce  the  number  of 
aircraft  flights  necessary. 

Based  on  measurements  of  undersea  noise  from  hurricane  Gert,  we  demonstrate 
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a  new  method  for  estimating  hurricane  wind  speed  and  power.  The  power-law  rela¬ 
tionship  and  high  correlation  between  noise  intensity  and  local  hurricaine  wind  speed 
show  that  estimates  of  hurricane  wind  speed  based  on  undersea  noise  measurements 
could  be  very  accurate  with  an  error  margin  of  only  5%.  An  undersea  acoustic  sensor 
could  then  effectively  act  as  an  acoustic  anemometer,  providing  an  accurate  estimate 
of  a  hurricane’s  wind  speeds  and  destructive  power.  Currently  many  ocean  acoustic 
systems  exist  that  could  be  used  for  meteorological  measurements,  such  as  the  NOAA 
sensor  used  in  this  work,  and  additional  systems  could  be  deployed  from  ships,  air¬ 
craft,  or  near  shore  at  relatively  low  cost. 
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Chapter  4 


Modeling  Microseism  Generation 
by  Inhomogeneous  Ocean  Surface 
Waves  in  Hurricane  Bonnie  Using 
the  Non-linear  Wave  Equation 


4.1  Introduction 

Hurricanes  generate  seismic  noise,  commonly  referred  to  as  microseisms,  in  the  0.1  to 
0.6  Hz  frequency  range.  Here  we  describe  the  microseisms  generated  by  the  spatially 
inhomogeneous  waves  in  a  hurricane.  Using  the  ocean  surface  directional  wave  spec¬ 
trum  in  hurricane  Bonnie  [142,  88],  we  hindcast  the  microseismic  field  and  compare 
it  with  seismic  measurements  from  Florida.  Previously  hurricane  surface  directional 
wave  spectra  had  not  been  adequately  measured  or  modeled  so  that  researchers  had 
to  rely  on  assumed  spectra  in  order  to  model  hurricane  microseisms  [78] .  This  analysis 
is  useful  because  microseisms  are  a  primary  cause  of  noise  in  seismic  measurements 
[74,  107,  105]  that  raise  the  detection  threshold  for  monitoring  earthquakes  [135] 
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and  tsunamis.  Historically  microseismic  measurements  have  also  been  used  to  track 
hurricanes  [107,  45]  although  that  task  is  now  accomplished  using  satellites. 

We  present  an  analytic  expression  for  microseism  generation  by  spatially  inhomo¬ 
geneous  waves  typical  in  a  hurricane  based  on  the  non-linear  wave  equation  where 
a  second-order  seismo-acoustic  field  is  generated  by  a  source  distribution  which  de¬ 
pends  on  the  first-order  ocean  surface  wave  motion.  The  seismo-acoustic  field  at  a 
receiver  can  then  be  expressed  as  the  integral  over  the  source  distribution  multiplied 
by  the  waveguide  Green  function.  This  expression  is  ideal  for  hurricane  generated 
microseisms  since  it  can  be  used  to  calculate  the  acoustic  field  due  to  spatially  in¬ 
homogeneous  surface  waves.  Also,  this  expression  may  be  used  in  range-dependent 
waveguide  environments  as  is  the  case  when  a  hurricane  at  sea  generates  microseisms 
that  propagate  up  the  continental  margin  to  a  receiver  on  land. 

In  the  past,  microseism  generation  has  been  modeled  for  idealized  hurricanes  [78] 
where  the  ocean  surface  wave  field  within  the  hurricane  was  taken  to  be  spatially 
homogeneous  and  the  ocean  waveguide  to  be  range-independent.  Previous  microseis¬ 
mic  models,  however,  cannot  be  applied  in  typical  hurricane  scenarios  because  the 
surface  wave  fields  are  inhomogeneous  in  that  the  wave  height  spectra  in  different 
parts  of  the  hurricane  can  vary  both  in  magnitude  and  directionality.  In  some  cases 
this  is  because  spatial  homogeneity  is  assumed  over  infinite  [51,  16,  54,  77,  62]  or 
very  large  [78,  63,  64]  surface  wave  areas.  One  model  that  does  account  for  finite 
microseismic  source  generation  regions  is  limited  by  the  assumption  that  the  receiver 
is  at  the  center  of  the  source  area  [21,  22].  These  model  requirements  are  summarized 
in  Table  4.1.  Also,  these  previous  models  are  not  applicable  to  the  range-dependent 
environments  typical  in  hurricane  measurements  where  microseisms  generated  by  a 
hurricane  over  the  deep  ocean  are  measured  by  a  sensor  on  land.  Our  non-linear  wave 
equation  method  is  shown  to  agree  with  earlier  models  [54,  77]  if  we  make  the  same 
simplifying  assumptions  that  the  source  generation  region  is  spatially  homogeneous 
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Table  4.1:  Previous  microseism  models  and  their  assumptions 
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Longuet-Higgins  1950 
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Hasselman  1963 
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Hughes  1976 
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and  that  the  waveguide  can  be  modeled  as  an  infinite  halfspace. 

While  the  treatment  of  inhomogeneous  surface  wave  fields  and  range-dependent 
waveguides  presented  here  had  not  been  previously  considered,  the  medium  nonlinear 
wave  equation  used  here  has  a  parallel  derivation  to  the  wave  equations  or  hydrody¬ 
namic  equations  used  in  earlier  works.  Longuet-Higgins  [78],  Hasselmann  [51]  and 
Brekhovskikh  [16]  base  their  derivations  on  perturbing  Bernoulli’s  equation  while 
Hughes  [54]  and  Kibblewhite  [62]  begin  by  separately  perturbing  the  equations  of 
momentum,  continuity  and  state.  We  show  that  the  non-linear  wave  equation  used 
here  is  equivalent  to  perturbation  expansions  [54]  for  the  physical  parameters  typical 
in  hurricane  microseisms.  Lloyd  [77]  and  Cato  [21]  base  their  derivations  on  Lighthill’s 
equation,  however,  Lloyd  [77]  shows  that  both  perturbation  and  Lighthill  approaches 
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yield  the  same  end  result.  We  also  note  that  the  non-linear  wave  equation  can  be 
derived  as  a  second-order  approximation  of  Lighthill’s  equation  [138]. 

These  second  order  non-linear  or  perturbation-based  theories  should  not  be  con¬ 
fused  with  the  linear  theories  proposed  by  Banerji  [5]  and  Bowen  et  al  [14].  These 
linear  theories  claim  that  the  first-order  motion  of  a  surface  gravity  wave  creates  a 
first  order  pressure  fluctuation  on  the  sea  floor  regardless  of  how  deep  the  ocean  is. 
This  is  contrary  to  classic  surface  wave  theory  which  shows  that  first-order  wave  mo¬ 
tion  decays  exponentailly  with  depth  such  that,  in  deep  water,  the  first-order  pressure 
fluctuation  on  the  sea  floor  goes  to  zero  [6,  78,  27]  (also  see  App.  C). 

Based  on  the  wave-height  spectra  in  hurricane  Bonnie,  we  calculate  the  microseis- 
mic  source  levels  generated  by  the  nonlinear  interaction  of  the  ocean  surface  waves. 
Our  derivation  is  in  agreement  with  earlier  works  [86,  78],  which  show  that  micro¬ 
seisms  are  generated  by  the  non-linear  interaction  of  ocean  surface  waves  with  roughly 
the  same  wavelength  but  opposing  propagation  directions.  Recent  measurements 
[142,  132]  and  models  [88]  of  surface  directional  wave  height  spectra  in  hurricane 
Bonnie  show  the  opposing  surface  waves  necessary  to  generate  microseisms. 


4.2  Ocean  Surface  Gravity  Waves 

Hurricanes  are  characterized  by  high  winds  that  can  vary  quickly  with  position,  both 
in  direction  and  speed,  as  shown  in  Fig.  4.2(A).  In  addition  a  hurricane  may  move  at 
speeds  up  to  15  m/s  [33]  so  that  the  winds  at  any  location  can  also  change  with  time. 
These  spatially  and  temporally  varying  winds  generate  complex  ocean  surface  wave 
directional  spectra  (Fig.  4.2)  with  wave  heights  that  can  exceed  10  m  [142,  132,  88]. 

For  surface  gravity  waves  that  are  homogeneous  over  an  infinite  ocean  surface  area 
we  can  express  the  complex  surface  wave  height  as  the  linear  superposition  of  plane 
waves  [68]  where  k  is  the  sea  surface  wavenumber  and  a  is  the  corresponding  frequency 
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Figure  4-1:  (A)  Wind  speed  in  m/s  and  (B)  microseismic  source  level  {lo,  Zg  —  0) 
at  cj  =  1.26  rad/s  (/  =  0.2  Hz)  in  dB  re  Pa^s/m^  from  Eq.  (4.39)  at  1200  on  24 
August  1998  as  a  function  of  latitude  and  longitude.  The  arrow  indicates  the  direc¬ 
tion  hurricane  Bonnie  wais  moving.  The  letters  (a),  (b),  (c)  and  (d)  represent  features 
of  interest;  (a)  indicates  the  location  of  maximum  wind  speed,  (b)  indicates  the  eye 
where  wind  speed  is  zero,  and  (c)  and  (d)  indicate  peaks  in  the  microseismic  source 
level.  This  figure  shows  that,  while  a  hurricane  can  produce  significant  microseis¬ 
mic  source  levels  (B),  these  level  do  notSiirectly  follow  the  wind  speeds  (A)  in  the 
hurricane. 
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Figure  4-2:  The  wave-height  power  spectral  level  (lOlog Smm{K)/{rn^/{ra(f/m'^))) 
in  dB  re  m^/(rad^/m^)  at  the  locations  of  interest  (a),  (b),  (c),  and  (d)  given  in  Fig. 
4.2  at  1200  on  24  August  1998.  The  peak  in  the  spectra  is  at  roughly  k  =  0.04 
rad/m  or  a  =  0.63  rad/s  which  corresponds  to  an  acoustic  frequency  of  w  =  1.26 
rad/s  or  /  =  0.2  Hz.  At  some  locations  (b,  c,  and  d)  there  are  waves  propagating 
with  opposing  wavenumber  vectors  k,  while  at  other  locations  (a),  most  of  the  waves 
propagate  in  the  same  direction.  From  Eq.  (4.37)  we  expect  these  locations  with 
opposing  waves  to  produce  the  greatest  microseismic  source  levels  and  shown  in  Fig. 
4.2(B). 
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where  in  deep  water  cr^  —  qk.  Throughout  this  paper  the  surface  wavenumber  k  is 
expressed  either  in  cartesian  coordinates  as  or  in  polar  coordinates  as  (/c,  0). 

Also,  a  cartesian  coordinate  system  is  used  for  position  where  r  =  (p,  z)  =  (x,  y,  z) 
where  2  is  defined  downward  from  the  sea  surface.  The  complex  surface  wave  height 
can  then  be  written  as 


00 

Qp,  t)  =  II  A{K)e-^^+'*^  P(fK 

—  00 


(4.1) 


where  A(k)  is  the  surface  wave  height  spectrum. 

Since  the  complex  wave  height  C(P)^)  is  produced  by  the  contributions  of  many 
independent  random  physical  processes,  we  can  assume  by  Central  Limit  Theorem 
that  the  statistical  distribution  of  C(p,t)  is  Gaussian  [66].  In  addition  we  may  define 
the  reference  position  from  which  wave  height  is  measured  such  that  C(p,  t)  has  zero 
mean.  Given  this  assumption,  in  App.  D.l  we  show  that  A{k)  must  also  be  a  zero- 
mean  complex  Gaussian  random  variable  [66,  68]. 

Since  the  spectrum  A{k)  is  zero-mean  and  homogeneous  over  an  infinite  ocean 
surface  area,  the  second  moment  of  the  spectrum  can  be  written  as  (see  App.  D.l) [66, 
84] 


<  A(k)A*(k')  >=  S{k)6{k  -  k'),  (4.2) 

where  <>  denotes  the  expectation.  This  expression  shows  that,  for  stationary  pro¬ 
cesses,  the  wavenumber  components  k  and  k'  decorrelate.  Also,  since  the  spectrum 
A(/c)  is  Gaussian,  the  fourth  moment  can  then  be  written  as 
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5(/Cl  -  Kz)S{Ki)S{k2  -  K4)S{K2) 

+(5(k;i  -  k;4)5(/ci)^(k;2  -  K3)S{k2).  (4.3) 

In  a  hurricane  the  surface  waves  are  not  homogeneous  over  an  infinite  area,  but 
rather  the  surface  wave  spectrum  A{k)  changes  gradually  with  position  in  the  storm. 
To  characterize  the  spatially  varying  surface  wave  field  in  the  hurricane,  we  divide 
the  sea  surface  into  a  grid  made  up  of  finite  regions  m  with  surface  areas  Sm  over 
which  the  wave  spectrum  Am{K)  can  be  taken  to  be  homogeneous.  The  complex  wave 
height  in  a  particular  finite  region  m  then  can  be  expressed  as 


00 

Up,t)=  jj  (44) 

-oo 

In  App.  D.2  we  show  that  the  finite  size  of  the  regions  m  introduces  a  ‘windowing’ 
effect  when  calculating  the  moments  of  the  spectrum  However,  we  also  show 

that  this  windowing  effect  can  be  neglected  if  the  dimensions  and  Ly  of  the  region 
are  much  greater  than  the  wavelength  of  the  surface  wave.  For  the  hurricane  waves 
of  interest  in  this  paper,  the  surface  wavelengths  Xg  range  from  roughly  100  to  300  m 
requiring  that  the  dimensions  Lx  and  Ly  be  on  the  order  of  1  km  or  greater.  Provided 
this  condition  is  met  we  can  make  the  approximation. 


<  Am{K)A*^{K')  Smn{K)S{K  -  k') 


(4.5) 


where  Smn{i^)  is  the  cross-spectral  density  of  the  surface  wave  fields  in  regions  m  and 
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n. 

As  before  we  can  assume  by  Central  Limit  Theorem  that  Cm(p,  0  is  a  zero-mean 
complex  Gaussian  random  variable  so  that  Am{K)  is  also  a  zero-mean  complex  Gaus¬ 
sian  random  variable.  Because  of  this  we  can  write  the  fourth  moment  as 


<  Am{Kx)Am{K2)A*^{Ki)A*^{Ki)  >« 

(5(Ki  ^3)‘S'mn(^l)<^(^2  ^4)'5'mn(^2) 

-(-<5(Ki  -  Ki)Smn{K-\)K>^2  “  >^i)S,nn{l<-2) ■  (4.6) 

This  expression  for  the  fourth  moment  of  the  wave  height  spectra  will  be  used  later 
in  Sec.  4.4. 

From  Eq.  (4.4)  we  can  write  the  real  part  of  the  surface  wave  height  as 


1  °° 

U(p,  t)  =  ^jj  (4.7) 

—  oo 

We  can  also  express  the  real  surface  wave  particle  velocity  v  as  a  linear  superpo¬ 
sition  of  plane  waves  such  that 


v™(r,i)=  jj  +  + 

—  OO 

-  Yi.(A„»(«)e--‘+‘«-^  -  A;(K)e-‘-’«-^)e-''^d2^  (4  8) 

where  Zj,  iy  and  represent  unit  vectors  in  the  x,  y  and  ^  directions  respectively. 

Wright  et  al  [142]  and  Walsh  et  al  [132]  measured  the  spatial  variation  of  the  direc¬ 
tional  wave  spectra  Am{K.)  in  hurricane  Bonnie  using  an  aircraft-mounted  Scanning 
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Radar  Altimeter  (SRA).  These  measurements,  however,  are  limited  to  the  locations 
and  times  of  the  aircraft  flights  and  do  not  give  a  complete  picture  of  the  wave  field 
in  the  hurricane. 

To  quantify  the  surface  wave  field  at  any  location  and  time  in  hurricane  Bonnie, 
Moon  et  al  [88]  used  the  Wave  Watch  III  (WW3)  ocean  surface  wave  model  [126]. 
These  model  results,  shown  in  Fig.  4.2,  are  in  close  agreement  with  measurements 
from  Wright  et  al  [142]  and  Walsh  et  al  [132]  and  with  data  from  bouys  and  oceano¬ 
graphic  stations.  In  this  paper  we  use  the  surface  directional  wave  spectra  A„,(k;)  in 
hurricane  Bonnie  calculated  using  WW3  [88]. 

Experiments  by  Forristall  [42]  show  that  the  linear  wave  theory  of  Eqs.  (4.1, 
4.4,  4.7  and  4.8)  is  adequate  to  describe  surface  waves  even  in  high  sea  states  with 
significant  wave  breaking  as  in  a  hurricane.  The  use  of  linear  surface  wave  theory  can 
also  be  justified  since  the  typical  surface  wave  heights  in  a  hurricane  are  an  order  of 
magnitude  less  than  the  wavelength  [42,  142,  132,  88]. 


4.3  Nonlinear  Wave  Equation 

The  non-linear  wave  equation  [137,  28]  describes  the  second-order  acoustic  wave  field 
generated  by  first-order  fluid  motion.  For  the  long  surface  waves  in  a  hurricane 
viscosity  can  be  neglected  so  that  the  non-linear  wave  equation  can  be  written  as  [89] 


V^P2(r,t) 


1  d‘^P2{T,t)  ^  F 

pc^  dt'^ 

(vi(r,t)  •  vi(r,t)) 


dt^ 
p  d"^ 


-V^(vi(r,t)  •  vi(r,t)). 


(4.9) 


where  the  first  order  velocity  vi  and  pressure  pi  terms  on  the  right  hand  side  generate 
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a  second  order  acoustic  pressure  field  p2  on  the  left  hand  side,  c  is  the  sound  speed, 
p  is  density  and  F  is  the  coefficient  of  nonlinearity  [28,  89]. 

In  this  work,  we  consider  the  case  where  the  first  order  velocity  vj,  from  Eq. 
(4.8),  is  due  to  surfax^e  wave  motion.  The  first  order  pressure  p\  can  be  found  using 
conservation  of  momentum  where  —  Vpi  =  pd\\/dt.  Using  the  relationship  between 
first-order  velocity  and  pressure  we  can  compare  the  relative  magnitudes  of  the  source 
terms  in  Eq.  (4.9).  Given  typical  parameters  k  —  0.02  to  0.7  rad/m,  a  =  0.4  to 
0.8  rad/s,  c  ss  1500  m/s,  p  ~  lOOO/c^/m^,  and  F  ~  —2.6  [13]  we  find  that  the 
^V^(vi(r,  <)  ■  Vi(r,  t))  term  of  Eq.  (4.9)  exceeds  the  other  source  terms  by  several 
orders  of  magnitude.  Dropping  the  lesser  terms  yields 


V^P2(r,  t)  -  V^(vi(r,  t)  •  vi (r,  t))  =  ^(r,  t).  (4.10) 

This  equation  is  equivalent  to  the  perturbation  expression  used  in  the  previous  mi¬ 
croseism  derivation  of  Hughes  [54].  This  can  be  seen  by  taking  Eq.  (10)  of  Hughes 
[54]  and  assuming  irrotational  flow  (V(vi(r,<)  •  Vi(r,t))  =  2(vi(r,t)  •  V)vi(r,t))  and 
continuity  (V  •  (/9oVi(r,  t))  =  0).  The  assumption  of  irrotational  flow  is  justifiable  for 
linear  ocean  surface  wave  fields. 

Taking  the  Fourier  Transform  of  Eq.  (4.10)  we  find  the  frequency  domain  Helmholtz 
equation  for  the  second  order  field 


V^P2(r,w)  -f  A:^P2(r,u;)  =  q{v,u).  (4.11) 

We  will  begin  by  solving  for  the  pressure  field  generated  by  a  finite  volume  Vm 
with  a  source  distribution  which  depends  on  the  fluid  velocity  from  Eq.  (4.8). 
The  volume  Vm  extends  vertically  from  just  below  the  ocean  surface  at  =  0  to  a 
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depth  below  the  surface  wave  region  (z  >  Xg  =  21: /k).  The  volume  extends  over  the 
horizontal  area  with  dimensions  Lx  and  Ly  centered  at  =  {xm,  Vm)  as  defined 
in  Sec.  4.2.  We  can  solve  the  Helmholtz  Equation  (Eq.  (4.11)),  using  the  frequency 
domain  Green  function  g{Tr,Ta,u)  such  that, 


^’2,m(l'r,w)  =  JJJ  g{Tr,T„u)qm(T„U>)<i‘ 


(4.12) 


where  and  =  {xr,yr,Zr)  are  the  source  and  receiver  locations 

respectively. 


Xg  OO 

A,m(rr,w)  j  dZs  jj  W{pa  -  p^)g{Tr,Ta,  w)qm{Ts,  w)cF p^,  (4.13) 

0  -OO 

where  the  window  function  w(p,  —  p^)  is  unity  for  —Lx/2<Xa—Xm<  Lx  12  and 
—Ly/2  <  Va  -  Vm  <  Ly/2  and  zero  elsewhere.  The  double  integral  over  pa  represents 
the  field  at  a  receiver  due  to  a  distribution  of  sources  q  over  a  finite  area  defined  by 
the  window  function  w.  In  effect  this  finite  area  acts  as  a  horizontal  planar  array.  If 
the  receiver  Fr  is  in  the  far  field  of  the  finite  area  (jr^  —  r„i|  >  L\/Xa  and  |rr  —  r,„l  > 
L'^lXa  where  Xa  is  the  acoustic  wavelength  and  where  =  {p^,Za)  =  (a^m)  J/m, -Zi)), 
we  can  make  the  plane  wave  approximation. 


g(Tr,Ta,w)7^  g{Vr,Tm,w)e  (4.14) 

where  the  acoustic  wavenumber  is  defined  in  cartesian  coordinates  as  k  =  (kr,  fc^)  = 
{kx,  ky,  kz)  where  ufc—  ^kl  +  k^  +  k^  =  k.  In  free  space  this  would  simplify  to 


94 


(/(rr,rs,u;)  =  — — ^ ^  — — - re'*!^’-  Pm)  (4.15) 

47r|p, -p,|  47rlp, -p^l 

however,  we  will  use  the  more  general  form  given  in  Eq.  (4.14).  Inserting  Eq.  (4.14) 
into  Eq.  (4.13)  yields 

Aj  oo 

^2,m(rr,'u^)  «  J  dZsg{rr,Tm,w)  JJ  w{p^  -  p^)e-'^'-^P‘-Pm) q^{T„  w)(f  p^,  (4.16) 

0  —00 

Again,  the  double  integral  over  p,  analogous  to  a  horizontal  planar  source  array 
with  broadside  directed  down  toward  the  ocean  bottom.  Following  this  analogy  we 
can  define  this  double  integral  as  the  beamformed  output  [129] 

OO 

BmiK  Zs)  =  JJ  w{p^  -  w)(fp^,  (4.17) 

—  OO 

The  beamformed  output  in  Eq.  (4.17)  can  be  evaluated  for  any  horizontal  wavenum¬ 
ber  vector  kr,  however,  only  those  wavenumber  vectors  where  |kr|  <  w/c  lead  to  a 
propagating  acoustic  field.  This  can  be  illustrated  if  we  consider  the  microseismic 
field  generated  by  a  horizontal  plane  of  sources  near  the  sea  surface  as  described  in 
Eq.  (4.17).  In  previous  works  [78]  it  has  been  shown  that  microseisms  generated  at 
the  ocean  surface  are  first  transmitted  downward  to  the  sea  floor  where  they  then 
propagate  along  the  bottom  as  Rayleigh  waves.  This  vertical  propagation  from  the 
sea  surface  to  the  sea  floor  can  be  modeled  to  first  order  using  the  free-space  Green 
function 
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(4.18) 


where  d  represents  the  water  column  depth.  If  |kr|  exceeds  ujc  then  becomes 
imaginary  leading  to  an  exponential  decay  with  depth  d.  For  example  if  |kr|  =  2oj/c 
then  kz  =  iy/^/c  so  that,  for  cj  =  1.3  rad/s,  c  =  1500  m/s  and  an  ocean  depth 
d  =  1000  m,  the  field  is  attenuated  by  ^  0.2  or  13  dB.  Greater  depths  will 

lead  to  even  more  attenuation.  Because  of  this  evanescent  attenuation  we  will  only 
consider  the  acoustic  field  where  |kr|  <  u/c. 

Assuming  the  acoustic  pressure  field  given  in  Eq.  (4.16)  is  temporally  stationary, 
the  power  spectral  density  may  be  written  as  [91] 


(rr,  -  a;')  =  <  P2.m(rr,  w) A*n(rr,  w')  > 

\g  \g 

^  J  J  dZgdz',g{Tr,Tm,w)g*{rr,Tn,w') 

0  0 

X  <  5,„(k,2,,,„).B;(k',2,,„)  >,  (4.19) 

where  <>  represents  the  expectation  and  r„  =  {Pn>K)  —  (^n,  2/n,  ■^a)-  In  the  next 
section  we  derive  expressions  for  the  second  moment  <  Zs^rn)S*{k' ,  Zg^n)  >  of 

the  beamformed  output  and  the  power  spectral  density  due  to  ocean  surface 

waves. 

The  power  spectral  density  5p^,n  in  Eq.  (4.19)  represents  the  field  at  a  receiver 
from  a  single  pair  of  finite  volumes  Vm  and  I4.  In  the  next  section  we  determine  the 
microseismic  field  generated  by  a  hurricane  where  the  surface  wave  field  is  inhomo¬ 
geneous  and  extends  over  a  large  region  hundreds  of  kilometers  across.  To  do  this  we 
divide  the  hurricane  region  into  finite  volumes  Kn,  as  shown  in  Fig.  4.3,  and  then 
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sum  their  contributions  to  find  the  total  power  spectral  density  of  the  received  field 
written  as 


Sp{Tr,Uj)  = 

m,n 


(4.20) 


4.4  Power  Spectral  Density  due  to  Ocean  Surface 
Gravity  Waves 


Taking  the  ocean  surface  wave  velocity  as  the  first-order  field  we  can  determine  the 
source  term  from  Eqs.  (4.11)  where  g(rg,aj)  is  the  Fourier  transform 


1  °° 

9m(r„i.i)  =  ^  j  qm{r„t)e"^dt 


(4.21) 


of 


9m(ra,  t)  —  2^  (Vl,m(rsi  0  '  0)’ 


(4.22) 


Given  Eq.  (4.8)  for  the  first  order  velocity  we  find 


oo  oo 

V2(vi,^(r„t)  •  vi,^(r„t))  =  jj  If 

— OO  — oo 

X  Fi  {ki  ,  K2) 

-f24^(K0^;;,(»C2)e--“V^'-^*F2(Ki,K2) 
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Hurricane  Surface  Waves 


Figure  4-3:  Geometry  of  the  hurricane  wave  field  and  ocean  waveguide  (not  to  scale). 
In  our  model  the  waveguide  environment  may  be  range  dependent  and  in  this  paper 
we  consider  the  example  of  an  upslope  propagation  from  the  deep  North  Atlantic 
to  Florida.  The  range  and  ocean  depth  parameters  R  and  d  are  given  in  Table  4.6. 
The  depth  of  the  receiver  below  the  ocean  bottom  of  162  m  corresponds  to  the  depth 
below  the  earth  surface  of  the  actual  seismometer  in  Florida.  The  compression  wave 
speeds  Cp  are  1500  and  5200  m/s  in  the  water  and  bottom  respectively.  The  shear 
wave  speed  c*  in  the  bottom  is  3000  m/s.  The  densities  p  in  the  water  and  bottom 
are  1.0  and  2.5  g/cm  respectively. 


98 


(4.23) 


^■A*^{Kx)A,n{l^2)e^  ‘e  '^•F2{Kx,K2)\ 
xe-('“+'‘^)^*d2/cirf2/c2, 

where  o-+  =  [o\  +  <72),  o~  =  (<ti  —  02),  =  (ki  +  /C2),  and  k~  =  (/ci  —  K2).  Note 

that  the  Laplacian  operator  (V^)  in  Eq.  (4.23)  leads  to  the  direction  cosine  functions 


^lx^2x  ^ 


K\K2  K\K2 

=  2ki«2[1  -  cos(0i  -  02)][-l  +  cos(0i  -  ^2)],  (4-24) 


and 


F2(Ki,  K2)  =  [(k  f  -  («*  f  -  («J,  f 


1  + 


^li^2g  K,\yK'2y' 


K1K2  K1K2 

=  2«i/C2[-1  +  cos(0i  -  02)][1  +  cos(0i  -  ^2)],  (4.25) 


where  «+  =  (ki*  +  K2x),  «j[  =  («iy  +  K2„),  =  (/cn  -  K2x),  and  Ky  =  (xiy  -  K2y). 

Substituting  Eq.  (4.23)  into  the  definition  of  qmi^a^t)  from  Eq.  (4.22)  and  taking 
the  Fourier  transform  yields, 


00  00 

=  -UI II 

—00  —00 

X  [v4m(«l)>lm(«2)Fi(Ki,  K2)6{u  - 

+>1;;,(/Ci)/1;;j(k2)Fi(ki,  k;2)(5(c*^  +  cT+)e“‘'^^'^* 

+Am{Ki)A*^{K2)F2{Ki,K2)6{uJ  - 
+A*^{Ki)Am{K2)F2{Ki,  K2)6{u  +  af  )e~'^  P‘ 

(4.26) 
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By  inserting  Eq.  (4.26)  into  the  definition  of  the  beamformed  output  Bm(k,  Zg)  from 
Eq.  (4.17)  we  obtain 


•®m(^)  Zg) 


0\0i 


-urn 

“00  -00 

X  |>im(«l)^m(#C2)-P’l(«l,«2)(5(w  -  0+) 

COO  OO  V 

f  j  w{pg- 

OO  — OO  ' 

■h^m(«l)^m(«^2)i^l(«l,K2)<^(t*^  +  (^^l  +  ^2)) 

/  00  00  V 

^00  “OO  ' 

+  ^m(«l)^m(K2)^2(«l,K2)(5(t*^  -  {<^1  “  ^^2)) 

y  00  00  V 

^00  “OO  ' 

+>i;„(«l)>l^(K2)i^2(«l,«2)<5(t^  +  {oi  -  (^2)) 

Coo  00  V  - 

00  —00  ' 


(4.27) 


Note  that  the  wave  height  spectra  Am  and  A„  in  Eq.  (4.27)  can  be  brought  out  of 
the  spatial  integral  since  they  are  constant  over  the  ocean  surface  area  Sm  defined  by 
the  window  function  w.  Integrating  over  p,  then  leads  to 


OO  OO 

^m(k,^,)  =  jj  ^1^2 

—  OO  —00 

X  [^„(«;i)44„(/t2)/^i(«i,«2)^(w  -  cr'^)(iy(lq.)  *  (5(kr  -  /c+)) 

+24;;,(/Ci)24;;,(«;2)^i(Ki,«;2)^(w  +  cr'^)(M^(kr)  *(5(kr  +  «+)) 

+44„(Ki)44;;,(«2)i^2(Ki,K2)'5(a;  -  a“)(iy(kr)  *  (5(kr  -  k“)) 
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+  a  )(VK(k^)  *  5(kr  +  K  ))] 

xe-('“+'‘^)"*d2^irf2K2,  (4.28) 

where  *  represents  the  two-dimensional  convolution  written  as, 


OO 

F{K)  *  G{K)  =  II  FiK)G{K  -  k;)d2k;  (4.29) 

—  OO 


and  where 


l^(k.)  = 


2sin(^)2sin(^) 


(4.30) 


Note  that,  when  integrating  over  p,,  the  p^  term  in  the  window  function  w{pg  —  p^) 
of  Eq.  (4.27)  introduces  a  phase  factor  Pm.  This  cancels  with  the  other  phase 
term  in  Eq.  (4.27).  From  this  we  see  that  the  source  beamformed  output 

given  in  Eq.  (4.28)  has  no  phase  dependence  on  location  p^. 

The  delta  functions  in  Eq.  (4.28)  show  that  the  horizontal  component  k^  of  the 
acoustic  wavenumber  vector  is  generated  by  and  equals  either  the  sum  or  difference 
K~  of  the  surface  wavenumber.  In  Section  4.3  we  discussed  how  the  non-propagating 
field  where  |kr|  >  ^  is  not  significant  and  may  be  ignored.  Therefore,  when  integrating 
over  K2  we  need  only  consider  a  range  of  integration  f2  corresponding  to  propagating 
acoustic  waves.  For  the  sum  terms  in  Eq.  (4.28)  corresponds  to  the  integration 
region  where  |k'''|  <  ^  and  for  the  difference  terms  Q.  corresponds  to  |k;~|  < 

Here  we  will  define  the  source  area  dimensions  Lx  and  Ly  to  be  ‘acoustically  small’ 
{Lx,Ly  «  Aa).  The  acoustic  wavelengths  of  the  microseismic  field  in  water  range 
from  roughly  6  to  1 1  km  requiring  that  the  dimensions  Lx  and  Ly  be  on  the  order  of  1 
km  or  less.  Note  that  this  requirement  on  Lx  and  Ly,  combined  with  the  requirement 
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that  Lj;,Ly  >>  Xg  from  Sec.  4.2,  means  that  Xg  «  Lx,Ly  «  Aq.  If 

the  function  W (k^)  ~  L^Lj,  in  Eq.  (4.30)  within  the  range  of  acoustically  propagating 

wavenumbers  Iq.  <  u/c.  With  this  approximation  Eq  (4.28)  can  be  written  as 


IT  °° 

-oo  n 

x[^m(«;i)>Im(«2)-f’l(Kl,«2)5(‘*^  -  0^^) 
+>i;,(K,)>i;;,(K2)^i(Ki,/c2)<5(ti^  +  <t+) 

+Am{Ki)A*^{K2)F2{Ku  K2)S{uJ  -  CT“) 
+>i;;,(«i)yl,„(K2)F2(Ki,  /C2)5(a;  +  a")] 

(4.31) 

The  power  spectral  density  of  the  pressure  in  Eq.  (4.19)  contains  the  second 
moment  or  variance  of  the  beamformed  output  BmO^,  Zg)  which  from  Eq.  (4.31)  can 
be  expressed  as, 


*2,  T  j  '2, 

— oo  h  -oo  fi 

X  <  [>1„(«;i)/1,„(k2)Fi(ki,K2)5(w  -  {ai  +  ^2)) 
+>i;;,(«;i  )>!;;;,  (/C2)Fi(ki,«;2)^(u;  +  {oi  +(72)) 
+24,„(/ci)^;;,(«2)F2(«;i,K2)<^(w  -  ((Ti  -  (72)) 

+yi;,(Ki)/l,„(«2)F2(Kl,K2)5(w  +  ((7i  -  (72))] 

X  [>1;(k;3)>1;(k4)Fi(k3,  «;4)5(a;  -  ((73  +  (74)) 

+  >In(«3)>In(K4)Fi(K3,  K;4)(5(w  +  {as  +  (74)) 
+/1;(k;3).4„(k4)F2(/C3,«4)<5(w  -  ((73  “  (T4)) 
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+A„{K3)A*^{k4)F2{K3,  K4)6{uJ  +  ((73  -  CT4))]  > 

(4  32) 

Since  the  wave  height  spectra  Am  and  A^  are  uncorrelated  zero-mean  Gaussian 
random  variables  as  discussed  in  Sec.  4.2,  we  can  substitute  Eq.  (4.6)  into  Eq.  (4.32) 
and  integrate  over  kz  and  /C4  to  obtain 

n  rr 

<  Bm{KZs)B^(k',z'J  >  =  - 32^  i/  jj  (^l<^lSmn[t^\)Smn{t^2) 

-00  n 

X  (kj,  /C2)<5(n;  —  a'^)6{u}'  —  a^) 

+Fi{ki,  K2)6{u  +  a'^)5(uj'  +  o^) 

+F2{ki,K2)6{u)  —  a~)5{uj'  —  a~) 

-f  F2  («!,  K2)^(w  -  o~)5{u'  -j-  a~) 

-I-F|(«i,  K2)<5'(w  +  a~)6{u'  -  a~) 

-f F2^(ki,  K2)S{u  +  a~)S{u)'  -H  a~) 

(4  33) 

Then  by  integrating  K2  in  Eq.  (4.33)  over  the  integration  region  f2,  corresponding  to 
<  c  I**"”!  —  c’ 


<Bm{k,z,)B:,{k',z,)>  =  2p^nLlLie  11  ay, 

—00 

X  [5,nn(«i)5mn(-«i)(^(‘^  -  2ai)5{u'  -  2cti) 
+‘S'mn(^l)5mn(  — «^l)<^(t^  +  2ai)6{uj'  +  2ai) 

+45,„„(«i)5,„„(«i)(5(u;)(5(u;')] 

xe-2'‘>{^'+^«)d2Ki.  (4.34) 
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Since  ^  <<  k  this  integration  leads  to  the  approximations  sa  — K2  in  the  sum 
terms  and  ki  «  /C2  in  the  difference  terms. 

The  term  4S',„„(/ci)5TOn(Ki)^(cL;)^(a;')  of  Eq.  (4.34)  leads  to  a  constant  {u  =  u'  = 
0)  value  in  the  power  spectral  density  that  is  irrelevant.  Also,  for  the  remainder  of 
the  paper,  we  suppress  the  subscripts  of  /ci  and  to  simplify  notation. 

The  wave  height  cross  spectral  density  can  be  converted  to  the  frequency- 

angle  {(7-6)  domain  Smn{(7,6)  by  evaluating  the  Jacobian  where  =  qk  =  p|/c| 
[127].  This  yields  the  relation  5mn(o’, 0)  =  |^/c5mn(/c)  =  ^Smn{K).  Making  this 
substitution  in  Eq.  (4.34)  and  integrating  over  o  yields. 


<Bm{Kz,)B*^{k\z',)> 


l^nLlLlkV  ,,  ,  ^ 

- - 6{uj  -u}) 


512p2 


2ir 

•>/ S™(|.»)Smn(|,«  +  Jr)de(4.35) 


Note  that  the  delta  functions  in  Eq.  (4.34)  introduce  a  frequency  doubling  effect  in 
Eq.  (4.35)  where  the  frequency  u  of  the  acoustic  wave  is  twice  the  frequency  a  of  the 
surface  gravity  wave. 

Substituting  Eq.  (4.35)  into  Eq.  (4.19)  and  integrating  over  u'  leads  to  the  power 
spectral  density  of  the  pressure  field 


'5Pm.n(rr,w)  = 


p^KLlLlk'^u^ 
512p2 


// 

0  0 


dzgdz'g 


Xg{Tr,Tm,w)9*{Tr,Tr„w)e~^  6 
2jr 

X  J  Smn{^,0)Smn{^,6  +  Tr)de 
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(4.36) 


0  0 


Tm,w)g*{Tr,Tr„w) 


where  we  define  the  ‘microseismic  source  cross-spectral  density’  as 


I  <Brn(k,z,)B:{k\z',)>(Lj' 
—00 

p^-KL'^Lyk'^u}^  -uPz,  -w^x'. 


512^^2 


-e  ^9  e 


2n 

X  1 5„„(|,^)5^„(|,0  +  7r)d^. 


(4.37) 


Summing  the  contribution  from  all  finite  volumes  gives  us  the  total  received  power 
spectral  density 


SpiTr,U)) 


f  j 


512^2 


a  “‘-9 


0  0 


xp(rr,r„,w)5*(rr,r„,u;)e  V  e 


27r 


X  /  Smn{^,e)Srani^,e  +  'n)d9 


Xg  Xg 

I  I  dZgdz'^g{Tr,Tm,w)g*{Tr,rn,w) 
m,n  Q  Q 

x5b„„(w,x,). 


(4.38) 


Equation  (4.38)  provides  an  analytic  expression  for  microseisms  generated  by  inho¬ 
mogeneous  ocean  surface  wave  fields,  as  opposed  to  previous  formulations  where  spa¬ 
tial  homogeneity  is  assumed  over  infinite  [51,  16,  54,  77,  62]  or  very  large  [78,  63,  64] 
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surf2ice  wave  areas.  Also  the  Green  functions  g{Tr,Ttn,w)  and  g*{Tr,Tn,w)  in  Eq. 
(4.38)  may  be  calculated  for  any  arbitrary  range-dependent  or  range-independent 
ocean  waveguide  using  standard  propagation  models.  Later  we  will  use  a  formulation 
for  range-dependent  Rayleigh  wave  propagation  to  calculate  the  Green  functions  and 
model  microseismic  propagation  in  a  typical  North  Atlantic  waveguide  environment. 

While  Eq.  (4.38)  is  applicable  to  inhomogeneous  surface  wave  fields  and  arbitrary 
ocean  waveguides,  it  can  also  be  applied  to  more  simple  homogeneous  surface  wave 
fields  and  Green  functions.  In  fact,  in  App.  E  we  find  that,  given  an  infinite  half¬ 
space  and  assuming  the  surface  wave  spectrum  is  range-independent  {Smn{‘^/‘2,6)  = 
S{u/2,d)),  the  power  spectral  density  of  the  pressure  field  in  Eq.  (4.38)  simplifies  to 
the  results  derived  by  Hughes  [54]  and  by  Lloyd  [77]. 

4.5  Microseismic  Source  Levels  in  Hurricane  Bon¬ 
nie 

In  this  section  we  show  how  Eq.  (4.37)  can  be  used  to  find  the  microseismic  source 
cross-spectral  density  generated  by  the  surface  wave  field  in  a  real  hurricane.  In  these 
examples  we  use  the  wave  spectra  from  hurricane  Bonnie.  Through  this  analysis  we 
demonstrate  the  relationship  between  the  wind  speeds,  surface  wave  spectra,  and 
microseismic  source  spectra  in  a  hurricane. 

Beginning  on  22  August  1998,  hurricane  Bonnie  traveled  along  the  east  coast  of 
the  Bahamas,  Florida,  and  South  Carolina  [3].  During  this  time  microseisms  were 
recorded  by  seismometers^  at  the  Disney  Wilderness  Preserve  in  Florida  (28.1°N, 
81.4°W).  Here  we  model  the  microseisms  generated  by  the  complex  surface  wave  field 

'The  IRIS  seismometers  are  part  of  the  Global  Seismic  Network  (GSN)  and  aire  instcilled,  main¬ 
tained  and  operated  by  the  United  States  Geological  Survey  (USGS)  Albuquerque  Seismological 
Laboratory  (http://aslwww.cr.usgs.gov  and  http://www.liss.org).  The  GSN  is  a  cooperative  scien¬ 
tific  facility  operated  jointly  by  the  Incorporated  Research  Institutions  for  Seismology  (IRIS),  the 
USGS,  and  the  National  Science  Foundation  (NSF). 
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of  hurricane  Bonnie. 


Evaluating  the  microseismic  source  cross-spectral  density  Sb^„  of  Eq.  (4.37)  re¬ 
quires  knowledge  of  the  surface  wave  height  spectrum  A,n-  Moon  et  al  [88]  calculate 
the  surface  wave  height  spectrum  for  the  wind  speeds  measured  in  hurricane  Bonnie 
using  the  Wave  Watch  III  (WW3)  wind-wave  modeling  program  [126].  They  also 
show  how  their  results  are  in  close  agreement  with  aircraft  based  measurements  of 
surface  wave  height  spectrum  by  Wright  et  al  [142]  and  Walsh  et  al  [132]. 

Figure  4.2(A)  shows  the  surface  wind  speed  in  hurricane  Bonnie  on  24  August 
1998.  Illustrated  is  the  typical  hurricane  structure  with  the  high  wind  speeds  of  the 
‘eye  wall’  surrounding  the  ‘eye’  at  the  hurricane’s  center.  From  these  wind  speeds, 
WW3  is  used  to  determine  the  wave  height  spectra  ^^(/c)  as  defined  in  Eq.  (4.5) 
(Fig.  4.2).  In  Fig.  4.2  we  see  that  hurricanes  generate  waves  that  propagate  in  many 
different  and  often  opposing  directions.  For  example  in  Fig.  4.2(c)  we  see  a  peak  in 
the  surface  wave  spectrum  (in  red)  with  waves  traveling  North  ((K^./Cy)  «  (0.0, 0.4)), 
West  ((/Ci,«y)  «  (—0.4, 0.0)),  and  South  ((Kx,Ky)  »  (0.0, —0.4)). 

From  the  wave  height  spectra  calculated  using  WW3,  the  microseismic  source 
cross-spectral  density  Sb„„  niay  be  calculated  from  Eq.  (4.37).  We  define  the  ‘mi¬ 
croseismic  source  level’  as 


r  ^  ^  -  in  w  Im'*  ^ 

10  log  (  ^2^2 

X  !/ 


Im^ 
IPa^ls^ 


(4.39) 


expressed  in  dB  re  Pa^s/m'*  where  the  somewhat  arbitrary  length  scales  and  Ly 
have  been  factored  out.  Figure  4.2(B)  shows  7/b^^(u;,  0)  at  w  =  1.26  rad/s  (/  =  0.2 
Hz)  on  24  August  1998,  the  day  where  we  found  the  source  level  to  be  highest.  We 
see  that  the  peak  of  Lb„„  at  location  (c)  is  not  at  the  same  location  as  the  maximum 
wind  speed  (a).  This  is  because,  while  location  (a)  h8is  large  wave  heights  (si  37 
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dB  re  m^/(rad^/m^)),  the  waves  are  all  propagating  primarily  to  the  South  without 
any  opposing  waves.  At  location  (c),  however,  there  are  opposing  waves  (propagating 
both  North  and  South)  which  generate  microseisms  as  expressed  in  Eq.  (4.37).  Also 
the  source  level  is  higher  in  the  low-wind-speed  eye  (b)  than  at  the  maximum  wind 
speed  location  (a)  due  to  the  opposition  in  the  surface  waves  (Fig.  4.2(b)). 

The  arrows  in  Fig.  4.2  represent  the  direction  hurricane  Bonnie  was  moving  and  we 
see  that  there  is  a  small  peak  in  the  source  level  well  ‘behind’  the  hurricane  at  location 
(d)  even  though  the  wind  speed  there  is  relatively  low.  This  is  again  due  to  the 
opposing  waves  at  this  location  (Fig.  4.2(d)).  This  illustrates  the  complex  relationship 
between  wind  speed  and  wave  spectra,  where  the  wave  spectra  is  a  function  not  only 
of  wind  speed  but  also  of  hurricane  geometry  and  translation  speed  [88],  and  the 
complex  relationship  between  wave  spectra  and  microseismic  source  level,  where  the 
level  depends  on  the  nonlinear  interaction  of  opposing  waves. 

4.6  Received  Microseismic  Levels  in  Florida  Due 
to  Hurricane  Bonnie 

In  this  section  we  model  the  propagation  of  the  microseisms  generated  by  hurricane 
Bonnie  through  the  North  Atlantic  waveguide.  The  model  results  are  then  compared 
with  seismic  data  gathered  in  Florida. 

To  determine  the  microseisms  received  at  the  sensor  in  Florida  the  Green  functions 
I’m,  w)  and  g*{Tr,  r„,  w)  of  Eq.  (4.38)  are  calculated  using  the  adiabatic  Rayleigh 
wave  propagation  model  derived  in  App.  F.  The  bathymetry  between  the  hurricane, 
located  over  the  roughly  5  km  deep  Hatteras  Abyssal  plain  on  24  August,  and  the 
seismometer,  located  in  Florida  is  characterized  by  a  gentle  upslope.  We  will  model 
this  up-sloping  environment  with  the  simplified  geometry  shown  in  Fig  4.6.  The 
sound  speeds  and  densities  for  this  model  waveguide  are  based  on  typical  values  for 
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Table  4.2;  Parameters  for  Hurricane  Bonnie  on  24  August  1998 


Time 

Hurricane  Center  Position 
Lat  (°N)  Lon  (°W) 

Range  from 
Sensor  (km) 

Ocean  Depth  at 
Hurricane  Center  (km) 

0000 

24.8 

71.8 

1028 

5.1 

0600 

25.2 

72.1 

983 

5.5 

1200 

25.6 

72.4 

939 

5.5 

1800 

26.1 

72.8 

881 

5.2 

the  deep  North  Atlantic  ocean  [83]  and  bottom  [73]. 

With  this  environment  the  Green  functions  ^(rr,rm,w)  and  g*{Tr,T„,w)  are  cal¬ 
culated  using  Eq.  (F.23).  The  sea-floor  depth  at  the  center  of  hurricane  Bonnie  is 
given  in  Table  4.6,  as  well  as  the  range  between  the  seismometer  and  the  hurricane 
center.  The  seismometer  in  Florida  is  buried  in  the  ground  to  a  depth  of  162  m  as 
shown  in  Fig.  4.3. 

Given  the  source  cross-spectral  density  from  Eq.  (4.37)  and  Green  functions  from 
Eq.  (F.23),  the  power  spectral  density  of  the  field,  from  Eq.  (4.38),  received  by  the 
seismometer  in  Florida  is  calculated.  Figure  4.6  shows  the  power  spectral  density 
of  the  horizontal  velocity  (Fig.  4.6(a-d))  and  vertical  velocity  (Fig.  4.6(e-h))  of  the 
earth’s  crust  at  the  Florida  seismometer  based  on  our  model  (in  red)  and  on  measured 
data  (in  blue). 

The  dashed  blue  lines  in  Fig.  4.6  represent  portions  of  the  measured  data  that  are 
corrupted  by  non-hurricane  related  ambient  noise.  The  level  of  this  ambient  noise  is 
determined  by  measuring  the  noise  levels  during  the  week  before  and  the  week  after 
the  hurricane  passed.  The  un-corrupted  microseism  signals  are  taken  to  be  those  that 
exceed  the  ambient  by  at  least  3  dB  while  anything  below  that  is  considered  to  be 
corrupted. 

The  theoretical  results  and  measured  data  show  reasonable  agreement  with  peaks 
in  the  spectrum  at  roughly  /  =  0.2  Hz  and  a  peak  level  between  50  and  60  dB  re 
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Range  from  hurricane  center  (km) 


Figure  4-4;  The  ocean  depth  (solid  line)  between  hurricane  Bonnie  and  the  seis¬ 
mometer  in  Florida  at  noon  on  24  August.  Also  shown  is  the  oceein  depth  for  the 
idealized  up-sloping  environment  used  to  calculate  the  Green  functions  in  Section  4.6 
(da.shed  line).  The  scale  of  the  figure  makes  the  actual  slope  appear  to  change  rapidly, 
however,  the  maximum  slope  of  the  ocean  floor  is  roughly  3°. 


no 


Frequency  (Hz)  Frequency  (Hz)  Frequency  (Hz)  Frequency  (Hz) 


Figure  4-5:  Horizontal  (a-d)  and  vertical  (e-h)  velocity  spectra,  as  a  function  of 
seismic  frequency  (/  =  2Tru),  modeled  (red)  and  measured  (blue)  at  the  seismometer 
in  Florida  at  four  times  on  24  August  1998.  The  dashed  blue  lines  represent  portions 
of  the  data  that  where  corrupted  by  non-hurricane  related  ambient  noise.  Note  that 
the  peak  in  the  spectra  is  at  roughly  /  =  0.2  Hz.  This  frequency  corresponds  to 
the  peak  in  the  wave-height  power  spectral  level  at  /t  —  0.04  rad/m  (surface  wave 
frequency  /  =  27r cr  =  0.1  Hz)  seen  in  Fig.  4.2  and  in  plots  (i-1).  The  surface  wave- 
height  power  spectral  level  (101ogS'mm(«^)/(^^/(»'od^/^^)))  in  plots  (i-1)  is  taken  at 
the  center  of  the  storm  and  averaged  over  wave  propagation  direction. 
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(nm/s)^/Hz.  Note  that  the  peak  of  the  wave  height  spectrum  in  Fig.  4.2  is  at  roughly 
K  —  0.04  rad/m  which  equates  to  a  surface  wave  frequency  a  ==  0.63  rad/s.  From  the 
frequency  doubling  effect  discussed  in  Sec  4.4,  this  gives  a  peak  in  the  acoustic  field 
at  =  1.26  rad/s  or  f  —  0.2  Hz  which  we  see  in  the  received  field  in  Fig.  4.6. 

4.7  Conclusions 

Here  we  present  an  analytic  model,  based  on  the  non-linear  wave  equation,  to  de¬ 
scribe  the  microseisms  generated  by  a  hurricane.  This  model  is  ideal  for  hurricane 
generated  microseisms  since  it  can  be  used  to  calculate  the  acoustic  field  due  to  spa¬ 
tially  inhomogeneous  surface  waves.  Also,  this  model  may  be  used  in  range-dependent 
waveguide  environments  as  is  the  case  when  a  hurricane  at  sea  generates  microseisms 
that  propagate  up  the  continental  margin  to  a  receiver  on  land.  This  modeling  is 
useful  because  microseisms  are  a  primary  cause  of  noise  in  seismic  measurements 
[74,  107,  105]  that  raise  the  detection  threshold  for  monitoring  earthquakes  [135]  and 
tsunamis. 

Based  on  the  ocean  surface  directional  wave  spectrum  in  hurricane  Bonnie[l42,  88], 
we  hindcast  the  microseismic  source  levels  generated  by  the  nonlinear  interaction  of 
the  ocean  surface  waves.  We  then  model  propagation  of  the  microseismic  field  from 
hurricane  Bonnie  in  the  North  Atlantic  to  a  seismometer  in  Florida  to  hindcast  the 
measured  signal.  We  find  that  these  results  compare  reasonably  well  with  actual 
seismic  measurements. 
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Chapter  5 


Conclusion 


We  have  shown  that  the  wind-generated  noise  in  the  10  Hz  to  10  kHz  range  received 
by  a  single  underwater  acoustic  sensor  in  a  hurricane  can  be  well  approximated  by 
sea-surface  contributions  so  local  that  wind  speed  and  surface  source  intensity  can  be 
taken  as  nearly  constant.  Two  terms  with  empirically  and  analytically  determined  de¬ 
pendencies  may  be  used  to  estimate  wind  speed  from  measured  ambient  noise  spectral 
level  (1)  a  universal  ambient  noise  source  term  and  (2)  a  local  waveguide  calibration 
term.  At  low  frequencies,  current  evidence,  including  measurement  of  the  underwater 
noise  from  hurricane  Gert,  demonstrate  a  simple  power-law  relationship  exists  be¬ 
tween  noise  intensity  and  wind  speed  so  that  the  log  of  wind  speed  may  be  estimated 
accurately  from  spectral  ambient  noise  level  by  linear  leeist  square  estimation.  The 
power-law  relationship  and  small  RMSE  between  noise  intensity  and  local  hurricane 
wind  speed  shows  that  estimates  of  hurricane  wind  speed  based  on  underwater  noise 
measurements  could  be  very  accurate  with  an  error  margin  of  only  5%.  At  higher 
frequencies,  a  non-linear  relationship  is  expected  but  we  show  that  it  should  be  possi¬ 
ble  to  make  unambiguous  low-variance  wind  speed  estimates  from  broad-band  noise 
measurements.  An  underwater  acoustic  sensor  could  then  effectively  act  as  an  acous¬ 
tic  anemometer,  providing  an  accurate  estimate  of  a  hurricane’s  wind  speeds  and 
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destructive  power.  Currently  many  ocean  acoustic  systems  exist,  such  as  the  PMEL 
sensor  used  in  this  work,  that  could  be  used  for  meteorological  measurements  and 
additional  systems  could  be  deployed  from  ships,  aircraft,  or  near  shore  at  relatively 
low  cost. 

We  also  develop  a  new  theoretical  model  for  the  0.1  to  1  Hz  microseisms  generated 
by  the  nonlinear  interaction  of  sea  surface  waves.  This  model  is  suitable  for  hurricane 
generated  microseisms  since  it  accounts  for  surface  wave  fields  that  change  with  posi¬ 
tion  and  for  complex  and  potentially  range-dependent  environments.  The  microseisms 
generated  by  hurricane  Bonnie  and  received  by  a  seismometer  in  Florida  is  calculated 
and  shown  to  agree  favorably  with  measured  data.  This  analysis  may  be  useful  be¬ 
cause  microseisms  are  a  primary  cause  of  noise  in  seismic  measurements [74,  107,  105] 
of  earthquakes[135]  and  tsunamis.  It  may  also  be  useful  for  determining  the  isotropy 
of  the  microseismic  field  necessary  for  surface  wave  tomography  [117]. 
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Appendix  A 


Noise  Correlation  from  Random 
Surface  Sources 


Previous  models  for  sea-surface  noise[70, 95, 19]  contain  approximations  or  parameter- 
izations  that  are  not  valid  for  the  rapidly  spatially- varying  source  levels  of  a  hurricane, 
particularly  in  the  case  where  the  hydrophone  is  near  or  under  the  hurricane.  Because 
of  this  an  alternative  expression  for  calculating  the  spatial  cross-spectral  density  of 
the  noise  field  is  necessary.  The  geometry  for  modeling  the  spatial  cross-spectral  den¬ 
sity  from  uncorrelated  noise  sources  at  the  surface  of  a  stratified  ocean  waveguide  is 
shown  in  Fig.  2-2. 

The  pressure  field  created  by  each  surface  source  is  given  by  the  solution  to  the 
Helmholtz  equation 


0(r,/)  =9(ro,/)p(r|ro,/)  (A.l) 

where  </)(r,  /)  is  the  pressure  at  r  given  a  source  spectral  amplitude  ^(ro,  /)  at  Tq  and 
p(r|ro,  /)  is  the  waveguide  Green  function.  The  total  noise  field  0s(r,  /)  is  the  sum 
of  the  fields  radiated  from  each  source 
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M  M 

<^5(r,  f)  =  Y.  /)  =  /)^(r|rm,  /)  (A.2) 

m  m 

The  spatial  correlation  of  the  total  field  between  two  receivers  ri  and  r2  can  then 
be  written  as 


A!(ri,r2, /',/")  =  <MriJ'mr2,n> 

=  Z  L  <  f')Q*K,  f")  > 

m  n 

xy(ri|rm,/')p*(r2lr„,/")  {A.3) 

If  the  sources  have  zero  mean  and  are  uncorrelated,  then  <  q{Tm,  f')Q*{Tn,  f")  >= 
^m,n  <  9(rm) /09*(*'n) /")  >,  where  Sm,n  is  the  Kronecker  delta  function  and  the 
correlation  simplifies  to 


fi(ri,r2, /',/")  =  E  <?(■■-. /V(rm.r)> 

m 

X5(ri|rm,/')^*(r2|rm,/")  (A.4) 

Assuming  that  the  source  function  q  and  the  received  field  (ps  can  be  taken 
to  follow  a  stationary  random  process  in  time,  at  least  over  the  measurement  pe¬ 
riod,  then[92]  <  (t)s{ri,  f)(l)*s{^2,  f)  >=  C{ti,T2,  f)S{f  -  /')  and  <  q{f)q*{f')  >  = 
Sqq{f)^{f  —  f)  where  5,,(/)  is  the  power-spectral  density  of  q  and 


C{Ti,T2,f)  =  Y^Sgg{Tm,f)g{Tl\T^,f)g*{T2\TmJ)  (A.5) 

m 

is  the  cross-spectral  density  of  (psi^i)  and  0s (r2). 
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For  dense  source  concentrations,  this  sum  can  be  expressed  as  an  area  integral 


C{Ti,r2j)  =  j  d^Po^^^^^^9{ri\To,f)g*{T2\ToJ)  (A.6) 

where  AA  is  equal  to  or  greater  than  the  coherence  area  of  the  random  source  distri¬ 
bution  and  corresponds  to  the  smallest  differential  area  summable. 

The  variance  of  the  source  amplitude  equals  the  integrated  power  spectral  density 
Sqq{f)df.  Since  this  variance  is  asymptotically  approximated  by  the 
sample  variance  <  \q{t)\'^  >~  ^  Jq  \q{t)\^dt  for  large  measurement  windows  T  and 
since  Parseval’s  Theorem  has  ^  Jq  \q{t)\'^dt  =  ^  \Q{f)\^df,  we  may  deduce  that 

s„(f)  «  (A.7) 

which  relates  Sgq  to  practical  measurements. 

Equation  (A.6)  is  similar  to  an  intermediate  expression  derived  by  Perkins,  Ku- 
perman,  Ingenito,  Fialkowski  and  Glattetre  (Eq.  (6)  of  Ref.  [95])  derived  for  a 
significantly  different  physical  scenario  as  noted  in  section  2.3.  Equation  (A.6)  can 
be  used  to  model  the  spatial  correlation  of  the  noise  field  for  uncorrelated  surface 
generated  noise  when  the  source  distribution  or  waveguide  is  range  dependent.  When 
the  source  distribution  and  environment  are  range  independent,  Eq.  (A.6)  reduces  to 


C{Ti,T2,!)  =  I  Mfc, 

Xg{kr,Zi,Zo)g*{kr,Z2,Zo)Jo{kr\pi  -  p2\)  (A.8) 

following  the  Kuperman  and  Ingenito  approach[70]  where  g{kj.,  Zi,Zo)  is  the  wavenum¬ 
ber  transform  of  5(ri|ro,  /). 
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Appendix  B 


Asymptotic  Variance  for  Wind 
Speed  Estimates 


The  variance  of  a  wind-speed  estimate  is  evaluated  numerically  in  Section  2.4.3  for 
some  specific  scenarios.  A  general  analytic  expression  is  derived  here  for  the  asymp¬ 
totic  variance  of  the  wind  speed  estimate  for  large  sample  size  //  using  Fisher  informa¬ 
tion.  For  the  expected  intensity  I  of  an  acoustic  measurement  with  a  signal-to-noise 
ratio  or  time-bandwidth  product  //,  the  inverse  Fisher  information  or  Cramer-Rao 
Lower  Bound  (CRLB)  of  a  wind  speed  estimate  V  is  given  as[79,  80] 


2\  -1 


(B.l) 


which  is  the  asymptotic  variance.  [92] 

At  low  frequencies  the  relationship  between  intensity  and  wind  speed  can  be  ex¬ 
pressed  using  Eq.  (2.22),  which,  when  inserted  into  Eq.  (B.l),  yields 
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At  higher  frequencies,  where  attenuation  due  to  bubbles  becomes  important,  the 
relationship  between  intensity  and  wind  speed  follows  Eq.  (2.19),  substituting  Eqs. 
(2.15),  (2.21)  and  (2.23),  so  that 


V'^ 

V O,rasymptotic(^)  —  VL  da  (^•^) 

10io9(e)Sv) 

For  the  attenuation  a{V,  f)  described  by  Weston[139]  in  Eq.  (2.4),  the  CRLB  becomes 


^  aaymptotici^  )  —  N 


_ _ 

p(n~6A6xlO-‘’LV^y/fy 

_ 

/j(n-1.69xlO-8r,V3/)2 


/  <  l.bkHz 
f  >  1.5kHz 


(B.4) 


For  V  =  Vmax  where  Vm^x  is  given  in  Eq.  (2.24),  the  CRLB  goes  to  infinity  indicating 
that  an  unbiased  estimate  of  wind  speed  is  not  possible  for  that  wind  speed  and 
frequency.  This  problem  can  be  overcome  by  broadband  intensity  measurements. 
For  a  wind  speed  estimate  given  intensity  measurements  at  multiple  frequencies,  Eq. 
(B.4)  becomes 


V  ar, 


asymptotic 


(E/  ‘^{n  -  6.46  X  lO-'^LV^y/Jf'^ 
(E/  ^{n  -  1.69  X  10-«LV"V)") 


-1 


-1 


/  <  1.5kHz 

(B.5) 

/  >  1.5kHz 
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which  remains  finite  and  can  be  made  small  by  increasing  the  time-bandwidth  product 

fi. 
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Appendix  C 


Review  of  “Microseism  and 
infrasound  generation  by  cyclones” 
by  Bowen,  Richard,  Mancini, 
Fessatidis  and  Crooker 


It  has  long  been  known  that  microseisms  are  generated  by  the  non-linear  interaction 
of  occean  surface  waves  [78,  51,  16,  54].  In  2003,  however,  Bowen,  Richard,  Mancini, 
Fessatidis  and  Crooker  [14]  published  an  alternative  explanation  for  the  generation 
of  microseisms  by  sea  surface  waves. 

The  theory  they  proposed  is  based  on  the  argument  that  “the  sum  of  the  pressure, 
the  kinetic  energy  per  unit  volume,  and  the  potential  energy  per  unit  volume  [of  the 
surface  gravity  waves]  must  have  the  same  value  at  all  points  in  the  [water]  column, 
even  if  any  of  these  variables  are  varying  in  time.  The  immediate  consequence  of 
this  is  that  the  pressure  time  dependence  at  the  base  of  the  column  of  water  [the  sea 
floor]  is  determined  by  the  time  dependence  of  the  surface  gravity  waves  at  the  top 
of  the  column,  and  this  effect  is  independent  of  depth.”  In  effect  they  propose  that 
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the  potential  field  generated  by  linear  surface  waves  does  not  decay  with  depth  but 
rather  extends,  unattenuated,  to  the  sea  floor.  This  contradicts  classic  linear  wave 
theory  where  the  linear  surface  wave  field  decays  exponentially  with  depth  [66].  The 
purpose  of  this  appendix  is  to  point  out  the  error  in  their  proof  in  Sec.  Ill  that  leads 
to  the  non-decay  of  the  surface  wave  field,  without  commenting  on  the  vadidity  of  the 
rest  of  the  paper 

To  begin,  the  authors  give  the  equations  for  surface  wave  height  in  Eq.  (14) 


rj{x,  y,t)  =  A  sin  {kx  —  ujt) 
and  vertical  velocity  in  Eq.  (15) 


v{x,  y,  t)  =  ojA  cos  {kx  —  cot). 

They  then  seek  to  determine  the  pressure  at  the  sea  floor  using  Bernoulli’s  equation 
which  they  write  in  Eq.  (12)  as 


Pa  +  P9v{x,  y,  t)  -I-  ^/9u(x,  y,  T],  =  P{x,  y,  -h)  -  pgh 

where  Pq  is  atmospheric  pressure,  p  is  water  density,  and  z  =  —  h  is  the  ocean  bottom. 
Combining  these  three  equations  they  derive  the  pressure  at  the  sea  floor  to  be 


P{x,  y,  —h)  =  P^  -I-  pg{h  4-  A  sin  {kx  —  ujt))  +  ^puj^A^  cos  {kx 


JL) 


in  Eq.  (16),  which  seems  to  show  that  there  is  a  fluctuating  pressure  on  the  sea  floor 
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which  does  not  decay  with  depth. 

The  error  in  this  proof  is  in  the  use  of  Bernoulli’s  equation  in  Eq.  (12).  This 
particular  form  of  Bernoulli’s  equation  is  only  applicable  to  steady  state  flows,  not 
to  the  time-dependent  fluid  flow  in  a  surface  wave  [66].  In  other  words  Eq.  (12)  was 
written  incorrectly  and  none  of  the  terms  in  it  should  depend  on  time  t.  Because  of 
this  misuse  of  Bernoulli’s  equation,  the  result  given  in  Eq.  (16)  cannot  be  justifled. 


125 


126 


Appendix  D 


Spectral  Properties  of  Ocean 
Surface  Waves 

D.l  Homogeneous  Surface  Wave  Fields 

In  Eq.  4.1  we  define  the  complex  surface  wave  height  in  the  form  of  an  inverse  Fourier 
transform  such  that 


OO 

C(p, 0)  =  jj  A{K)e*^'P(fK  (D.l) 

— OO 

The  surface  wave  height  spectra  A{k)  can  then  be  written  as  the  Fourier  transform 


1 

^  Jj  C{p,0)e-''^  P(fp  (D.2) 

— OO 

We  assume  that  surface  wave  height  has  zero  mean  (<  C(P)O)  >=  0)-  Taking  the 
expectation  of  Eq.  (D.2) 
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(D.3) 


1  °° 

<A(k)  >= -^  jj  <  C(p,0)  > 

— oo 

we  can  see  that  the  wave  height  spectra  must  also  have  zero  mean  (<  A{k)  >—  0). 
The  second  moment  <  A{k)A*{k')  >  can  also  be  derived  using  Eq.  D.2  such  that 


<  >1(k)>4*(/c')  > 


oo  1  °° 

<  3^  II  II  C(p',0)e<-'l^cPp'  > 

— OO  -oo 

00  oo 

T^JI  II  <C(P,0)C(p',0)>e-’«-^e‘'"'-^'dW(D.4) 

— OO  — oo 


If  we  define  p"  =  p  —  p'  we  can  write  this  as 


00  00 

<  A(k)A‘{k')  >  =  II  II 

— OO  — oo 

oo  oo 

=  1^#  fi(p")e-'''^'dV"j//  e-‘''  Pe-'‘'  P<i^p 

— oo  — oo 

=  S{k)S{k  —  k')  (D-h) 

where  R{p")  =<  C(p.o)C(p'.o)  >  is  the  wave  height  correlation  and  the  power 
spectral  density  S{k)  is  its  Fourier  transform.  This  expression  shows  that,  for  homo¬ 
geneous  surface  wave  fields,  the  different  wavenumber  components  are  uncorrelated. 

In  Sec.  4.2  we  also  assume,  based  on  the  Central  Limit  Theorem,  that  the  complex 
wave  height  C(P)^)  is  a  Gaussian  random  variable.  Since  the  spectrum  A{k)  is  the 
Fourier  transform  of  the  complex  wave  height  (Eq.  (D.2),  it  must  also  be  a  Gaussian 
random  variable.  The  fourth  moment  of  a  Gaussian  random  variable  can  be  written 
in  terms  of  its  second  moments  such  that 
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<  v4(Ki)>1(/C2)-4*(/C3)>1*(k4)  >  = 

5(Ki  -  K3)S{ki)6{K2  -  K4)S{k2) 

+5{Ki  -  K4)S{Ki)S{k2  -  K3)S{k2).  (D.6) 

We  can  also  examine  the  cross  correlation  between  two  surface  wave  processes  Cm 
and  Cn  where 


00 

Cm(p,  =  jj  (D.7) 

— oo 

end 


CXJ 

4(p.  0)  =  jj  A„{K)e''^  f’ipK 


(D.8) 


The  surface  wave  height  spectra  can  then  be  written  as  the  Fourier  transform 


OO 

^m(«)  ^  jj  Cm(p,0)e“’''-^dV 


(D.9) 


and 


1  °° 

A„(k)  =  -^JJ  Lip,0)e~‘^'P(fp 


(D.IO) 


Again  we  assume  that  surface  wave  heights  and  the  wave  height  spectra  have  zero 
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mean  so  that  the  second  moment  <  AmiK)A*„{K')  >  can  be  written  as 


OO  1  °° 

^  II  II  C(p',0)e''''-^rfV 

-OO  — OO 

OO  OO 

=  II  II  <  Cm(p,0)c;;(p',0)  > 

—  OO  — OO 
OO  OO 

=  ik-JI  II 

—  OO  — OO 

OO  OO 

=  ib  II  II 


> 


=  <5„n(«)<5(»c  -  k') 


(D.ll) 


where  Rmnip")  =<  Cm(p, 0)C^(p', 0)  >  is  the  wave  height  cross  correlation  and  the 
cross  spectral  density  5nin(K)  is  its  Fourier  transform. 

Again  we  assume  that  the  spectra  and  A„(ic)  are  Gaussian  random  vari¬ 

ables  so  that  the  fourth  moment  can  be  written  as 


A,7,(/Ci)A,7i(/C2)A^(k3)A„(K4)  >  — 


(5(/Ci  -  K3)Smn{Ki)6{K2  -  K4)Smn{l^2) 


-|-(5(Ki  -  K4)Smn{Kl)S{K2  -  #t3)‘5,„n(«2)- 


(D.12) 


D.2  Inhomogeneous  Surface  Wave  Fields 

Now  we  consider  the  case  where  the  surface  wave  height  is  defined  over  a  finite  area. 
The  complex  surface  wave  height  in  a  finite  region  centered  at  the  origin  can  be 
written  as 
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4(p.o)  =  i(p,t)w(p) 


(D.13) 


where  the  window  function  w  is  unity  for  —Lx/2  <  x  <  Lx/2  and  —Ly/2  <  y  <  Ly/2 
and  zero  elsewhere.  Taking  the  Fourier  transform  (Eq.  (D.2))  of  both  sides  of  Eq. 
(D.13)  yields 


OO 

A^{k)  =  A{k)  *  W{k)  =  jj  A{k')Wo{k  -  k')Sk'  (D.14) 

— OO 


The  second  moment  of  A^  can  then  be  written  as 


OO  OO 

IJ  W{k-  K')cfK'  II  W{k" 

— OO  -OO 


X  <  A{k')A*{k"')  > 


K'")d^K"' 

(D.15) 


which,  from  Eq.  (D.5)  becomes 


<  A^{k)AI,{k")  > 


OO  OO 

II  W{K-K')cfK’  II  W*{k"  -  K"')d^K"'S{K')S{K'  -  k'") 

—OO  —OO 

OO 

II  Win- k’)W*{k"  -  K')S{K')d'^K'  (D.16) 

—  OO 


As  the  dimensions  Lx  and  Ly  become  much  larger  than  a  wavelength  Xg  the 
functions  W  begin  to  approximate  delta  functions  such  that 
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(D.17) 


<  A^{k)A*^[k")  >  TH  j  S{k  —  k')5{k"  —  K')S{K')(fK' 

«  6{k  —  k'')S{k)  —<  yl(K)>l*(K")  > 


We  can  now  consider  the  more  complicated  case  of  the  cross  power  spectral  density 
between  the  surface  wave  heights  in  two  finte  regions  m  and  n. 


(w,m{Py  0)  =  C(P,  0)w(p  -  p^) 


(D.18) 


and 


Cu,,n(p,0)  =  C(p,0)w(p-p„)  (D.19) 

which,  after  taking  the  Fourier  transform,  become 


*  (W(K)e  ’'^  P-n) 


(D.20) 


and 


=  ^„(k)  *  {W{K)e-^^^  Pn)  (D.21) 

As  before,  we  can  take  the  second  moment  which  becomes 
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00 
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Again,  as  the  dimensions  Lx  and  Ly  become  much  larger  than  a  wavelength  A<, 
the  functions  W  begin  to  approximate  delta  functions  such  that 


<  >«  -  k")S„„(k)  (D.23) 

This  expression  is  similar  to  the  previous  Eqs.  (D.5)  and  (D.17)  with  the  addition 
of  an  exponential  phase  term.  Since  the  delta  function  is  zero  for  all  k'  —  k'  ^  0  and 
since  the  phase  term  is  unity  when  k'  —  k'  =  0  this  expression  can  be  reduced  to 


<  >«  s{k 


=<  A^(/c)A;(k")  > 


(D.24) 
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Appendix  E 


Microseisms  in  a 

Range- independent  half- space 


To  compare  the  expression  in  Eq.  (4.38)  with  derivations  by  others  we  can  simplify 
our  solution  for  the  case  of  a  range-independent  surface  wave  field  over  an  infinite 
ocean  half-space.  The  Green  function  for  a  source  near  the  free  surface  of  an  infinite 
ocean  half-space  can  be  written  as  a  dipole 


^ifcRm 

27rRm 

where  Rm  =  |rr  —  r^|  is  the  distance  between  source  and  receiver  positions  and 
Fr  respectively  and  Pm  =  Atan(|p^  —  p^\/\zr  —  Zm\)  is  the  angle  from  vertical.  Prom 
Eq.  (4.38)  we  write  the  power  spectral  density  of  the  pressure  field,  substituting  Eq. 
(E.l)  as 


9i^ri  w) 


kZgCOs{^Pm)  1 


(E.l) 


Sp{Tr.,U}) 
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X  j  j  dzsdz[e  ^9  ’  e 


0  0 


^  X!  ;r-^^^sCo«{/?m)  v-5-fc^icos(^n)  (E-2) 

27rilm  27r/t„ 


For  densely  spaced  and  r„  the  sums  can  be  approximated  as  integrals  on  and 

p's 


Sp{Tr,ljj) 


2  2  / 


^  j  J  dz^dz'.e  29  °  e  ^9 
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oo  oc 

>>// 


“OO  “OO 

oo  oo 


cPp.  e 


LjLy  27ri? 


.ikR 

^  kzsCos{P) 


r  f  dp o*  6 

/  /  z:t; 

—  /V*!  ^/■V>  “ 


(E.3) 


For  /czr  >>  1  the  integral  over  can  be  approximated  as  ^.)+»’r/2 

integral  over  p,  as  so  that 


27r 


Sp{Tr,u) 


P^Tlk'^U}^  f 
512^2 


I  s(^,e)s(^,e  +  T,)de 


^  J  J  dzsdz'^zsz'ge  e  e'''^’e 
0  0 


(E.4) 


Integrating  over  Zg  then  yields 


Sp{Tr.,u) 


p^'nk'^ug^ 

32 


27r 


f  s('^,e)s{‘^,e  +  -w)de 


(E.5) 
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since  ^  >>  A:  at  the  frequencies  considered  here.  This  result  is  identical  to  Eq.  (33) 

0^3 

of  Hughes  [54]  and  Eq.  (35)  of  Lloyd  [77]  where  S{o,9)  equates  to  -^X{ki)G{9)  from 
Hughes  and  to  ^/(«)  from  Lloyd  (including  the  ^  correction  of  Hughes  derivation 
discussed  by  Lloyd). 
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Appendix  F 


Adiabatic  Propagation  of 
Generalized  Rayleigh  Waves  in  a 
Range  Dependent  Ocean 
Environment 


To  determine  the  microseisms  received  by  a  sensor  either  in  the  ocean  or  in  the  earth’s 
crust,  the  Green  function  of  Eq.  (4.38)  must  be  calculated.  It  has  been  shown  both 
theoretically  [78,  105]  and  experimentally  [47,  75,  45,  15,  11,  46]  that  microseisms 
propagate  as  Rayleigh  waves  along  the  sea  floor. 

[122]  shows  how  the  depth  or  thickness  of  the  ocean  layer  affects  Rayleigh  wave 
phase  speed  and  how,  for  large  depths  or  frequencies,  there  may  be  multiple  propa¬ 
gating  Rayleigh  modes.  [105,  106]  and  [35]  later  show  that  the  Rayleigh  wave  field 
can  be  expressed  as  a  sum  of  normal  mode  contributions.  Unlike  ‘cl8issic’  Rayleigh 
waves  which  are  defined  to  propagate  along  a  vacuum/elastic  boundary,  [39]  coin  the 
term  ‘generalized  Rayleigh  waves’  to  describe  Rayleigh  waves  that  propagate  along 
an  elastic  boundary,  such  as  at  the  sea  floor,  under  a  finite  thickness  fluid  layer  like 
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the  ocean.  The  modal  expressions  of  [105,  106]  and  [35]  are  only  applicable,  how¬ 
ever,  to  environments  where  the  ocean  depth  is  constant.  Varying  water  depth  at  the 
continental  margin  can  affect  the  Rayleigh  wave  propagation. 

To  solve  the  range-dependent  propagation  problem,  [56]  divide  the  waveguide  into 
range-independent  segments,  each  with  its  own  numerically  calculated  eigenvalues  and 
eigenfunctions,  and  apply  the  adiabatic  normal  mode  approximation  to  evaluate  the 
seismic  and  acoustic  propagation  between  segments.  The  computational  nature  of 
their  approach,  however,  does  not  provide  the  physical  insight  given  by  the  anaJytic 
results  of  the  earlier  range-independent  studies  [122, 105,  106,  39,  35].  Here  we  provide 
an  analytic  solution  for  the  range-dependent  case  where  a  fluid  ocean  of  varying  depth 
overlays  an  elastic  ocean  floor,  as  opposed  to  the  computational  solution  method  given 
by  Arvello  and  Uberall  [56]. 

We  consider  the  simple  case  where  a  monopole  source  is  located  in  a  homogeneous 
fluid  medium  overlaying  an  elastic  half  space.  Unlike  previous  range-independent 
studies,  here  the  depth  or  thickness  of  the  fluid  layer  may  change  with  range  from 
the  source.  For  the  example  in  this  paper  we  model  an  up-sloping  environment  where 
depth  decreases  with  range,  however,  other  depth  profiles  may  also  be  considered. 

It  should  be  noted  that,  in  addition  to  modeling  microseisms,  this  range-dependent 
Rayleigh  wave  model  also  has  application  in  other  seismic  research  where  a  sensor  on 
land  measures  signals  generated  by  sources  at  sea.  For  example  undersea  earthquakes 
[124]  and  underwater  explosions,  like  those  monitored  under  the  Comprehensive  Nu¬ 
clear  Test  Ban  Treaty  [87,  110],  also  produce  Rayleigh  waves  on  the  sea  floor.  The 
propagation  of  these  Rayleigh  waves  is  also  often  measured  to  infer  the  geologic  char¬ 
acteristics  of  the  ocean  floor  using  surface- wave  tomography  [8,  103,  119]. 

We  begin  by  considering  the  range-independent  form  of  the  normal  mode  solution 
for  a  source  at  a  depth  zq  in  a  constant-depth  fluid  layer  overlaying  an  elastic  halfspace 
|56] 
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(F.l) 


n 


=  Y.'^^A^s)u2,n{zT)H'^\knp), 

n 


(F.2) 


and 


^2{P,Z)  =  Y,Ui^n{Zs)v2,n{Zr)Ho\knp),  (F.3) 

n 

where  p=  |p|  is  the  horizontcd  distance  from  source  to  receiver  and  Zr  and  z,  are  the 
receiver  and  source  depth  respectivelly.  The  variable  (^i  represents  the  displacement 
potential  due  to  compressional  waves  in  the  fluid  medium,  and  02  and  02  represent 
the  compression  and  shear  potentials  in  the  solid  medium.  The  variable  k^  is  the 
horizontal  component  of  the  wavenumber  and  tt2,n  and  V2,n  are  the  modal  eigen- 
fuctions  for  each  mode  n.  The  potentials  are  related  to  the  horizontal  displacement 
dr,  vertical  displacement  d^,  vertical  normal  stress  and  vertical  shear  stress  p^r  by 
[106,  39] 


dr  drdz 

,  90  9^0  uP'  , 

a?  + 

o  dvj 

Pzz  =  AV^0  +  2/X— 
oz 

,dq  dw. 


(F.4) 
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The  Lame’s  constants  A  and  ^  are  related  to  the  compressional  wave  speed  in  the 
fluid  ui  and  the  compression  and  shear  wave  speeds  in  the  solid  a2  and  02  by 


=  Ai/pi 

~  (^^2  +  2//2)/P2 

0l  =  /^2/p2,  (F-S) 

where  pi  and  p2  are  the  densities  in  the  fluid  and  elastic  media  respectively  [106,  39]. 

Pierce  [98]  has  shown  that  for  slowly-varying  range-dependent  waveguides  the 
modal  sums  of  Eqs.  (F.1-F.3)  can  be  written  as 


Jo  '^r,{p')dp' 
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and 


i^2{p,z) 
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(F.8) 


142 


based  on  the  adiabatic  mode  approximation.  This  approximation  requires  that  the 
change  in  the  waveguide  environment  as  a  function  of  range  p  is  negligible  over  a 
wavelength  scale  and  that  there  is  no  coupling  between  modes.  In  these  expressions 
both  the  modal  eigenfunctions  ui,„,  U2,n  and  U2,„  and  the  horizontal  wavenumber 
are  allowed  to  change  with  range  p  from  the  source. 

For  an  iso-speed  fluid  layer  of  depth  H{p)  overlaying  an  elastic  half  space  the 
general  solutions  for  the  mode  functions  «2,n  and  can  be  written  in  terms  of 
sinusoids  and  exponentials  [106,  35]  such  that  Eqs.  (F.6-F.8)  become 


4>i{p,z) 
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sin(^„(p)^r) 


(F.9) 


and 


02  (P,  Z) 
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yjknip) 
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where 
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klip)  =  Kip) 
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vM  = 


Kip)  = 


T  =  kM 
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^2  “  Kip) 


clip)  _  . 
,  clip) 
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(F.12) 


represent  the  vertical  components  of  the  wavenumber  vector  for  each  mode  n.  The 
phase  speed  for  each  mode  is  defined  as  Cnip)  =  u)/Kip)  where  u  is  the  frequency  in 
radians/sec. 

The  amplitudes  $i,  $2  and  ^2  of  the  mode  functions  in  Eqs.  (F.9-F.11)  must  be 
chosen  such  that  the  boundary  conditions 


iPzz)f  —  iPzz)s 
iPzr)s  —  0 

(^2)/  —  id'z)s 


>  at  z  =  Hip) 


(F.13) 


are  satisfied  at  the  fluid  solid  interface  z  =  Hip).  The  subscripts  /  and  s  indicate 
that  the  variable  is  to  be  evaluated  in  the  fluid  or  solid  layer  respectively.  These 
boundary  conditions  yield  a  system  of  three  equations  for  the  unknown  amplitudes 
^>1,  $2  and  ^2  whereby  any  two  of  the  unknowns  may  be  solved  in  terms  of  the 
third.  The  third  amplitude  is  then  normalized  to  some  convenient  value.  Later  we 
will  choose  a  normalization  such  that  our  solution  is  consistent  with  that  of  Press 
and  Ewing  [106]  for  constant  bathymetery  Hip)  =  H.  The  solution  to  this  set  of 
equations  only  exists,  however,  if  the  determinant 
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(F.14) 


+  kl)  sini^riH)  (--^2(7?^  +  A:^)  +  2/i2r?n)e“’'""  -2/i2A:^e-^"" 

0  (2C^  + 

^„cos(^„//)  T/„e“’^" 

equals  zero  [106,  39]  where  the  dependence  of  A:„,  77„,  and  H  on  range  p  has 

been  suppressed  to  simplify  notation.  Equation  (F.14)  reduces  to  the  equation 


Pi  ^'‘^n(p) 
P2  /32  ^n(p) 


tan(^„(p)//(p)) 


-  4/!:^(p)r/„(p)C„(p)  -  (2A:^(p) 


(F.15) 


where  the  roots  A:„(p)  correspond  to  the  horizontal  wavenumbers  for  each  mode  n. 
For  the  case  of  constant  bathymetry  H{p)  =  H,  Eq.  (F.15)  is  identical  to  those  given 
in  earlier  works  [105,  106,  39],  however,  for  this  more  general  range-dependent  case 
the  roots  kn{p)  can  now  vary  with  range  p.  Unlike  a  ‘classic’  Rayleigh  wave  which  has 
only  one  root  or  propagating  value,  ‘generalized’  Rayleigh  waves  may  have  multiple 
roots  kn{p)  depending  on  frequency  and  fluid  layer  thickness  as  shown  in  Eq.  (F.15) 
for  the  range-dependent  case  and  as  shown  by  [105]  for  the  range-independent  case. 
Note  that  Eq.  (F.15)  has  no  roots  for  kn{p)  >  u/v\. 

Given  the  roots  A:„(r)  from  Eq.  (F.15),  we  can  express  the  potentials  as 


g-*/o'’fcn(p'W 
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yknKp) 


145 


Mp^z) 


2  =i. 

- 4  / - P  4 

Hip)  V  p 


^i2{kM,H{p)) 


X 


\J^nip) 


(F.17) 
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where  the  solutions  for  the  modal  amplitudes  are 


^i{kn{rho),H{p))  = 
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(F.19) 


MkM,H{p))  = 

^,■^(2  -  <?M/g)kMH{p) 
\/<Hp)/vt-  lD(k„ip),H{p)) 


and 


<S/,(k„{p),Hip)) 


Dik„(p),H(p)) 


(F.21) 


146 


and  where 
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(F.22) 


Equations  (F.16-F.22)  provide  an  analytic  model  for  the  propagation  of  ‘generalized’ 
Rayleigh  waves  in  a  range  dependent  environment.  Note,  however,  that  the  mode 
amplitudes  have  been  normalized  such  that  $i(A:„(p), /f(p)),  ^2{kn{p) ,  H {p))  and 
^i(^n(p))  ^(p))  are  the  same  as  and  ii/\{kn,H)  of  Press  and 

Ewing  [106]  for  the  case  of  constant  bathymetry  H{p)  =  H  (Equation  (38)  of  [106]  is 
missing  a  factor  of  2  which  was  later  corrected  in  Eq.  (4-184)  of  Ewing,  Jardetzky, 
and  Press[39]).  The  Green  function  used  in  this  paper  can  be  written  as 


47r 

g{Tr,T„w)  =  -jP^^i\Pr^Psl^r,Zs) 


(F.23) 


where  the  y  accounts  for  the  normalization  adopted  by  Press  and  Ewing  [106]. 
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Appendix  G 


Review  of  “Emergence  rate  of  the 
time-domain  Green’s  function  from 
the  ambient  noise  cross-correlation 
function”  by  Sabra,  Roux  and 
Kuperman 


In  the  paper  “Emergence  rate  of  the  time-domain  Green’s  function  from  the  ambient 
noise  cross-correlation  function”  [113]  the  authors  attempt  to  show  how  the  waveguide 
Green  function  can  be  determined  from  measurements  of  the  ambient  noise  correlation 
function  using  the  following  equation 


- - -G{Ta,n,T)  +G{ra,Tb,-T)  (G.l) 

dr 

Sabra,  Roux  and  Kuperman  cite  two  of  their  previous  papers  [112,  111]  to  justify  this 
equation.  Reference  [111],  however,  only  considers  the  case  of  homogeneous  noise  in 
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free  space  and  not  the  case  of  typical  surface  generated  noise  in  an  ocean  waveguide. 
Reference  [112]  does  consider  surface  noise  in  a  waveguide  and  demonstrates  that  Eq. 
(G.l)  is  not  true.  Figures  6  and  7  of  Ref.  [112]  show  that  the  right  and  left  hand 
sides  of  Eq.  (G.l),  while  exhibiting  some  similar  modal  structure,  are  different.  It 
should  be  noted  that  in  these  figures  Sabra,  Roux  and  Kuperman  normalized  both 
sides  of  the  equation  such  that  the  maximum  value  is  unity  thus  hiding  significant 
differences  in  amplitude  which  we  will  demonstrate  shortly. 

To  do  this  we  calculate  the  Green  function  and  noise  correlation  in  a  typical 
shallow  water  scenario  and  show  that  Eq.  (G.l)  is  not  valid  in  that  the  right  hand 
side  and  the  left  hand  side  are  not  even  approximately  equal.  To  begin  the  waveguide 
Green  function  can  be  written  in  the  frequency  domain  as  a  sum  of  normal  modes 

^(ru,  n,  =  ^m{Za)n*^{Zb)Hi^Hkm\Pa  -  Pb\)  (G-2) 

m=l 

where  Ta  =  (paiZa)  and  r*  =  {p^„Zb).  For  surface  generated  noise  in  the  ocean,  the 
correlation  can  be  written  as 

'KQ  |u^(2^g)|  y'fni.Za')U^(^Zb) 

2p‘^k^ 

following  the  convention  of  Kuperman  and  Ingenito  [70].  In  this  formulation  — 
Km  +  ioLm  and  Zg  is  a  surface  source  depth  close  to  {zg  <<  A)  the  ocean  surface. 

We  calculate  both  the  Green  function  and  the  noise  correlation  over  the  10  to  2000 
Hz  frequency  range  using  mode  shapes  Um  and  wavenumbers  km  from  the  KrakenC 
normal  mode  code.  In  this  example  the  waveguide  is  110  m  deep  isovelocity  water 
layer  (c=1500  m/s,  p  —  1000  kg/m^,  a  =  6  x  10“®  dB/A)  over  a  halfspace  (c=1800 
m/s,  p  —  1800  kg/m^,  a  =  1  dB/A).  The  receivers  a  and  h  are  at  a  depth  of  50  m 
and  20  km  apart.  Through  Fourier  synthesis  we  can  then  find  the  time-domain  Green 
function  and  correlation. 
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With  the  time-domain  Green  function  and  correlation  we  calculate  the  right  hand 
side  (RHS)  and  left  hand  side  (LHS)  of  Eq.  (G.l)  as  shown  in  Figs.  G-1  and  G-2.  The 
time  series  shown  in  these  figures  start  at  r  =  13.33  sec,  which  is  the  time  required 
for  sound  to  travel  from  point  Ta  to  point  rj.  Figure  G-1  compares  the  RHS  and 
LHS  over  the  first  tenth  of  a  second  in  both  log  and  linear  domain.  First  we  notice 
that  while  the  peaks  of  both  the  RHS  and  LHS  occur  at  the  same  time  they  are  not 
equal  as  suggested  by  Eq.  (G.l).  These  peaks  correspond  to  the  modal  arrivals  with 
the  lower  order  modes  coming  first  («  13.35  sec)  and  the  higher  order  modes  coming 
later.  We  see,  however,  that  early  on  the  peaks  of  the  RHS  are  higher  but  that  they 
decays  with  time.  Figure  G-2  shows  the  comparison  over  a  longer  (2  sec)  time  history 
and  we  see  that  the  RHS  continues  to  decay  while  the  LHS  remains  fairly  constant. 

From  these  figures  we  see  that  the  RHS,  which  depends  on  the  Green  function, 
‘favors’  the  low  order  early  arriving  modes,  while  the  LHS,  which  depends  on  the 
correlation,  favors  the  higher  order  later  modes.  This  difference  can  be  explained  if 
we  compare  Eq.  (G.2)  for  the  Green  function  with  Eq.  (G.3)  for  the  correlation. 
Both  consist  of  a  sum  over  modes  Um{za)  and  Um{zb),  however,  the  correlation  has 
additional  term  which  effectively  weights  the  sum  over  m. 

This  dependence  on  m  can  be  simulated  using  a  normal  mode  code.  Figure  G-3 
for  an  isovelocity  waveguide  shows  that  varies  significantly  with  m  where 

Zs  <<  A.  We  see  that  the  weighting  tends  to  minimize  low  order  modes  in  the 
correlation  and  accentuates  high  order  modes.  This  explains  why  the  LHS  in  Figs. 
G-1  and  G-2  is  lower  than  the  RHS  for  the  early  arriving  low-order  modes  and  higher 
for  the  later  arriving  high-order  modes. 
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Figure  G-1:  The  RHS  and  LHS  in  log  (upper  plot)  and  linear  (lower  plot)  domains 
from  13.33  to  13.4  seconds.  We  see  that  the  peaks  in  the  curves  occur  at  the  same 
times  but  that  the  amplitudes  of  the  peaks  differ  significantly.  We  also  see  that  the 
RHS  peaks  are  higher  for  the  earlier  low-order  modes  while  the  LHS  peaks  are  higher 
for  the  later  high-order  modes.  Note  that  the  RHS  in  the  lower  plot  has  been  scaled 
by  a  factor  of  five  so  that  it  can  be  compared  against  the  LHS. 
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Figure  G-2:  The  RHS  and  LHS  (upper  plot)  over  a  longer  time  duration.  We 
continue  to  see  the  general  trend  where  the  peaks  of  the  RHS  are  higher  for  low-order 
early  arrivals  and  lower  for  later  high-order  arrivals.  This  can  also  be  seen  in  the 
lower  plot  which  shows  the  difference  between  the  RHS  and  LHS. 
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Figure  G-3;  The  weighting  term  as  a  function  of  mode  number  m  at  several 

OCm^m 

frequencies  in  an  isovelocity  waveguide.  The  source  depth  Zg  «  X  typical  of  ocean 
surface  noise. 
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